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BaCeo.85Gdo.i503-x  (BCGO)  has  high  proton  conductivity  but  does  not  exhibit  any 
of  the  necessary  electronic  conductivity  required  for  use  as  a hydrogen  separation 
membrane.  To  address  this  problem,  a two  phase  composite  Pd-BCGO  was  synthesized 
in  our  lab  to  form  a mixed  ionic  electronic  conductor  (MIEC).  Composites  of  35vol%Pd- 
65vol%BCGO  (BCG35Pd)  and  5vol%Pd-95vol%BCGO  (BCG5Pd)  were  synthesized  by 
solid-state  reaction  and  electroless  deposition,  respectively. 

At  temperatures  lower  than  550  "C,  the  hydrogen  flux  of  BCG5Pd  was  higher  than 
that  of  BCG35Pd.  The  hydrogen  flux  of  BCG5Pd  also  increased  with  time  at  650  °C  and 
continued  to  increase  for  20  hours  until  hydrogen  was  shut  off  This  was  an  unusual 
behavior  compared  with  BCG35Pd.  The  second  phase  Pd  did  not  exhibit  such  a trend 
under  isothermal  conditions,  indicating  the  existence  of  a tertiary  phase. 

Analysis  was  performed  using  X-ray  diffraction  (XRD)  and  infrared  (IR) 
spectroscopy.  The  results  showed  that  the  acidic  electroless  bath  decomposed  part  of 
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BCGO  and  that  a phosphorous  impurity  was  absorbed.  BCGO  and  the  phosphorus 
impurity  reacted  at  high  sintering  temperature  to  form  a BaO-rich  compound 
Ba^Ce{PO^)^  and  leave  excess  Ce02.  The  compound  Ba^Ce(PO^)^  was  most  likely  the 
tertiary  phase. 

Single-phase  Ba^Ce{PO^)^  was  synthesized  successfully  at  1 100“C  by  the 
coprecipitation  method.  The  conductivities  of  Ba^Ce(PO^)^  at  600°C  in  dry  and  wet 

4vol%  hydrogen  were  4.04x10’^  and  1.46xlO''*S/cm,  respectively,  less  than  those  of 
BCGO.  The  proton  transference  number  was  around  0.8  at  600“C  in  hydrogen 
concentration  cells.  Therefore,  the  electronic  transference  number  was  about  0.2. 
Ba^Ce(PO^)^  is  a mixed  conductor  with  dominant  proton  conduction.  However,  neither 
its  proton  nor  electron  conductivity  is  high  enough  for  hydrogen  separation.  Composites 
of  Ce02- Ba^Ce(PO^)^  were  also  prepared.  However,  the  electronic  conductivity 
provided  by  the  hopping  between  Ce^"^  and  Ce'*^  was  also  not  high  enough  to  allow 
hydrogen  separation  in  the  composite  CeOi- Ba^Ce(PO^)^ . 

The  results  of  Raman  spectroscopy  and  magic-angle  spinning  nuclear  magnetic 
resonance  (MAS-NMR)  showed  the  local  structures  of  Ba^Ce(PO^)j  under  different 

atmospheres.  Free  hydrogen  can  be  oxidized  forming  protons  and  incorporated  into  the 
lattice.  In  the  presence  of  water  vapor,  bound  hydrogen  is  produced  and  oxidized  forming 
protons.  The  formed  protons  attach  onto  the  O atoms  of  PO4  group  and  are  incorporated 
into  the  lattice.  The  defect  structure  of  Ba^Ce(PO^)^  was  proposed.  The  conductivities 

of  Ba^Ce(PO^)^  as  a function  of  oxygen  partial  pressure  or  water  partial  pressure 
supported  the  proposed  defect  structure. 
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CHAPTER  1 
INTRODUCTION 

1.1  Application  of  Proton  Conducting  Material 

Due  to  dwindling  natural  energy  resources  and  increased  concern  over  the 
environment,  more  attention  is  being  paid  to  clean  and  efficient  fuels  such  as  hydrogen. 
For  example,  the  efficiency  of  an  internal  combustion  engine  is  20%,  but  the  efficiency  of 
a fuel  cell  is  above  33%  when  running  on  pure  hydrogen.  Fuel  cells  produce  no  pollution 
and  the  only  byproduct  is  water  [1].  It  is  expected  that  hydrogen  consumption  will 
increase  dramatically  in  the  future  with  the  wider  application  of  fuel  cells. 

Inexpensive  and  abundant  sources  of  hydrogen  such  as  coal  gasification,  natural 
gas  reforming,  and  off  gas  streams  from  various  process  industries  usually  contain 
hydrogen  mixed  with  other  gases.  The  largest  quantities  of  hydrogen  are  produced  from 
syn  gas.  Hydrogen  separation  membranes  are  thin  materials  through  which  pure  hydrogen 
is  selectively  permeated  from  a hydrocarbon  based  gas  mixture.  They  are  important  to 
enable  a more  widespread  and  economic  use  of  hydrogen. 

There  are  many  kinds  of  hydrogen  separation  membranes.  One  type  is  a composite 
consisting  of  a porous  ceramic  support  and  a Pd  metal  film.  Hydrogen  is  permeated 
exclusively  by  the  metal  film  and  the  ceramic  provides  the  mechanical  strength  but  the  Pd 
cost  is  prohibitively  high  [2-5].  A second  type  is  the  crystalline,  inorganic  molecular 
sieve  membrane  for  hydrogen  separation.  Crystalline  frameworks  have  tunnel  pores  that 
are  capable  of  size  exclusion  separation  [6-9].  The  other  kinds  are  perovskites  or  related 
type  materials  which  are  mixed  conductors.  In  recent  years,  some  materials  with  different 
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structures  have  also  been  developed;  however,  the  hydrogen  flux  is  low  compared  with 
the  oxygen  membrane  and  more  research  is  needed  before  commercialization  is  practical. 

1.2  Proton  Conducting  Materials 

1.2.1  Proton  Conductor 

The  proton  has  no  electron  shell.  It  resides  within  an  associated  anion’s  electron 
shell.  The  small  size  of  protons  makes  it  difficult  to  jump  from  one  site  to  another 
unassisted.  Protons  usually  do  not  occupy  lattice  positions  in  oxides  but  attach  to  oxide 
ions,  forming  OH q group.  Proton  conduction  in  perovskites  was  first  discovered  in 
SrCeOa  [10].  Since  then,  a series  of  doped  perovskite  oxides  ABO3  (A  =Ba,  Sr,  B =Ce, 
Zr)  have  been  studied  [11-26].  Trivalent  dopants  at  B site  are  charge  compensated  by  the 
oxygen  vacancies  Vq  . Protons  replace  oxygen  vacancies  in  the  presence  of  water  vapor. 
They  also  can  be  incorporated  into  lattice  in  the  hydrogen-rich  atmosphere. 

The  proton  transport  usually  has  two  mechanisms  [27-29]:  the  transport  of  protons 
between  relatively  stationary  host  anions,  known  as  the  Grotthuss  or  ffee-proton 
mechanism,  and  the  transport  with  other  species,  termed  the  vehicle  mechanism. 

Proton  conducting  materials  can  be  used  as  an  electrolyte  in  a cell,  whose  driving 
force  is  an  external  voltage.  When  a gas  mixture  containing  hydrogen  is  applied  at  the 
anode,  hydrogen  will  be  oxidized  forming  protons.  These  protons  will  enter  the 
membrane  and  travel  to  the  cathode.  At  the  cathode  they  will  be  reduced  back  to 
hydrogen  molecules.  This  process  results  in  separation  of  hydrogen  from  other  species 
present  at  the  anodic  side  [30-31]. 

Ideal  hydrogen  separation  membrane  materials  should  be  mixed  ionic-electronic 
conductors  (MIEC),  which  not  only  have  high  proton  conductivity  but  also  have  high 
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electronic  conductivity  and  no  oxygen  ionic  conductivity.  In  addition,  the  proton  and 
electronic  conductivities  should  be  nearly  equal  and  sufficiently  high  (>  5x  10'  Scm'  ). 
Materials  must  also  exhibit  high  catalytic  activity  for  the  dissociation  and  recombination 
of  hydrogen  at  gas/solid  interfaces  [32],  In  that  case,  the  hydrogen  separation  process  will 
have  high  hydrogen  permeation  and  will  not  need  any  external  electrical  source. 

1.2.2  BaCeOs  Based  Proton  Conductors 

Barium  cerates  doped  with  lower  valent  cations  on  the  Ce  site  are  high  temperature 
proton  conductors,  which  can  be  used  as  the  electrolytes  in  a fuel  cell  or  hydrogen 
separation  membrane.  BaCCf^g^Gd ^ (BCGO)  exhibits  the  highest  proton  conductivity 

in  the  BaCeOs  family. 

Doping  of  on  the  Ce  site  of  BaCeOs  can  be  written  as 

Gd,0,  Ss£^o^2Gdce  +3V„-  +2Vl . (1-1) 

Kroger-Vink  notation  [33]  is  used  for  point  defects  in  the  above  formula. 

The  proton  incorporation  in  a hydrogen  atmosphere  follows  the  reaction: 

. K1  (1-2) 

The  proton  incorporation  in  wet  atmosphere  follows  the  reaction: 

H^O  + Vo  -F  01  = 20H^ . K2  (1-3) 

Reaction  ( 1 -2)  indicates  that  annealing  in  hydrogen  atmosphere  incorporates 
protons  and  an  equal  number  of  electrons.  Reaction  (1-3)  suggests  that  incorporation  of 
protons  by  annealing  in  water  vapor  will  change  the  ionic  conductivity  by  reducing  the 
number  of  oxygen  vacancies  and  creating  OH^  groups. 
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The  reaction  may  be  split  into  the  formation  of  ion  vacancies  and  the  formation  of 


protons  [34]: 

0^+2/!  =Ko+y02  K3  (1-4) 

H20  + 2h  +20o=  20H  0+^02-  K4  (1-5) 

The  internal  electron-hole  pair  formation  follows  the  reaction: 

null  = e +h  . K5  (1"6) 

The  electroneutrality  condition  for  Gd  doped  barium  cerate  is 

[h]  + 2[V^]  + [0H^  = [Gd,^]^{e].  (1-7) 


For  the  low  water  partial  pressure  condition,  three  different  Brouwer-type  charge 
neutrality  conditions  are  ( e' , Vo),{ Gd'ce , ) and  ( Gd'ce , /;  ) for  low,  intermediate  and 

high  oxygen  partial  pressure  regions,  respectively  [35].  For  the  high  water  partial 
pressure  condition,  three  different  Brouwer-type  charge  neutrality  conditions  are  ( e , 

OH ^ ),  ( Gd'ce , OH^ ) and  ( Gd'c^ , /;  ) for  low,  intermediate  and  high  oxygen  partial 
pressure  regions,  respectively  [35], 

When  the  concentration  of  electronic  charge  carriers  is  neglected,  the  solution  for 
the  concentration  of  proton  defects  is 

[OHo ]{Ph,o ) = Ph,o (^1  + ] / ^2 T’/z.o ) - 1)  [34].  (1-8) 

For  the  low  condition,  [OHq]  corresponds  to  Sievert’s  law: 

{OHo ]{Ph,o ) = Jj^2lOdce]PH,o  ■ [34]  ( 1 -9) 

The  proton  conductivity  of  BaCcQ^^Gd is  theoretically  proportional  to 
and  it  is  a good  proton  conductor  at  high  temperature.  However,  the  electronic 
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conductivity  of  BaCe^^^fid is  not  high  enough  for  hydrogen  separation. 
Krishnamurthy  [36]  addressed  this  problem  by  using  a two-phase  composite  membrane, 
which  consisted  of  BaCe^^fid ^ jjOj  ceramic  material  and  a metallic  component,  Pd.  In 
this  composite  BCGxPd  (x:  the  volume  percent  of  Pd  in  the  composite),  Pd  provided  an 
internal  conduction  path  for  electronic  transport  and  the  ceramic  matrix  facilitated  proton 
transport. 

Figure  1-1  shows  the  total  conductivity  of  BCGO  and  Pd-BCGO  composites. 
Composites  BCG35Pd  and  BCGlSPd  were  synthesized  by  solid  state  reaction.  BCGlPd 
composite  was  synthesized  by  wet  impregnation  process.  BCGSPd  composite  was 
synthesized  by  electroless  deposition  proeess.  In  the  composites  BCGlPd  and  BCGSPd, 
Pd  was  coated  on  the  surface  of  BCGO  powders  to  obtain  a continuous  path  of  Pd  in 
grain  boundaries  of  sintered  BCGO-Pd  composites.  In  the  composite  BCG35Pd,  the 
volume  fraction  of  the  component  Pd  was  0.35  and  there  existed  a continuous  network  of 
Pd  in  the  BCGO.  Under  this  condition,  Pd  was  said  to  have  established  a percolation  path 
in  the  matrix  BCGO  [37].  The  conductivity  of  BCG35Pd  increased  with  decreasing 
temperature  indicative  of  metallic  conduction.  This  also  indicated  that  an  electronic  path 
was  established  in  the  BCG35Pd  composite.  The  bulk  conductivity  of  BCGlPd  or 
BCGSPd  is  lower  than  that  of  BCGO  in  wet  10vol%  hydrogen,  which  is  shown  in  figure 
1-2  (a).  The  grain  boundary  conductivity  decreases  with  increasing  Pd  concentration  in 
the  composite,  which  is  shown  in  figure  1-2  (b).  This  may  be  caused  by  a grain  boundary 
phase  [36]. 

Figure  1-3,  1-4  and  1-5  show  the  hydrogen  flux  through  BCGO,  BCGlPd,  and 
BCGSPd,  respectively,  as  a function  of  time  at  650  °C.  Figure  1-6  shows  the  hydrogen 
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flux  as  a function  of  temperature  through  different  systems.  The  low  hydrogen  flux  of 
BCGO  was  expected  due  to  its  lack  of  electronic  conductivity.  The  hydrogen  flux  of 
BCG  1 Pd  shown  was  significantly  higher  in  contrast  to  the  lower  conductivity  of 
BCG  1 Pd.  This  was  because  the  electronic  conductivity  was  the  limiting  factor  in 
hydrogen  permeation  when  Pd  content  was  low  [36]. 

Hydrogen  flux  of  BCGSPd  increased  with  time  at  650  “C  in  Figure  1-5.  The 
hydrogen  flux  kept  this  increasing  trend  for  20  hours  until  the  feed  gas  was  shut  off. 
Hydrogen  flux  of  BCG5Pd  demonstrated  better  hydrogen  permeability  than  that  of 
BCG35Pd  when  temperature  was  lower  than  550  ”C  as  shown  in  Figure  1-6.  It  was 
unusual  compared  with  other  materials,  whose  hydrogen  flux  became  stable  within 
several  hours.  The  second  phase  Pd  did  not  exhibit  such  a trend  under  isothermal 
condition.  The  BCG5Pd  was  unique  in  its  synthesis  and  some  impurity  could  be 
introduced  into  the  composite  causing  formation  of  a tertiary  phase  at  high  sintering 
temperatures.  It  was  the  tertiary  phase  that  contributed  to  such  an  unusual  trend.  This 
tertiary  phase  was  expected  to  be  a good  mixed  ionic  electronic  conductor  or  at  least  a 
good  proton  conductor. 

To  identify  this  tertiary  phase  from  composite  BCG5Pd  and  synthesize  it  is  the 
main  objective  of  this  research  work.  Conductivity  and  hydrogen  permeability  of  this 
tertiary  phase  material  around  600  °C  needed  to  be  checked  to  assure  whether  or  not  it  is 
a good  proton  conductor  as  expected.  As  shown  in  figure  1-7  [38],  there  is  a gap  for 
proton  conductor  in  the  300-600  °C  range.  This  intermediate  temperature  range  is  very 
important  for  fuel  cell  and  hydrogen  separation.  If  this  tertiary  phase  material  is  a good 
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proton  conductor,  it  will  be  an  important  development  in  fuel  cells  and  hydrogen 
separation  research. 

BCGSPd  powders  synthesized  by  Krishnamurthy  were  characterized  by  X-ray 
diffraction  and  infrared  spectroscopy.  Characterization  results  in  chapter  3 suggest  that 
the  tertiary  phase  is  most  likely  the  barium  cerium  phosphate.  Doped  phosphates  are 
reported  as  proton  conductors;  however,  their  conductivities  are  not  as  high  as  those  of 
perovskites.  Therefore,  barium  cerium  phosphate,  as  the  tertiary  phase,  is  most  likely  a 
proton  conductor,  but  its  conductivity  behavior  is  unknown.  To  explore  this  hypothesis, 
barium  cerium  phosphate  was  synthesized  and  characterized  in  chapter  4.  Corresponding 
defect  structures  were  proposed  based  on  the  characterization  results.  Its  total 
conductivity  and  ionic  transference  number  were  measured  and  discussed  in  chapter  5 
and  6,  respectively.  Its  possible  use  as  hydrogen  separation  membrane  was  attempted  in 
chapter  7.  Finally,  conclusions  about  this  new  proton  conducting  material  were  drawn  in 
chapter  8. 

1.3  Experimental  Techniques 

1.3.1  Synthesis 

Coprecipitation  refers  to  a process  in  which  components  precipitate  simultaneously 
by  the  formation  of  mixed  crystals,  or  by  adsorption,  or  mechanical  entrapment.  It 
generally  produces  an  intimate  mixture  of  precipitates.  The  mixture  is  calcined  at  an 
elevated  temperature  to  produce  the  desired  chemical  composition  and  crystalline  phase. 

1.3.2  Common  Characterization 

Powder  X-ray  diffraction  (XRD)  is  an  essential  technique  for  phase  identification 
and  crystal  structure  refinement.  A scanning  electron  microscope  (SEM)  in  conjunction 
with  energy  dispersive  spectroscope  (EDS)  is  used  to  examine  microstructure  and 
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chemical  features.  A transmission  electron  microscope  (TEM)  is  used  to  characterize  the 
mierostructure  at  higher  magnifieation. 

In  addition,  vibrational  and  nuclear  magnetic  resonance  (NMR)  speetroscopes  are 
very  important  methods,  especially  in  the  study  of  solid  phosphate.  These  methods 
eharacterize  the  local  environments  around  phosphates. 

1. 3.2.1  Infrared  (IR)  and  Raman  spectroscopy 

Each  atom  has  three  degrees  of  freedom  due  to  its  motions  along  the  Cartesian 
coordination  axes.  Therefore,  a molecule  with  N atoms  has  3N  degrees  of  freedom,  which 
describe  the  translation,  rotation  and  vibration  of  molecules.  Translation  refers  to  the 
motion  of  the  entire  molecule  through  the  space  and  it  has  3 degrees  of  freedom.  Rotation 
refers  to  the  rotation  of  the  entire  moleeule  and  it  has  3 degrees  of  freedom  for  nonlinear 
molecules  and  2 degrees  of  freedom  for  linear  moleeules.  Vibration  refers  to  the 
continuous  vibrations  of  all  atoms  in  the  molecules  respective  to  each  other.  It  has  3N-6 
degrees  of  freedom  for  nonlinear  molecules  and  3N-5  degrees  of  freedom  for  linear 
molecules. 

Molecule  vibrational  spectroscopy  is  relatively  complicated.  Quantum  mechanics 
only  allows  certain  vibrational  frequencies,  which  are  known  as  the  normal  molecule 
vibrational  modes.  There  are  several  types  of  motions  and  some  examples  are 

• Stretching  motion  between  two  bonded  atoms 

• Bending  motion  between  three  atoms  connected  by  two  bonds. 

Infrared  spectroscopy  characterizes  molecule  vibrations  by  measuring  the  light 

absorption,  whose  frequency  matches  the  vibrational  excitation  energy  of  a molecule 
from  low  energy  states  to  higher  energy  states.  For  a vibration  to  be  IR  active,  the 
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molecule  dipole  moment  must  change.  The  Raman  effect  is  due  to  the  interaction  of  a 
molecule  with  the  electromagnetic  field  of  the  incident  radiant.  For  a vibration  to  be 
Raman  active,  the  molecule  polarizability  must  change  [39]. 

Infrared  spectroscopy  and  Raman  spectroscopy  are  used  to  study  the  local 
symmetry  of  the  phosphate  groups.  Characteristic  group  frequencies  of  orthophosphates 
and  diphosphates  are  mentioned  in  many  references  [40-42].  Both  IR  and  Raman  spectra 
of  orthophosphates  have  no  band  at  700-900cm'\  The  IR  spectrum  of  a tetrahedral 
symmetric  orthophosphate  anion  is  characterized  by  a strong  stretching  mode  at  around 
1000  cm‘^  and  a strong  bending  mode  at  500-600  cm''.  Diphosphates  are  distinguished 
from  orthophosphates  by  the  P-O-P  bridge  vibrations.  For  diphosphates,  the  asymmetric 
P-O-P  stretching  mode  at  900-950cm  ’ is  strong  in  the  IR  spectrum  and  the  symmetric  P- 
0-P  stretching  mode  around  750cm‘'  is  strong  in  the  Raman  spectrum.  However, 
vibrations  are  no  longer  characteristic  of  the  anion  alone  if  the  M-0  (P)  bonds  in  the 
lattice  are  strong. 

Table  1-1  lists  the  infrared  and  Raman  stretching  frequency  ranges  for  phosphate 
anions  in  the  solid  state  [40].  Table  1-2  also  lists  the  infrared  and  Raman  characteristic 
group  frequencies  of  acid  phosphates  [41].  Table  1-3  lists  the  Fourier  transform  Raman 
shifts  and  assignments  of  Sr-doped  LaP04,  which  reveals  that  HPOl^  and  ^^ist  in 

Sr-doped  LaP04  under  wet  and  dry  conditions,  respectively,  in  addition  to  POl'  [42]. 

1. 3.2.2  Nuclear  Magnetic  Resonance  spectroscopy 

When  an  atom  is  placed  in  a magnetic  field,  its  electrons  will  circulate  about  the 
direction  of  the  applied  magnetic  field.  This  circulation  causes  a small  opposite  magnetic 
field  at  the  nucleus.  The  effective  magnetic  field  at  the  nucleus  is  therefore  generally  less 
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than  the  applied  field.  This  phenomenon  is  called  nuclear  magnetic  resonance  (NMR). 
NMR  spectroscopy  is  the  use  of  the  NMR  phenomenon  to  study  physical,  chemical  and 
biological  properties  of  materials  [40], 

The  electron  density  around  each  nucleus  varies  according  to  the  types  of  nuclei 
and  bonds  in  the  molecule.  The  opposing  field  and  the  effective  field  at  each  nucleus  will 
vary.  This  is  called  the  chemical  shift  phenomenon  and  related  with  local  environments. 
The  chemical  shift  of  a nucleus  is  the  difference  between  the  resonance  frequency  of  the 
nucleus  and  reference.  The  chemical  shift  is  given  the  symbol  5,  which  can  be  expressed 
as 

S = (1-10) 

where  u and  ^re  resonance  frequencies  of  the  nucleus  and  reference,  respectively. 

MAS-NMR  spectroscopy  was  applied  to  investigate  the  defect  structures  in  Sr 
doped  LaP04  [42].  It  was  carried  out  with  the  Bruker  MSL-300  spectrometer  with  an 
operating  frequency  of  300.13  MHz  for  and  121.50  MHz  for  The  spinning  speed 
was  10  kHz.  The  chemical  shifts  for  'H  and  were  given  with  respective  to 
tetramethylsilane  (TMS)  and  85%  phosphoric  acid,  respectively.  Table  1-4  shows  the 
chemical  shift  of 'H  and  ^'P  in  Sr  doped  LaP04  [42]. 

1.3.3  Electrochemical  Characterization 
1. 3.3.1  Conductivity  measurement 

When  a constant  electric  field  is  applied  to  material,  a steady-state  current  may 
result.  The  ratio  of  this  current  to  the  applied  electric  field  is  proportional  to  a material 
property,  which  is  known  as  conductivity.  Conductivity  is  defined  as  the  conductance  of 
a cube  of  material  of  unit  cross  section.  In  ceramics,  ionic  conduction  is  caused  by 
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migration  of  ions  in  an  electric  field  and  electronic  conduction  is  caused  by  the  migration 
of  mobile  electronic  carriers.  Both  ionic  and  electronic  conductivity  contribute  to  the  total 
conductivity.  The  total  conductivity  is  given  by 


where  is  the  electric  mobility,  c,  is  the  number  of  mobile  carriers  per  unit  volume,  z,- 
is  the  net  charge  of  the  moving  ion. 

Conductivity  can  be  measured  by  several  techniques  such  as  four-point  method, 

AC  impedance  measurement  and  I-V  relations.  Small  amplitude  AC  impedance 
measurement  can  be  used  to  determine  the  bulk,  grain  boundary  and  the  electrode 
characteristics.  Their  responses  can  be  described  as  a capacitor  connected  with  a resistor 
in  parallel.  Parameters  such  as  oxygen  partial  pressure,  sample  length  and  grain  size  can 
he  varied  to  identify  the  bulk,  grain  boundary  and  the  electrode  effect.  Raising  the 
temperature  will  enhance  the  contribution  of  the  process  with  the  higher  activation 
energy.  Variation  in  oxygen  partial  pressure  influences  the  electrode  parameter  and  can 
affect  the  bulk  conductivity,  depending  on  defect  equilibrium.  Increasing  the  sample 
length  increases  the  bulk  impedance.  Increasing  the  grain  size  reduces  the  impedance 
associated  with  the  total  grain  boundary. 

1.3.3.2  Transference  number  measurement 

The  transference  number  characterizes  the  contribution  of  each  mobile  carrier  to 
the  current  in  a material.  The  knowledge  of  for  all  components  i is  very  important  for 
understanding  the  electrolytic  process.  It  is  defined  as 


^ total 


(1-11) 


tt  = 


(1-12) 
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In  proton  conducting  ceramics,  total  conductivity  includes  ionic  conductivity  ( cr,„„ ) 
and  electronic  conductivity  (0-^,^^).  Ionic  conductivity  includes  contributions  from 
oxygen  vacancies  Vq  and  protons  OHq  . Electronic  conductivity  includes  contributions 

from  electrons  e'  and  electron  holes  h . Therefore,  the  contribution  of  ionic  conductivity 
and  electronic  conductivity  can  be  expressed  as 


(^,otal  =(^ion  +<^elec  =<^ion  + +0".^  + a _■  + a 


OHn  V„ 


(1-13) 


OH„ 


OHn 


' total 


(1-14) 


' total 


(1-15) 


h‘on  ^elec  ^ 


(1-16) 


where  is  the  ionic  transference  number,  is  the  electronic  transference  number. 


is  the  proton  transference  number  and  is  the  oxygen  ionic  transference  number. 

There  are  many  methods  to  determine  the  transference  number,  such  as 
measurement  of  permeation,  Hebb-Wagner  polarization  and  electromotive  force  (EMF) 
measurement.  In  EMF  measurement,  the  material  is  used  as  an  electrolyte  in  a galvanic 
cell.  If  it  is  a pure  ionic  conductor,  EMF  of  the  cell  is  given  by  the  equation  AG  = -nFE , 
where  AG  is  the  Gibbs  free  energy  of  the  cell  reaction,  F is  the  Faraday  constant  and  n 
is  the  number  of  Faraday  equivalents  which  flow  through  the  cell.  If  this  material  also  has 
electronic  conductivity,  the  experimental  value  of  the  EMF  is  smaller  than  the  theoretical 

£ 

one.  Therefore,  ionic  transference  number  is  obtained. 

^ theoretical 
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1. 3.3.3  Hydrogen  permeation 

Hydrogen  permeation  can  be  defined  as  a process  by  which  pure  hydrogen  is 
selectively  passed  through  a separation  membrane  from  a gas  mixture.  The  hydrogen 
permeation  process  is  a kinetic  process,  which  is  controlled  by  the  bulk  transport  and 
surface  exchange.  The  exchange  of  hydrogen  between  the  surface  and  the  gas  phases 
involves  adsorption,  dissociation,  charge  transfer,  surface  diffusion  of  intermediate 
species,  and,  finally,  incorporation  in  the  surface  layer.  During  bulk  transport,  the 
hydrogen  is  transported  through  the  membrane  in  the  form  of  protons,  rather  than 
molecules,  under  the  driving  force  of  a gradient  in  hydrogen  chemical  potential.  Protons 
are  charge  compensated  by  a simultaneous  flux  of  electrons. 

Generally,  the  hydrogen  permeation  flux  across  an  oxide  membrane  can  be 


calculated  using  the  Wagner  equation,  which  assumes  the  bulk  diffusion  to  be  the  rate- 
limiting  step  [35]. 


OH, 


\ RT 


'total' 


OH„ 


din  Pq  + 


RT 

2F^ 


where  L is  the  thickness  of  the  membrane,  R is  the  universal  gas  constant;  F is  the 
Faraday  constant  and  dlnP  is  the  chemical  potential  gradient  across  the  membrane. 
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. Total  conductivity  of  BCGO  and  Pd-BCGO  composite  (from  ref  [36]). 
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Figure  1-2.  Conductivity  of  various  systems  in  wet  10vol%  hydrogen,  (a)  Bulk 
conductivity  and  (b)  grain  boundary  (from  ref  [36]). 
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Figure  1-3.  Hydrogen  flux  of  BCGO  with  time  at  650  “C  ( from  ref  [36]). 
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4.  Hydrogen  flux  through  BCGlPd  with  time  at  650  C (from  ref  [36]). 


Figure  1-5.  Hydrogen  flux  through  BCG5Pd  with  time  at  650  C (from  ref  [36]). 
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Figure  1-6.  Hydrogen  flux  as  a function  of  temperature  through  different  system  (from  ref 
[36]). 
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1000/T 

Figure  1-7.  Conductivity  in  some  proton  conductors  as  a function  of  temperature  (from 
ref  [38]). 
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Table  1-1 . Infrared  and  Raman  stretching  frequency  ranges  for  phosphate  anions. 


Anion 

Infrared  (cm'*) 

Raman  (cm'*) 

POl' 

1000-1070,  strong,  asymmetric,  PO4 

930-970,  strong,  symmetric,  PO4 

PiO^- 

720-750,  weak,  symmetric,  POP 
900-950,  strong,  asymmetric,  POP 
1070-1120,  strong,  asymmetric,  PO3 

720-750,  strong,  symmetric,  POP 
1010-1020,  strong,  symmetric,  PO3 

Data  from  ref  [40]. 


Table  1-2.  Infrared  and  Raman  functional  group  frequencies. 


Functional 

Region 

Intensity 

Group 

-1 

cm 

IR 

Raman 

POl~ 

475 

medium 

1180-1000 

strong 

strong 

1000-900 

strong 

580-540 

very  week 

weak-medium 

415-380 

weak-medium 

PiO^- 

1220-1060 

strong 

medium-strong 

1060-960 

weak-medium 

strong 

980-850 

medium 

weak 

770-700 

medium-weak 

weak-medium 

600-500 

medium-strong 

weak-medium 

530-400 

weak-medium 

355-315 

weak 

HPOl~ 

2900-2750 

medium 

medium-weak 

2500-2150 

weak 

strong 

1900-1600 

weak 

1410-1200 

weak 

1220-1100 

medium-strong 

1150-1000 

strong 

strong 

1110-925 

medium-weak 

medium 

920-825 

medium-weak 

weak 

580-450 

weak 

430-390 

weak 

H^PO~, 

2700 

weak 

medium-weak 

2400-2200 

weak 

strong 

1700 

weak 

1250 

strong-medium 

1200-950 

strong 

medium-strong 

950-850 

strong 

strong 

580-540 

weak 

medium 

450-350 

medium-strong 

Data  from  ref  [41]. 
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Table  1-3.  Reported  Raman  shifts  of  Sr  doped  LaP04. 


Functional 

Group 

Raman  (cm‘) 

HPO]- 

1135  and  672 

P-O-H  in  plane  and  out-of-plane 
deformation  mode 

P2O'- 

1036,  764  and  728 

po\- 

1073, 1065,  1055  and  1025 

asymmetric  stretching  mode 

991  and  967 

symmetric  stretching  mode 

619,  589,571,559  and  537 

O-P-0  F2-bending  mode 

466,416  and  395 

O-P-0  E-bending  mode 

Data  from  re 

'[42]. 

Table  1-4.  Chemical  shift  of  NMR  and  NMR  in  Sr  doped  LaP04. 

High  resolution  H MAS-NMR 

Water  molecules 5.4  and  1.2  ppm 

Proton  in  the  crystals  12.4  ppm,  relaxation  time  is  shorter,  intensity 

increases  with  doping  concentration 

High  resolution  P MAS-NMR 

Orthophosphate  ions  POl~  -4.4  ppm 

Protons  near  phosphorus  -30.6  ppm,  relaxation  time  is  shorter,  intensity 

increases  with  doping  concentration 

Pyrophosphate  P2O7  -7.2  and -8.4  ppm 


Data  from  ref  [42]. 


CHAPTER  2 

REVIEW  OF  PROTONIC  CONDUCTORS 

2.1  Introduction 

As  pointed  out  in  the  first  chapter,  if  the  tertiary  phase  in  composite  BCGSPd  is  a 
good  proton  conductor  between  300  and  600  °C,  it  would  be  an  important  development  in 
fuel  cells  and  hydrogen  separation.  BCGSPd  was  synthesized  by  a unique  process, 
electroless  deposition.  In  this  process,  BCGO  particles  were  suspended  in  an  electroless 
solution  containing  palladium  chloride  and  ammonium  hydroxide.  Hypophosphite  was 
used  as  the  reducing  agent  to  deposit  Pd  on  the  surface  of  BCGO  particles  from  an 
electroless  bath.  The  coated  particles  were  separated  from  solution  by  filtration  and 
washed  several  times  to  remove  salts.  However,  some  impurity  was  still  absorbed  during 
the  process.  Based  on  information  about  the  synthesis  process  of  BCGSPd,  the  tertiary 
phase  is  more  likely  to  be  a phosphate.  Phosphates  with  or  without  structure  protons  have 
been  studied  recently  as  proton  conductors  at  intermediate  temperatures.  In  this  chapter, 
the  research  status  of  proton  conduction,  defect  structure,  transference  number  and 
hydrogen  permeation  of  phosphates  is  reviewed. 

2.2  Proton  Conduction  in  Phosphates 
Most  recent  research  about  proton  conduction  in  phosphates  focused  on  doped 
lanthanum  phosphate,  LaP04.  According  to  JCPDS  32-0493,  LaP04  crystallizes  in  space 
group  P2]/n  (14)  with  lattice  parameter  a=6.8366  A,  b=7.0769  A,  c=6.5095  A and  p= 
103.23°.  It  has  a low  symmetry  monoclinic  monazite  structure.  Figure  2-1  approximately 
illustrates  the  relative  position  of  elements  and  the  bond  lengths  in  the  structure  [43].  The 
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monazite  structure  is  composed  of  PO4  tetrahedra  and  La09  edge  sharing  polyhedra, 
which  is  clearly  shown  in  figure  2-2  [44],  In  the  reported  lanthanum  phosphate  studies, 

La  sites  were  usually  substituted  by  alkaline  earth  metals  or  by  a composition  of  rare 
earth  elements.  The  substitution  leads  to  the  formation  of  oxygen  vacancies  for  charge 
balance.  The  protons  are  introduced  into  the  lattice  through  the  equilibrium  between 
oxygen  vacancies  and  water  vapor.  The  conductivities  of  doped  LaP04  were  usually 
measured  in  air  or  oxygen.  Their  conductivities  were  higher  under  the  wet  atmosphere 
than  under  the  dry  atmosphere,  indicating  dominant  proton  conduction  in  the  presence  of 
water  vapor.  Doped  LnP04  (Ln=  Pr,  Nd  and  Sm)  and  YPO4  were  also  investigated  [45- 
51].  No  literature  reported  proton  conductivity  data  measured  in  a hydrogen  atmosphere. 
Compared  with  perovskite  type  oxides,  conductivities  of  doped  phosphates  were  lower 
when  temperature  was  higher  than  300  °C.  The  following  are  some  details  of  recent 
research. 

2.2.1  Proton  Conduction  in  Doped  LaP04 

Norby  studied  proton  conduction  in  Ca  and  Sr  substituted  LaP04  in  air,  wet  air  and 
air  with  D2O  [26-27].  The  magnitude  of  the  isotope  effect,  although  smaller  than  the  ratio 
of  1.37  expected  for  classical  free  proton  jump  migration,  was  high  enough  to  identify  a 
major  contribution  from  proton  conductivity.  The  proton  conductivity  went  through  a 
maximum  as  a function  of  temperature.  This  maximum  occurred  at  approximately  950 
and  850  "C,  for  Ca  and  Sr  substituted  sample  in  wet  air,  respectively.  At  800  °C,  the 
proton  conductivity  was  about  6xl0'^  and  3x10“'*  S cm"'  for  Ca  and  Sr  substituted  LaP04 
in  wet  air,  respectively.  The  activation  energy  for  proton  mobility  was  approximately  0.8 


eV. 
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Gallini  took  impedance  spectra  on  the  2.5  mol%  and  5 mol%  Sr  doped  LaP04  in 
wet  and  dry  oxygen  [45],  The  conductivity  of  5 mol%  doped  LaP04  was  around  2x10'^ 
S/cm  at  600  °C  in  the  wet  oxygen,  higher  than  that  of  5 mol%  doped  LaP04.  The  bulk 
resistance  dominated  above  400  °C  and  grain  boundary  resistance  dominated  below  it. 
Water  vapor  was  incorporated  to  hydrolyze  the  group  into  two  HPO}~  groups.  The 
entropy  and  enthalpy  changes  of  the  incorporation  were  -90  J/mol  K and  -83  KJ/mol, 
respectively. 

Amezawa  investigated  electrical  conduction  in  LaP04  doped  with  1 mol%  of 
alkaline  earth  metal  M (M=  Mg,  Ca,  Sr  and  Ba)  at  500-925  “C  [46-47].  Conductivities  of 
1 mol%  Sr,  Ca , and  Ba  doped  LaP04  in  wet  oxygen  at  500  ”C  are  5x10'^,  4x10'^  and 
2x10'^  S/cm,  respectively.  The  conductivity  of  Mg  doped  LaP04  in  wet  oxygen  was 
lower,  only  2x10'^  S/cm  at  600  °C.  Ca,  Sr,  and  Ba  doped  LaP04  conducted  protons 
dominantly  under  the  wet  atmosphere  but  Mg  doped  LaP04  did  not.  The  activation 
energy  for  proton  conduction  in  the  Ca,  Sr,  Ba  doped  LaP04  at  500-600°C  are  0.95,  0.98 
and  0.79  eV,  respectively.  Their  conductivity  differences  originated  from  the  solubility 
differences  of  the  alkaline  earth  metals  in  the  LaP04.  The  solubility  limits  of  Ca,  Sr  and 
Ba  in  LaP04  are  4.2%,  1.9%  and  0.4%,  respectively. 

Amezawa  also  investigated  the  conductivities  of  Sr  doped  LaP04  with  different 
sintering  temperatures  and  different  dopant  concentrations  (Sr  = 0-7  mol%)  [42].  Sr 
doped  LaP04  sintered  at  1350  ®C  had  a second  phase,  Sr2P207,  and  lower  conductivity 
than  that  sintered  at  1200  "C.  The  conductivity  of  (x/2)/Sr2P207-(l-x)LaP04  (x==  0-0.6) 
showed  a maximum  when  x=  0.02-0.03,  which  nearly  equaled  the  solubility  limit  of  Sr  in 
the  LaP04  [48].  Sr2P207  appeared  as  a second  phase  with  sintering  temperature  1200  "C 
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when  x>0.02.  The  proton  activation  energy  of  (x/2)/Sr2P207-(l-x)LaP04  (x=  0.02)  was 
0.86  eV  at  500-600  “C. 

Kitamura  investigated  the  electrical  conduction  in  1 mol%  Sr  doped  (Lai.xSmx)P04 
(x=^  0.3,  0.5  and  0.7)  [49].The  conductivities  of  1 mol%  Sr  doped  (Lai-xSmx)P04  (x=0.3, 
0.5,  0.7)  are  at  the  same  magnitude  -5  at  600  “C  in  wet  oxygen  atmosphere  and  did  not 
increase  much  with  the  increased  Sm  content.  The  activation  energies  are  0.975,  0.91, 
0.92,  0.94  and  0.96  eV  for  x=0,  0.3,  0.5,  0.7  and  1,  respectively. 

2.2.2  Proton  Conduction  in  Other  Phosphates 

Ca  doped  YPO4  has  a tetragonal  xenotime-type  structure.  Amezawa  investigated  its 
electrical  conductivity  under  0.4-8  kPa  of  water  partial  pressure  and  0.01-100  kPa  of 
oxygen  partial  pressure  [50].  The  conductivity  of  the  3 mol%  Ca  doped  YPO4  was  4xl0'^ 
S/cm  at  500°C  in  wet  conditions.  The  proton  activation  energy  was  estimated  as  0.84  eV 
at  500-700  ”C. 

LnP04  (Ln=La,  Pr,  Nd  and  Sm)  has  a monazite  structure.  Kitamura  measured  the 
electrical  conductivity  of  LnP04  and  1 mol%  Sr  doped  LnP04  [51].  The  conductivities  of 
1 mol%  Sr  doped  LnP04  under  wet  oxygen  atmosphere  were  10'^  ^ S/cm  at  500  °C  and 
were  about  10  times  higher  than  those  of  undoped  LnP04.  The  conductivities  of  undoped 
LnP04  were  almost  independent  of  humidity.  The  conductivities  of  1 mol%  Sr  doped 
LnP04  were  drastically  increased  under  wet  atmosphere  compared  with  those  under  dry 
atmosphere.  Conductivities  of  1 mol%  Sr  doped  LnP04  (Ln=  La,  Pr,  Nd  and  Sm)  were 
almost  independent  of  rare  earth  elements.  The  activation  energies  at  500-600  °C  were 
estimated  to  be  0.98,  0.91,  0.92  and  0.96  eV  for  Ln=La,  Pr,  Nd  and  Sm,  respectively. 


24 


Baranov  reported  conductivity  of  single  crystal  Pb^(,PO^)2  increased  20-40  times 
the  original  value  after  annealing  in  wet  nitrogen  [52],  Transmission  IR  spectra  of 
Pb2,(,PO^)2  after  annealing  in  wet  nitrogen  or  air  showed  a broad  band  at  3250  and  a 
doublet  at  2750-2800  cm'’,  which  were  characteristic  of  the  stretching  vibration  of  two 
different  types  of  hydrogen  bonding. 

Haile  reported  conductivity  of  a new  compound  /3-Cs2{HSO^)2{H  x{P,S)0^)  in 
ambient  air  and  dry  argon  from  room  temperature  to  the  thermal  decomposition 
temperature  190  °C  [53-54].  The  compound  went  through  a superprotonic  transition  at 
125  °C.  Its  conductivity  jumped  from  2x10'^  to  6x10'^  S/cm.  The  activation  energy  was 
0.70  eV  at  low  temperature  and  0.37  eV  at  high  temperature.  Its  low  temperature 
conductivity  was  about  one  magnitude  greater  than  that  of  a-Cs2{HSO^)2{H2PO^)  [55] 
and  about  two  orders  of  magnitude  greater  than  that  of  CsHSO^  [56].  She  also  studied 
Cs^{HSO^)2(H2PO^)2  [57]  and  Cs2(HSO^){H2PO^)  [58].  The  conductivity  of 
Cs2(HSO^){H2PO^)  in  the  high  temperature  phase  at  1 10  “C  was  3x10'^  S/cm  and  the 
activation  energy  for  proton  transport  was  0.37  eV. 

Otomo  carried  out  conductivity  measurements  of  polycrystalline  CS//2PO4  and 
CsH 2 PO^/ silica  composite  from  150  to  250  “C  under  dry  and  wet  argon  atmosphere  [59]. 
The  melting  point  of  CsH2PO^  was  345  "C.  A sudden  ionic  conductivity-jump  took  place 
at  around  230  "C.  D-H  isotopic  exchange  suggested  that  C5//2PO4  was  a proton 
conductor.  The  conductivity  maintained  the  constant  value  1x10'^  S/cm  during  20  h in 
wet  atmosphere  but  it  decreased  from  1 xlO"^  S/cm  with  time  in  dry  atmosphere  at  250  °C. 
In  Cs// 2 FO4 /silica  composite  systems,  the  conductivity  was  influenced  by  hydrophilicity 
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and  the  specific  surface  area  of  silica.  The  hydrophilic  silica  accelerated  the  conductivity 
at  the  temperature  approximately  below  the  phase  transition. 

Alberti  also  mentioned  other  phosphates  exhibited  proton  conducting  behavior  only 
when  hydrated.  The  proton  conductivity  of  HUOjPO^  AHjO  at  room  temperature  was 
5x10'^  S/cm  and  activation  energy  was  0.35  eV  [60-61].  The  proton  conductivity  of 
a-Zr{HPO^)2  nH^O  was  IxlO"^  S/cm  and  its  activation  energy  was  0.3  eV  at  room 
temperature  [62].  Its  conductivity  increased  to  1.5x10  '*  S/cm  with  97%  relative  humidity 
at  100  °C  [63].  The  proton  conductivity  of  y-Zr{P0^){H2P0^)-2H20  was  3x10''*  S/cm 
and  its  activation  energy  was  0.24  eV  at  room  temperature  [64].  Its  conductivity 
maintained  3x10"'*  S/cm  with  95%  relative  humidity  at  100  "C  [65].  The  proton 
conductivity  of  polybenzimidazole-xH3P04  was  4x10'^  S/cm  with  11%  relative  humidity 
at  190  °C  [66]. 

2.3  Proton  Incorporation  Mechanism  in  Phosphates 

The  proton  incorporation  mechanism  in  phosphates  is  still  not  clearly  understood. 
Norby,  Amezawa  and  Haile  proposed  the  following  defect  structures  based  on  their 
experimental  results. 

2.3.1  Norby’s  Hypothesis 

Due  to  the  complexity  of  the  proton  incorporation,  Norby  proposed  two 
mechanisms.  He  initially  assumed  that  dopant  Ca  or  Sr  dissolved  on  the  La  site  in  LaP04 
and  the  native  compensating  defect  was  an  oxygen  vacancy  [26-27]. 

The  electroneutrality  could  be  written  as 

Wo^+[Hi]  = [M\„]=  constant  (M=  Ca,  Sr).  (2-1) 

The  external  equilibrium  between  oxygen  vacancies  and  protons  can  be  written  as 
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Vo+H^O{g)  = 2H,+Ol 

(2-2) 

(2-3) 

Oxygen  vacancies  were  considered  to  be  dominant  at  high  temperature  and  the 
electroneutrality  in  (2-1)  was  simplified  as 

2[Fo  ] = [Mi J = constant  (M=  Ca,  Sr).  (2- 1 a) 

Equation  (2-3)  can  then  be  written  as 

= (2-3a) 

Proton  concentration  [ //;  ] should  be  proportional  to  based  on  equation  (2-3a). 
Experimental  results  showed  that  the  conductivities  of  Ca  or  Sr  doped  LaP04  were 
approximately  p'/,lo  dependent  at  800  “C  at  the  high  water  partial  pressure. 

Because  the  experimental  conductivity  dependence  of  water  partial  pressure  is 
lower  than  the  theoretical  dependence  deduced  from  his  first  hypothesis,  Norby  came  up 
with  an  alternative  defect  model  comprised  of  protons  and  oxygen  interstitials,  which  are 
formed  through  the  reaction: 

H^O{g)  = 2H,+o:  (2-4) 

Kh=[Hif[o:]p-„\o.  (2-5) 

The  proton  concentration  is  twice  that  of  the  oxygen  interstitial  concentration  based 
on  equation  (2-4).  Equation  (2-5)  can  then  be  written  as 

[//,:]  = (2/:,  (2-5a) 

Proton  concentration  should  be  proportional  to  P^hIo  based  on  this  model.  It 
matched  the  conductivity  behavior  better  than  the  first  model.  However,  Norby  thought 
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this  model  required  higher  defect  concentration  than  those  from  Ca  or  Sr  substitution  but 
he  did  not  explain  it  in  detail  [27], 

Actually,  both  of  the  two  models  have  limitations.  He  did  not  point  out  the  existing 
form  of  protons  in  the  phosphate  and  proposed  the  defect  structures  only  with 
conductivity  information.  Corresponding  structural  characterizations  should  be  done  to 
assure  his  defect  model.  He  considered  that  oxygen  vacancies  dominate  at  high 
temperature  in  his  first  hypothesis.  However,  it  is  not  absolutely  true  in  the  presence  of 
water  vapor.  The  proton  conduction  sometimes  still  dominates.  There  is  also  a possibility 
of  p-type  electronic  conduction  dominance.  His  second  hypothesis  matched  his 
experimental  results  better  but  it  did  not  indicate  the  doping  effect  on  the  possible  proton 
concentration  in  the  lattice.  If  we  consider  that  the  second  defect  structure  is  right,  then 
the  electroneutrality  could  be 

p + [H,]  = 2[0:]  + [M'^^]  + n.  (2-6) 

However,  he  considered  that  the  proton  concentration  was  twice  that  of  the  oxygen 
interstitial  concentration  based  on  equation  (2-4).  Generally,  the  oxygen  interstitials  are 
not  easily  incorporated  into  the  lattice  due  to  their  large  size.  Therefore,  the  amount  of 
protons  produced  would  not  be  as  much  as  the  model  suggested. 

2.3.2  Amezawa’s  Hypothesis 

Later,  Amezawa  came  up  with  a better  hypothesis  than  Norby.  He  sintered  a pellet 
of  Sr  doped  LaP04  (Sr=  0-7  mol%),  ground  it  into  powder  and  heat  treated  it  at  900  “C 
for  10  h in  an  unspecified  atmosphere  under  two  different  conditions.  One  condition  is  1 
kPa  of  water  vapor  and  the  other  is  vacuum  [42].  He  used  FT-Raman  and  MAS-NMR 
spectroscopy  to  examine  local  structures  in  Sr  doped  LaP04.  Those  results  revealed  that 
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hydrogen  phosphate  ions  HPOl~  and  pyrophosphate  ions  P20^~  existed  in  Sr  doped 
LaP04  under  these  two  conditions,  respectively,  in  addition  to  orthophosphate  ions  PO\~ . 
No  pyrophosphate  ion  PiO^"  was  detected  in  the  water  vapor  treated  sample  by  Raman  or 
NMR  spectroscopy.  The  chemical  shift  assigned  to  hydrogen  phosphate  ions  HPOl'  also 
appeared  in  the  5 mol%  and  7 mol%  Sr  doped  LaP04  after  heat  treatment  under  vacuum. 
However,  their  intensities  decreased  compared  with  sample  after  heat  treatment  under 
water  vapor.  He  explained  the  results  as  decreasing  concentration  of  protons  during  heat 
treatment  under  vacuum,  but  did  not  mention  the  origin  of  proton  in  the  sample. 

He  proposed  a defect  structure  model  for  Sr  doped  LaP04  based  on  the 
characterization  results.  According  to  the  model.  La  substitution  with  Sr  leads  to 
formation  of  Pj02~  as  intrinsic  positive  defects.  Protons  are  introduced  into  the 
phosphate  forming  hydrogen  phosphate  groups  through  the  equilibrium  between  P^O*" 
and  water  vapor.  Figure  2-3  is  drawn  to  approximately  illustrate  the  chemical  structure  of 
phosphate,  pyrophosphate  and  hypophosphate.  The  proton  transfer  between  phosphates  is 
also  shown.  Mobile  protons  can  transfer  to  the  nearest  neighboring  tetrahedron  or  to  the 
next  oxygen  in  the  same  tetrahedron.  The  reorientation  of  the  PO4  tetrahedron  facilitates 
the  proton  transfer. 

The  specific  relationship  between  defect  species  from  Amezawa  ‘s  paper  are 
reproduced  in  equations  (2-7)  to  (2-11)  below; 


(2-7) 


-<.PiOi)2Po,  +-H20{g)-^{HPO,)po, 


(2-8) 


K,  = [{HPO,  )2Po,  Ph^o 


(2-9) 
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\iPiOi)2Po,  +\o^[g)^{PO,Ypo,  +h 


K^-p[{P20■,)■,poJ""Po^' 


(2-10) 

(2-11) 


where 


(PO^Ypq^  : phosphate  on  phosphate  site 

{PiOjijpo^ : oxygen  vacancy  in  double  phosphate  site 

{HPOYpo^  ■ interstitial  proton  in  phosphate. 

An  issue  with  these  formulas  is  that  half  oxygen  atoms  have  no  physical  meaning. 
Therefore,  the  formulas  were  multiplied  by  a factor  of  2 and  rewritten  as  follows: 


Sp2 PjO-j  —>  + (/j O 7 ) 2PO4 


(2-7a) 


{P2O-!  )2PO,  ^ 2^(S)  pQ^ 


(2-8a) 


(^2^7)2^04  2^2!^)  2(PC>4 );.(3^  + 2h 


The  electroneutrality  conduction  can  be  written  as 


[{HPOYpo,  ] + 2[(P207 )2po  ] + P = [Sru,]  + n. 


(2- 10a) 


If  P2O2  and  HPOl  are  the  dominant  positive  defects  and  electronic  defects  are 
only  present  in  minor  concentration,  the  electroneutrality  conduction  can  be 
approximated  by 


[{HPO 4 ) pQ^  ] + 2[{P lOy) 2Po^  ] - ] constant. 


(2-12) 


He  mentioned  two  cases.  In  the  first  case,  P2O2  is  the  dominant  defect  and  charge 


neutrality  can  be  written  as 


2[(/2<9?  )2/>04  ] - ] • 


(2-13) 


Based  on  formulas  (2-9)  and  (2-13),  the  proton  concentration  can  be  described  as 
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[{HPO,),o,  ] = K,  {{Sr^yiy^  . (2-14) 

Based  on  formulas  (2-11)  and  (2-13),  the  electron  hole  concentration  can  be 
described  as 

P = Ky[Sr[„]l2f^p^^^ . (2-15) 

In  this  case,  proton  concentration  is  supposed  to  be  proportional  to  and  p-type 
electronic  concentration  is  supposed  to  follow  pq^  . 

In  the  second  case,  hydrogen  phosphate  groups  are  the  dominant  defects  and  charge 
neutrality  can  be  written  as 

{{HPO,)po,]  = [Sr^^].  (2-16) 

Based  on  formulas  (2-9)  and  (2-16),  the  oxygen  vacancy  concentration  can  be 
described  as 

mOn)2Po,  ] = ]'  Ph\o  ■ (2- 1 7) 

Based  on  formulas  (2-11)  and  (2-17),  the  electron  hole  concentration  can  be 
described  as 

p = (2-18) 

In  this  case,  p-type  electronic  concentration  is  supposed  to  follow  pq^  . Proton 

concentration  is  supposed  to  be  independent  of  water  pressure  and  oxygen  partial 
pressure. 

However,  for  Sr  doped  LaP04,  oxygen  partial  pressure  and  water  partial  pressure 
dependence  of  conductivity  is  much  smaller  than  expected.  He  considered  it  was  caused 
by  the  change  from  hydrogen  phosphate  dominance  to  oxygen  vacancy  (pyrophosphate) 


dominance. 
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Amezawa’s  hypothesis  not  only  pointed  out  the  existing  form  of  hydrogen  in 
phosphates  but  also  pointed  out  the  relationship  between  defect  species  based  on 
characterization  results.  He  predicted  the  oxygen  partial  pressure  and  water  partial 
pressure  dependence  of  conductivity  using  his  proposed  defect  structure.  Amezawa’s 
hypothesis,  however,  addressed  only  the  intermediate  and  high  oxygen  partial  pressure 
region,  ignoring  the  low  oxygen  partial  pressure  region.  He  considered  oxygen  vacancies 
or  protons  instead  of  electron  holes  dominating  in  the  high  oxygen  partial  pressure 
region.  His  NMR  and  Raman  characterization  results  did  not  indicate  the  existence  of 
pyrophosphate  the  wet  atmosphere  treated  sample.  Therefore,  the  first  case  did 

not  exist.  If  we  consider  the  undetected  pyrophosphates  were  due  to  the  limitation  of  the 
characterization  technology,  then  the  results  also  could  not  exclude  the  possibility  of 
electron  holes  as  the  dominant  species.  The  contribution  from  electron  holes  at  high 
temperature  should  not  be  omitted.  In  this  case,  we  can  consider  the  electroneutrality 
condition  as 

P=[Sr’La]-  (2-19) 

The  proton  and  oxygen  vacancy  concentration  could  be  written  as 

{{HPO,)po,  ] = [Sr'u.  (2-20) 

[iPlOl)2PO,  ] = ^2'  [Sr La  ? Poi^  ■ (2-21 ) 

This  model  also  can  explain  the  conductivity  behavior  of  Sr  doped  LaP04.  The 
electron  hole  concentration  will  not  increase  much  with  oxygen  partial  pressure  based  on 
Eqs.  (2-19),  which  might  explain  the  conductivity  dependence  better. 

He  predicted  the  conductivity  dependence  of  oxygen  partial  pressure  or  water 
partial  pressure  based  on  his  proposed  model.  Actually,  the  conductivity  is  not  only 


32 


related  with  charge  carrier’s  concentration  but  also  related  with  carrier’s  mobility.  The 
proposed  model  only  predicted  the  concentration  dependence  of  oxygen  partial  pressure 
or  water  partial  pressure. 

2.3.2  Haile’s  and  Mizuno’s  Hypothesis 

Haile  [58]  revealed  the  high  temperature  phase  of  Cs2{HSO^){H2PO^)  is  cubic  with 
a CsCl  structure  by  using  high  temperature  X-ray  powder  diffraction  and  infrared 
spectroscopy.  Cs  atoms  are  at  the  comers  of  a simple  cubic  unit  cell  and  the  XO4  group 
(X=P  or  S)  is  at  the  center.  As  a result,  oxygen  atoms  must  be  distributed  over  many  sites 
and  oxygen  site  occupancy  is  less  than  one.  Accordingly,  the  PO4  group  must  exhibit 
multiple  orientations.  The  possible  rapid  PO4  group  reorientation  can  facilitate  proton 
transfer.  Mizuno  studied  proton  dynamics  in  phase  II  of  CsHSO^  by  NMR  between 
210  and  470  K [67]  and  proton  diffusion  in  the  superprotonic  phase  of  CsHSO^  (phase  I) 
by  ’H  NMR  relaxation  [68].  The  NMR  results  demonstrate  that  proton  transfer  between 
two  neighbouring  SO4  tetrahedra  takes  place  much  easier  than  the  reorientation  of  the 
SO4  tetrahedra.  Thus,  the  reorientation  limits  the  rate  of  the  proton  transport  in  phase  II 
and  I.  The  electronic  structure  and  proton  transfer  mechanism  of  CsHSO^  phase  I were 
studied  by  the  first-principles  methods  [69].  The  calculated  result  showed  that  the  band 
gap  is  6.03  eV,  indicating  that  it  is  electronically  insulating.  The  reorientation  of  the 
sulfate  tetrahedron  can  take  place  very  frequently.  Sometimes  the  proton  transfers  to  the 
nearest  neighboring  tetrahedron  and  from  there  it  continues  to  transfer  either  to  the  next 
oxygen  ion  or  to  the  next  hydrogen  bond.  Sometimes  the  proton  rotates  to  a local 
minimum  and  from  there  it  continues  to  transfer  to  the  next  nearest  neighboring 


tetrahedron. 
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2.4  Transference  Number  in  Phosphates 

The  transference  number  characterizes  the  contribution  of  each  charge  carrier  to 
the  current  in  a material.  There  are  few  instances  in  the  literatures  where  the  proton 
transference  number  in  phosphates  is  reported. 

Gallini  showed  for  the  first  time  with  the  EMF  measurements  of  water  vapor 
concentration  cells  that  2.5  mol%  and  5 mol%  Sr  doped  LaP04  were  predominantly 
proton  conductors  below  700  “C  in  wet  oxygen  [45].  The  cell  was  exposed  to  wet  oxygen 
while  the  water  vapor  partial  pressure  was  different  on  both  sides.  The  obtained  proton 
transport  number  was  approximately  0.62  at  600  "C.  The  cell  was  then  exposed  to  an 
oxygen  partial  pressure  gradient  with  dry  gas  on  one  side  and  wet  gas  on  another  side. 
The  water  partial  pressures  on  both  sides  were  monitored.  The  obtained  total  ionic 
transport  number  is  0.59  at  600  “C.  It  was  not  significantly  different  from  the  proton 
transference  number  of  0.62.  Based  on  this,  Sr  doped  LaP04  was  considered  by  Gallini  to 
be  a proton  conductor. 

Amezawa  studied  proton  and  native  conduction  in  5 mol%  Sr  doped  LaP04  by 
thermoelectric  power  measurements  [46].  The  proton  transport  number  of  5 mol%  Sr 
doped  LaP04  was  calculated  to  be  unity  in  wet  air  at  600  °C  when  water  pressure  was 
larger  than  1 kPa  and  at  700  "C  when  water  pressure  was  larger  than  3 kPa.  He  used  the 
same  electroneutrality  conditions  as  in  the  defect  structure  studies  [42],  therefore,  the 
obtained  proton  transference  number  may  be  slightly  high. 

2.5  Hydrogen  Permeation  in  Phosphates 

When  a material  is  used  as  hydrogen  separation  membrane,  it  will  be  exposed  to 
atmospheres  in  the  low  oxygen  partial  pressure  region.  Its  hydrogen  permeability  is 
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closely  related  with  its  total  conductivity  and  transference  number  in  that  oxygen  partial 
pressure  region.  The  material  will  have  a high  hydrogen  permeation  flux  if  its  proton 
conductivity  and  electronic  conductivities  are  sufficiently  high  and  nearly  equal. 
Unfortunately,  no  one  mentioned  the  conductivity  behavior  and  transference  number 
results  of  phosphates  in  the  low  oxygen  partial  pressure  region,  especially  in  hydrogen 
atmosphere.  Thus  far,  no  literature  has  reported  hydrogen  permeability  of  phosphates. 

The  reported  conductivity  data  and  transference  number  results  of  phosphates  all 
focused  on  the  wet  air  or  wet  oxygen  atmosphere.  However,  the  conductivity  and 
transference  number  changed  with  oxygen  partial  pressures.  Therefore,  the  hydrogen 
permeation  in  phosphates  could  not  be  predicted  without  information  about  their  total 
conductivity  and  proton  transference  number  in  hydrogen  atmosphere. 

2.6  Summary 

Conductivity  and  transference  number  measurement  of  doped  phosphates  indicated 
that  they  are  proton  conductors  in  wet  oxygen  or  air  atmosphere.  The  conductivities  of 
reported  phosphates  were  still  lower  than  that  of  perovskite  oxide  system  above  300  °C. 
No  proton  conductivity  data  of  phosphates  in  hydrogen  atmosphere  was  reported.  There 
was  also  no  reported  hydrogen  permeability  of  phosphate  so  far.  The  best  proposed  defect 
structure  of  phosphate  still  has  limitations  and  only  covered  the  intermediate  and  high 
oxygen  partial  pressure  region.  As  proton  conductors,  conductivity  behaviors  of 
phosphates  and  corresponding  defect  structures  should  be  studied  in  hydrogen 
atmosphere. 

Based  on  reviewed  literatures,  we  expected  the  desired  new  phosphate  was  a proton 
conductor.  Its  conductivity  as  functions  of  temperature,  oxygen  partial  pressure  or  water 
partial  pressure  will  be  studied  and  compared  with  literature.  Its  corresponding  defect 
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structure  will  be  proposed  based  on  characterization  results.  More  attention  will  be  paid 
to  conductivity  behavior  and  defect  structure  in  the  low  oxygen  partial  pressure  region, 
especially  hydrogen  atmosphere.  Proton  transference  number  and  hydrogen  permeation 
of  the  desired  phosphate  will  be  measured  after  finding  an  appropriate  sealant  at 
intermediate  temperature  range. 
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Figure  2-1.  The  crystal  structure  of  LaP04  with  monazite  structure  (from  ref  [43]). 


Figure  2-2.  Polyhedral  diagrams  of  the  monazite  structure  (from  ref  [44]) 
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HPO4 


Figure  2-3.  Phosphate,  diphosphate  and  hydrogen  phosphate  structure  and  proton  transfer 
illustration  using  “MS  modeling”  software  and  photoshop. 


CHAPTER  3 

PHASE  IDENTIFICATION 

3.1  Introduction 

As  stated  in  the  first  chapter,  the  hydrogen  flux  of  composite  BCG5Pd  was  higher 
than  that  of  BCG35Pd  when  the  temperature  was  lower  than  550°C.  Additionally,  this 
hydrogen  flux  did  not  achieve  an  equilibrium  value  even  after  20  hours  at  650  “C.  The 
second  phase  Pd  did  not  exhibit  such  a trend  under  isothermal  conditions.  Krishnamurthy 
in  our  lab  [36]  synthesized  the  composite  BCGSPd  powder  by  electroless  deposition, 
which  is  a unique  process  and  may  introduce  some  impurity.  The  impurity  may  react  with 
composite  BCGSPd  to  form  a tertiary  phase  at  high  sintering  temperature.  The  produced 
tertiary  phase  may  contribute  to  such  an  unusual  trend.  If  we  could  identify  this  tertiary 
phase,  synthesize  it  and  prove  it  as  a good  proton  conductor,  it  will  be  an  important 
development  in  fuel  cell  and  hydrogen  research  field.  Therefore,  composite  BCGSPd  was 
characterized  by  X-ray  diffraction  and  infrared  spectroscopy.  The  electroless  deposition 
process  of  composite  BCGSPd  was  also  repeated  and  absorbed  impurity  was  measured  by 
an  inductively  coupled  plasma  atomic  emission  spectrometer  (ICP-AES).  The  pH  effects 
on  the  stability  of  BCGO  in  the  electroless  bath  as  well  as  the  impurities  absorbed  were 
studied. 

3.2.  Experiment 

3.2.1  Electroless  Deposition 

BCGO  powder  was  synthesized  by  solid-state  reaction  and  sieved  to  less  than  45 
pm.  For  the  electroless  deposition  at  pH=  3,  the  procedure  was  performed  the  same  as 
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Krishnamurthy  [36],  For  the  electroless  deposition  at  pH-  10  [70],  12g  BCGO  powder 
was  suspended  in  the  electroless  solution  bath  consisting  of  palladium  chloride  (1.8  g), 
ammonium  hydroxide  (24  mL),  ammonium  chloride  (3.6  g),  formaldehyde  (37  wt%,  16 
mL)  and  buffer  (pH=  10,  80  mL).  Sodium  hypophosphite  (8  g)  was  added  to  the  bath 
under  rapid  stirring.  In  both  pH  conditions,  the  color  of  the  electroless  baths  changed  to 
black  immediately  when  sodium  hypophosphate  was  added.  Color  of  the  bath  became 
shallow  green  and  black  particles  precipitated  to  the  bottom,  which  suggested  that  the 
bath  was  exhausted.  Powders  were  then  filtered  through  3.0  um  filter  paper  and  washed 
with  deionized  water  several  times.  Then  they  were  dried  and  calcined  at  900,  1100  and 
1475  °C  in  air  for  10  h and  finally  treated  at  900  °C  in  argon  for  5 h. 

3.2.2  Solid  State  Reaction 

A mixture  of  BCGO  and  NH4H2PO4  was  milled,  dried  and  calcined  at  1475  °C  for 
10  h and  argon  treated  at  900  °C  for  5 h.  For  conductivity  measurement,  sintered  powders 
were  uniaxially  pressed  at  35  MPa,  cold  isostatic  pressed  at  200  MPa,  sintered  at  1450  "C 
in  air  for  1 0 h and  treated  at  900  °C  in  argon  for  1 0 h. 

3.2.3  Characterization 

The  actual  concentration  of  coated  Pd  and  P impurity  absorbed  under  different  pH 
conditions  were  measured  by  ICP-AES  (Perkin-Elmer  Plasma  3200)  with  a detection 
limit  range  of  less  than  1 part  per  million.  1 g of  electroless  deposited  powder  was  dried, 
accurately  weighted,  and  then  transferred  into  a 500  ml  volumetric  flask.  10  ml  37%  HCl 
was  added  to  totally  dissolve  the  mixture  and  then  diluted  with  distilled  deionized  water. 

XRD  (PHILIPS-APD-3720)  was  used  to  analyze  the  phase  distribution  in  the 
powders.  IR  (INTCOLET  Magna  760  Bench  spectrometer)  was  used  to  investigate  the 
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existence  form  of  the  tertiary  phase.  The  measured  wave  number  range  was  400-4000  cm' 
' . The  reference  was  KBr  powder. 

The  conductivity  measurement  followed  Krishnamurthy’s  procedure  [36].  Gold 
paste  was  brushed  on  top  of  the  sintered  pellet  of  the  mixture  of  BCGO  and  NH4H2PO4, 
dried  at  125  "C  for  30  min.  The  sample  was  then  fired  at  850  “C  for  1 h.  The  gold  pastes 
acted  as  electrodes  for  the  sample.  Thin  gold  wires  were  attached  to  the  gold  electrodes  to 
make  contact  with  leads  of  the  analyzer.  The  conductivity  was  measured  by  a Solatron 
1260.  The  atmosphere  was  wet  10vol%  hydrogen,  using  temperature  range  of  400-800 
°C. 

3.3  Results  and  Discussion 

3.3.1  XRD  Results 

Figure  3-1  shows  XRD  result  of  composite  BCG5Pd  synthesized  by 
Krishnamurthy.  It  is  compared  with  that  of  composite  BCG35Pd  synthesized  by  solid- 
state  reaction.  Peak  position  of  Pd  and  BCGO  is  marked.  For  composite  BCG5Pd  before 
sintering,  in  addition  to  Pd  and  BCGO,  other  peaks  also  appeared.  However,  their 
intensities  were  so  low  that  they  can  be  considered  as  machine  noise.  The  intensity  of 
impurity  peaks  increased  after  sintering.  Impurity  peaks  at  20=33.2°  and  56.3°  belonged 
to  CeOj  (200)  and  (311).  The  first  peak’s  intensity  was  30%  relative  to  the  second  peak. 
This  could  be  used  as  a sure  indication  that  CeO^  was  present  in  the  material  under 
examination  [71].  However,  determining  the  presence  of  other  possible  products  was 
complicated  by  the  fact  that  many  peaks  of  possible  products  overlap  each  other  and  their 


intensities  were  all  low. 
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In  Krishnamurthy’s  process  [36],  BCGO  powders  were  suspended  in  an  electroless 
solution  bath  (pH=3)  consisting  of  PdCl2  and  ammonium  hydroxide.  Sodium 
hypophosphite  was  added  to  reduce  Pd  from  the  solution  under  rapid  stirring.  Powders 
were  then  filtered  and  washed.  Some  ions  were  not  easily  washed  off  and  absorbed  as 
impurities  on  the  surface  of  the  composite.  However,  XRD  is  not  very  sensitive  to  the 
surface  composition  change  compared  with  the  grain  interior  even  if  the  surface  is  well 
crystallized  [72].  Absorbed  impurities  began  to  react  with  the  interior  of  BCGO  during 
sintering  and  their  XRD  peak  intensity  increased.  CeOj  was  easily  determined  in  the 
sintered  composite  BCGSPd. 

As  we  know,  the  basicity  of  BCGO  is  advantageous  for  the  formation  of  proton 
charge  carriers,  but  makes  it  unstable  in  acidic  solution,  which  is  used  in 
Krishnamurthy’s  process  [36].  The  pH  condition  also  affects  the  amount  of  impurities 
absorbed.  Therefore,  electroless  deposition  process  of  BCGSPd  was  repeated  at  different 
pH  conditions. 

ICP-AES  analysis  results  showed  that  the  composite  was  composed  of  6.87  wt%  Pd 
and  2.89  wt%  P with  deposition  at  pH=  10  and  it  was  composed  of  7.38  wt%  Pd 
(4.31vol%)  and  2.32  wt%  P with  deposition  at  pH=  3.  Pd  concentration  was  less  than  the 
theoretical  value  9.07  wt%.  It  indicated  that  Pd  did  not  deposit  completely  onto  the 
BCGO  particle  surface  while  reduced  from  the  PdCb  solution.  Electroless  bath  with  pH= 
3 had  a better  deposition  effect  of  Pd  on  BCGO  and  less  phosphorus  adsorption.  Both  of 
these  two  electroless  bathes  decomposed  part  of  BCGO  to  form  Ce02.  This 
decomposition  is  more  severe  in  the  electroless  bath  with  pH=  3,  which  is  shown  in 
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figure  3-2.  From  the  stability  issue  of  BCGO,  it  is  better  to  use  a basic  electroless  bath  for 
composite  synthesis. 

Phosphorus  impurity  was  both  absorbed  onto  BCGO  under  two  pH  conditions.  It 
reacted  with  BCGO  to  form  some  tertiary  phase  at  higher  sintering  temperature.  Two 
samples  showed  little  difference  in  the  absorbed  phosphate  concentration  and  the 
intensity  of  Ce02  was  a little  higher  in  the  composite  with  4.31vol%Pd.  Ce02  was  an 
additional  unwanted  impurity  in  the  composite  and  made  the  tertiary  phase 
characterization  more  difficult.  However,  the  pH  condition  of  the  electroless  bath  is 
similar  between  BCG4.31Pd  and  BCGSPd.  We  chose  the  composite  with  4.31vol%  Pd 
(BCG4.31Pd)  to  do  further  tertiary  phase  characterizations  because  of  it’s  same  process 
with  Krishnamurthy’s. 

Figure  3-3  shows  the  XRD  results  of  BCG4.31Pd  with  different  sintering 
temperatures.  The  main  peak  of  BCGO  decreased  and  some  peaks  even  disappeared  after 
sintering  at  900  °C.  The  intensity  of  Ce02  peaks  increased  with  sintering  temperature. 
This  suggests  that  decomposition  of  BCGO  became  worse  due  to  the  reaction  with 
phosphorus  impurity  at  high  temperature.  In  addition  to  Ce02,  there  are  other  impurity 
peaks  and  their  intensities  increase  with  increasing  sintering  temperature.  Ba3Ce(P04)3, 
Ba3(P04)2  and  CeP04  are  all  possible  impurities.  As  pointed  out  in  chapter  1,  the  desired 
tertiary  phase  might  be  a good  proton  conductor.  The  reported  conductivities  of 
Ba3(P04)2  and  CeP04  were  low  [25]  and  therefore  may  be  just  unwanted  impurities. 
Ba3Ce(P04)3  might  be  the  desired  tertiary  phase.  The  existence  of  Ba3Ce(P04)3  could  be 
produced  by  the  solid  state  reaction  of  Ba3(P04)2  and  CeP04  at  high  sintering 
temperature,  however,  the  phase  diagram  of  system  CeP04-Ba3(P04)2  was  not  found. 
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Figure  3-4  shows  the  phase  diagram  of  system  YP04-Ca3(P04)2  [73].  Ca3YP04  is 
produced  when  the  temperature  is  higher  than  1215  °C.  Ba3Ce(P04)3,  Ba3(P04)2  and 
CeP04  may  co-exist  in  a similar  phase  diagram. 

Caballero  [72]  reported  that  residual  phosphorous  located  at  the  surface  of  the 
particles  reacted  with  BaTi03  during  sintering  to  form  BaO-rich  compounds  and  leave  a 
certain  Ti02  excess.  CeOj  excess  and  BaO-rich  compound  Ba^Ce{PO^)^  may  also  be 
produced  by  the  reaction  of  BCGO  with  other  phosphorous  containing  compounds. 
Therefore,  phosphorous  was  introduced  by  solid-state  reaction  of  BCGO  with  NH4H2PO4 
as  the  starting  material,  which  avoided  the  sodium  impurity  introduced  by  hypophosphite. 
Different  ratios  of  BCGO  to  NH4H2PO4  were  used  to  investigate  the  phosphorous 
concentration  effect  on  the  composition  of  the  product. 

In  Krishnamurthy’s  experiment,  BCG5Pd  was  sintered  at  1475  "C  in  air  and  900”  C 
in  argon  after  electroless  deposition.  Argon  treatment  was  used  to  prevent  the  formation 
of  PdO  [36].  It  is  unknown  whether  or  not  argon  atmosphere  has  an  effect  on  the 
composition  of  the  tertiary  phase.  Therefore,  the  exact  sintering  atmosphere  was  followed 
in  solid-state  reaction  of  BCGO  with  NH4H2PO4. 

Figure  3-5  shows  XRD  results  of  sintered  mixture  of  BCGO  and  NH4H2PO4  with  a 
sintering  temperature  of  1475  °C  and  dwell  time  10  h.  The  molar  ratios  of  BCGO  and 
NH4H2PO4  are  1:0.1,  1:0.3  and  1:1,  respectively.  Figure  3-6  shows  XRD  results  of 
sintered  mixture  of  BCGO  and  NH4H2PO4  with  molar  ratio  of  1 :0.1  and  1:1  with  or 
without  argon  treatment.  Compared  with  XRD  result  of  BCGO,  they  all  contained 
impurity  peaks.  There  are  few  differences  between  XRD  results  of  sintered  mixture  with 
or  without  argon  treatment,  suggesting  the  little  effect  of  atmosphere  on  the  product 
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composition.  Products  were  identified  as  Ba3(P04)2,  CeP04,  CeO^ , Ba^Ce{PO^  )j . The 
relative  intensity  of  CeOj  and  Ba^Ce(PO^)^  increased  with  initial  concentration  of 
NH4H2PO4. 

All  the  produced  mixtures  contained  one  common  composition  Ba-^Ce{PO^)^ , 

which  is  the  most  likely  the  desired  tertiary  phase.  Its  crystal  structure  is  cubic  with  a 
lattice  constant  of  a=  10.51 12  A according  to  JCPDS  43-0643. 

3.3.2  IR  Results 

Figure  3-7  (a)  shows  the  infrared  spectroscopy  of  sintered  composite  BCG5Pd. 
Figure  3-7  (b)  shows  the  infrared  spectroscopy  of  SrCe^_^Yfi^_g  and  SrZr^_^Yfi^g  [19]. 
Table  3-1  lists  their  band  assignment.  Perovskites  SrCe^_^Yfi^_g  and  SrZr^_^Yfi^__g  have 

broad  OH  stretching  bands  at  3200  cm'*.  OH  stretching  band  of  composite  BCGSPd  is 
obscured  by  a strong  and  broad  band  at  3430  cm'',  which  is  attributed  to  the  stretching 
vibrations  of  water.  The  band  at  1631  cm  ' of  composite  BCG5Pd  is  due  to  the  bending 
vibrations  of  water  molecules.  The  band  at  1035  cm''  of  composite  BCG5Pd  is  attributed 
to  the  tertiary  phase. 

Ba^Ce{PO^)^ , the  possible  tertiary  phase,  is  an  orthophosphate  with  tetrahedral 
symmetry  (Td).  Tetrahedral  symmetric  orthophosphate  anion  has  four  types  of  vibrations: 
symmetric  stretching  ( u, ),  symmetric  bending  ( ),  antisymmetric  stretching  ( ) and 
antisymmetric  bending  ( ).  They  are  illustrated  in  figure  3-8.  Symmetric  stretching  and 
symmetric  bending  are  IR  inactive.  IR  spectrum  of  a tetrahedral  symmetric 
orthophosphate  has  a strong  stretching  mode  at  around  1 000  and  a strong  bending  mode 
at  500-600  cm''.  However,  the  vibrations  of  molecules  fixed  in  the  lattice  are  not  only 
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affected  by  their  symmetry  but  also  by  the  lattice  symmetry,  which  may  shift  the  band 
position.  The  band  at  1035  cm  * of  composite  BCG5Pd  is  attributed  to  asymmetric 
stretching  of  PO2  group.  Asymmetric  bending  of  PO2  group  is  observed  at  456  cm  *. 

Figure  3-8  (a)  and  (b)  compares  IR  of  composite  BCG5Pd  with  the  high 
temperature  IR  of  Cs^i^HSO^^HjPO^) , which  was  reported  to  have  a cubic  structure  at 
high  temperature  [58].  Figure  3-8  (c)  illustrates  the  difference  between  the 
orthophosphates,  diphosphates  and  hydrogen  phosphates.  Diphosphates  are  characterized 
by  the  POP  bridge  vibrations  with  a strong  asymmetric  stretching  mode  around  900-950 
cm'*  in  ER  spectrum.  There  is  no  POP  vibration  observed  in  the  composite  BCG5Pd, 
suggesting  no  diphosphate  in  the  tertiary  phase.  Hydrogen  phosphates  are  characterized 
by  the  POH  vibrations.  Table  3-2  lists  the  wave  number  of  OH  vibrations  of 
Cs2{HSO^){H2PO^)  at  199°C,  which  are  at  1715,  2350  and  2945  cm'*.  These  OH 
vibrations  are  not  observed  in  the  composite  BCG5Pd,  suggesting  no  H2PO4  group  in  the 
tertiary  phase. 

Figure  3-9  shows  the  IR  spectra  of  the  sintered  mixture  of  BCGO  and  NH4H2PO4 
with  ratio  1:1  compared  with  those  of  composite  BCG4.31Pd  and  BCG5Pd.  Table  3-3 
lists  their  band  assignments.  The  stretching  vibration  mode  is  degenerated  in  the 
composite  BCG5Pd  but  split  in  the  composite  BCG4.31Pd.  The  stretching  vibration 
modes  of  BCG4.31Pd  are  observed  at  1060  and  1000  cm'*  while  that  of  the  sintered 
mixture  is  observed  at  1032  cm'*.  The  bending  vibration  mode  is  degenerated  in  the 
composite  BCG5Pd  but  split  in  tbe  sintered  mixture.  The  bending  vibration  modes  of  the 
sintered  mixture  are  observed  at  546  and  407  cm'*  while  that  of  composite  BCG4.31Pd  is 
observed  at  559  cm'*.  However,  POP  bridge  vibrations  are  observed  at  935  cm'*  in 
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composite  BCG4.31Pd  and  at  897  cm’’  in  the  sintered  mixture.  This  shows  that  small 
amount  of  diphosphates  were  produced  in  the  composite  BCG4.31Pd  and  sintered 
mixture. 

3.3.3  Conductivity  Measurement 

Conductivity  of  the  sintered  mixture  of  BCGO  and  NH4H2PO4  was  measured  to 
investigate  the  effect  of  the  desired  phase  Ba^Ce(PO^)^  on  the  total  conductivity.  Figure 

3-10  shows  the  impedance  plot  of  the  sintered  mixture  of  BCGO  and  NH4H2PO4  at  800 
”C  in  wet  10vol%  hydrogen  atmosphere. 

In  a typical  impedance  plot  of  a ceramic  system,  the  semicircles  belong  to  the 
electrode,  grain  boundary,  and  bulk,  from  low  frequency  to  high  frequency.  The 
corresponding  equivalent  circuit  is  composed  of  resistors  and  capacitors  in  an  ideal 
system.  In  the  system  of  the  sintered  mixture  of  BCGO  and  NH4H2PO4,  there  are  many 
produced  impurities.  It  is  difficult  to  use  a simple  electrical  circuit  to  simulate  the 
involved  physical  transport.  However,  information  about  the  bulk  conductivity  is  clearly 
demonstrated  by  semicircle  I.  The  bulk  electrical  resistance  values  were  obtained  from 
the  low  frequency  interception  of  the  semicircle  I on  the  real  axis  of  the  impedance  plot 
in  figure  3-10  (a).  Semicircle  II  is  more  complex,  including  electrode  effect  and  other 
contribution.  The  information  about  the  tertiary  phase  resistance  was  found  to  be 
impossible  to  obtain  from  this  impedance  plot. 

Figure  3-1 1 is  the  Arrhenius  plot  of  the  sintered  mixture  of  BCGO  and  NH4H2PO4. 
The  diameter  of  the  sample  was  25.92  mm  and  thickness  was  3.344  mm.  The  relative 
density  was  87%.  The  temperature  range  was  400-800  "C.  The  bulk  conductivity  was  on 
the  order  of  10'^  S/cm  at  500-800  “C,  which  is  less  than  that  of  BCGO  with  the  order  of 
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2 

10'  s/cm  (shown  in  figure  1-1).  The  sintered  mixture  did  not  have  as  good  a conductivity 
as  expected.  This  may  be  caused  by  the  additional  unwanted  impurity  phases  such  as 
Ce02,  Ba3(P04)2  and  CeP04.  The  conductivity  of  our  desired  phase  Ba^Ce{PO^  )j  was 
not  clearly  demonstrated  in  this  measurement.  Single  phase  Ba^Ce(PO^)^  should  be 
synthesized  and  electrochemically  characterized. 

3.4  Conclusion 

In  electroless  deposition,  exposure  to  an  acidic  bath  decomposed  part  of  BCGO  to 
form  Ce02  and  phosphorus  impurity  was  absorbed  and  covered  the  surface  of  BCGO 
particles.  They  reacted  at  high  sintering  temperature  to  form  a BaO-rich  compound 
Ba^Ce(PO^)^  and  leave  a certain  Ce02  excess.  Ba^Ce(PO^)^  was  most  likely  the 
desired  tertiary  phase  in  the  composite  BCGSPd. 

Solid-state  reaction  between  BCGO  and  NH4H2PO4  produced  Ba^Ce(PO^)^  and 
Ce02,  which  increased  with  increasing  amount  of  NH4H2PO4.  The  conductivity  of  the 
desired  phase  Ba^Ce(PO^)^  was  not  clearly  demonstrated  in  the  conductivity  results  of 
the  sintered  mixture.  Single  phase  Ba^Ce(PO^)^  will  be  synthesized  and 


electrochemically  characterized  in  the  following  chapter. 
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Figure  3-1.  XRD  patterns  of  composite  BCGSPd  and  BCG35Pd. 
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Figure  3-2.  XRD  results  of  deposited  composites  at  different  pH  conditions. 
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Figure  3-3.  XRD  results  of  composite  BCG4.31Pd  with  different  sintering  temperature 
compared  with  JCPDS  43-0643,  25-0028,  43-1002,  22-0074  and  32-0199. 
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Figure  3-4.  Phase  diagram  of  system  YP04-Ca3(P04)2  (from  ref  [73]). 
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Figure  3-5.  XRD  results  of  sintered  mixtures  of  BCGO  and  NH4H2PO4  with  different 
NH4H2PO4  concentration. 
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Figure  3-6.  XRD  results  of  mixtures  of  BCGO  and  NH4H2PO4  with  or  without  argon 
treatment. 
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Figure  3-7.  IR  spectra  of  composite  BCGSPd,  SrCe^_^Y^O^_g  and  SrZr^_^Yfi^_g  . (a) 
BCGSPd  and  (b)  SrCe^^^Y^O^_g  and  SrZr^_^Y^O^_g  (from  ref  [19]). 
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Figure  3-8.  IR  spectra  of  composite  BCG5Pd  and  Cs2{HSO^){H2PO^) . (a)  IR  spectra  of 
BCG5Pd,  (b)  IR  spectra  of  Cs2{HSO:^){H2PO^)  with  a cubic  structure  at  199  °C 
[58],  and  (c)  vibration  illustration  of  orthophosphate,  diphosphate  and 
hydrogen  phosphate. 


Figure  3-9.  IR  spectra  of  the  sintered  mixture  of  BCGO  and  NH4H2PO4  compared  with 
those  of  composite  BCG5Pd  and  BCG4.31Pd. 
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Figure  3-10.  Impedance  plot  of  the  sintered  mixture  of  BCGO  and  NH4H2PO4  at  800  °C 
in  wet  10%  hydrogen  atmosphere:  (a)  imaginary  vs.  real  impedance  (b) 
imaginary  impedance  vs.  frequency. 
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Figure  3-11.  Arrhenius  plot  of  the  mixture  of  BCGO  and  NH4H2PO4  in  wet  10% 
hydrogen. 
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Table  3-1  Band  assignment  of  SrCe^_^Yfi^_g , SrZr^_^Yfi^_g  , and  BCGSPd. 


Wave  number  (cm‘‘) 

assignment 

SrZr,_Jfi,_,[\9] 

3200 

oOH 

BCGSPd 

3430 

uHjO 

1631 

SH2O 

1035 

^asPOi 

456 

SP02 

V = stretching,  S = bending. 


Table  3-2  Band  assignment  for  CsjiHSO^YHjPO^) . 


Wave  number_(cm'‘) 

Qualitative 

Parent  mode  in  Ta  symmetry 

HighT 

Room  T 

assignment 

495 

474 

Bending 

Uj  (E)  OXO-asymmetric  bend  (IR  inactive) 

513 

496 

modes 

517 

599 

601 

04  (F2)  OXO-symmetric  bend 

620 

623 

641 

670 

968 

971 

u,  (Ai)  XO-symmetric  strech  (IR  inactive) 

1012 

X-OH  stretch 

1109 

1065 

X-0  strech 

U3  (F2)  XO-asymmetric  stretch 

1124 

1091 

1147 

1173 

1286 

X-O-H  in 
plane  bending 
mode 

1715 

1700 

u OH(C) 

2350 

2340 

u OH(B) 

2945 

2760 

V OH(A) 

Data  from  ref 

[58], 
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Table  3-3  Band  assignment  for  BCG4.31Pd  and  sintered  mixture  of  BCGO+NH4H2PO4 


Wave  number  (cm  ') 

assignment 

BCG4.31Pd 

3353 

uHjO 

1060 

1000 

^asP02 

935 

uPOP 

559 

SPO^ 

Sintered  mixture  of 
BCGO+NH4H2PO4 

3455 

vH^O 

1032 

^asP02 

897 

vPOP 

546 

407 

SPO2 

CHAPTER  4 

SYNTHESIS  AND  CHARACTERIZATION  OF  Ba,Ce(PO,), 

4.1  Introduction 

As  mentioned  in  previous  chapters,  Ba^Ce(PO^)^  was  most  likely  the  tertiary 
phase  we  obtained.  The  main  objective  of  this  chapter  is  to  describe  the  synthesis  and 
characterization  of  single-phase  Ba^Ce(PO^)^ . According  to  JCPDS  43-0643, 

Ba^Ce(PO^)^  belongs  to  cubic  system  with  space  group  I43d  (220)  and  lattice 
parameter  a=10.51 12  A.  It  can  be  obtained  by  annealing  a stoichiometric  mixture  of 
BaCeO^ , CeOj  and  NH ^H^PO^  at  1000  °C  for  one  day  and  at  1350  °C  for  10  days  in 
an  open  Pt  crucible  and  then  quenching  in  air;  however,  this  process  is  not  efficient. 

Since  the  solubility  products  of  5(33(^04)2  and  CePO^  are  similar  (3.2x10”^^  and 

lxl0’^\  respectively  [74]),  the  co-precipitation  method  was  attempted.  XRD  and  IR 
spectroscopy  were  used  to  characterize  the  precipitates  with  different  calcination 
temperatures.  XRD  patterns  were  compared  with  that  of  JCPDS  43-0643  to  assure  a 
single-phase  material  was  obtained.  Software  “MS  modeling”  and  “CaRIne 
Crystallography  V3.1”  were  used  to  illustrate  the  crystal  structure  of  Ba^Ce(PO^)^ . 
Raman  and  nuclear  magnetic  resonance  (NMR)  spectroscopy  were  used  to  investigate  the 
local  structures  of  Ba^Ce(PO^)^  under  different  atmospheres. 

As  mentioned  in  chapter  2,  defect  structures  of  doped  phosphates  were  proposed  to 
explain  their  conductivity  behavior.  They  have  limitations  and  only  apply  in  the 


60 


61 


intermediate  and  high  oxygen  partial  pressure  regions.  Nobody  has  studied  the  defect 
structure  in  the  low  oxygen  partial  pressure  region,  especially  in  hydrogen  atmosphere. 
However,  that  pressure  region  is  very  important  for  proton  conductors.  In  this  chapter,  the 
defect  structure  of  Ba^Ce{PO^)y  is  proposed  based  on  the  characterization  results.  The 

concentration  dependence  of  charge  carriers  on  the  oxygen,  hydrogen  or  water  partial 
pressure  was  predicted  according  to  the  defect  structure. 

The  chemical  stability  of  a hydrogen  separation  membrane  material  should  be  good 
in  the  atmosphere  of  its  intended  application.  Therefore,  chemical  stability  of 
Ba^Ce{PO^)^  in  dry  helium  was  tested  by  high  temperature  XRD  (HT-XRD).  XRD 
patterns  of  Ba^Ce(PO^)^  after  separate  treatments  in  hydrogen  and  air  were  compared. 
Finally,  thermal  expansion  behaviors  of  Ba^Ce(PO^)^  in  air  and  hydrogen  atmospheres 
were  investigated  with  a dilatometer. 

4.2  Experiment 

4.2.1  Synthesis 

A solution  of  Ba{NO^)2  and  CeiTVOjlj  -e/ZjO  was  mixed  with  an  appropriate  3:1 
molar  ratio.  The  mixture  was  stirred  and  {NH ^)2HPO^  was  added  drop  by  drop  into  the 
mixture  until  it  was  excess.  The  pH  was  controlled  to  be  certain  Ba^"^  and  Ce^"^  were 
precipitated  in  the  form  of  Ba2{PO^)2  and  CePO^ . The  dried  precipitates  were  calcined  at 
different  temperatures  with  a 3 “C/min  heating  rate  and  4 h dwell  time. 

4.2.2  Characterization 

The  precipitates  with  different  calcination  temperatures  were  characterized  by  XRD 
on  a PHILIPS-APD-3720  using  CuKa  radiation  and  by  IR  on  a INTCOLET  Magna  760 
Bench  spectrometer.  XRD  patterns  were  compared  with  that  of  JCPDS  43-0643.  The 
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synthesized  particles  of  Ba-^CeiPO^)-^  were  suspended  in  acetone  and  a drop  of  the 
ultrasonicated  suspension  was  transferred  to  a carbon  coated  copper  microgrid  for  TEM 
characterization.  Bright  field  images  and  selected  area  diffraction  patterns  of  particles 
were  taken  on  a JEOL200CX.  EMS  online  software  [75]  was  used  to  index  diffraction 
patterns. 

Ba^Ce{PO^)-i  powder  was  pressed  into  pellets,  cold  isostatic  pressed  at  200  MPa  and 
sintered  at  1200  °C  for  8 h.  Pellets  were  ground  into  powder  and  separately  treated  in 
different  atmospheres  (wet  and  dry  air,  nitrogen  and  hydrogen)  at  550  °C  for  24  h.  Raman 
spectra  of  treated  samples  were  measured  using  a Nicolet  MAGNA760  spectrometer 
equipped  with  a multi-channel  CCD  detector. 

For  NMR  experiments,  typical  experimental  parameters  were  as  follows:  a proton 
pulse  duration  of  10.0  ps;  a recycle  delay  of  100  s and  a spectral  width  of  60  Hz. 
Approximately  40-60  mg  of  sample  were  used  and  1024  transients  were  required.  All  'H 
magic  angle  spinning  NMR  spectra  in  this  work  were  acquired  on  a Bruker  Avance 
Spectrospin  400  MHz  spectrometer  operating  at  a frequency  of  400.13  MHz  using  a triple 
resonance  probe  with  a 4.0  mm  rotor.  Spectra  were  acquired  at  a spinning  rate  of  7 kHz 
to  ensure  all  rotational  sidebands  would  fall  outside  of  the  spectral  region  of  interest.  All 
’h  NMR  spectra  were  referenced  to  tetramethylsilane  (TMS).  High-resolution  ^’P  NMR 
was  not  applied  due  to  the  low  coneentration  of  phosphorus  (10"'^  /g)  in  the 
Ba^Ce{PO^)^. 

High  temperature  XRD  was  performed  on  the  Ba^Ce(PO^)^  powders  in  the  dry 
helium.  Ba^Ce{PO^)^  powders  were  separately  treated  in  hydrogen,  argon  and  air 


atmosphere  at  600  °C  for  10  h to  study  its  stability  by  XRD. 
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A Ba^Ce(PO^)^  pellet  was  sintered  at  1200  °C  for  8 h,  polished  to  a 15  mm  length 

and  put  into  a Theta  dilatometer  to  investigate  its  thermal  expansion  behavior.  The 
heating  and  cooling  rates  were  set  at  2 °C/min.  The  gas  flow  rate  was  adjusted  to  maintain 
20  kPa  over-pressure  in  the  chamber. 

4.3  Results  and  Discussion 

4.3.1  Synthesis 

Figure  4-1  shows  XRD  patterns  of  precipitates  with  different  calcination 
temperatures.  All  the  heating  rates  and  dwell  times  were  3 °C/min  and  4 h,  respectively. 
The  precipitate  at  room  temperature  is  a mixture  of  Ba^{PO^)2  and  CePO^ . Peaks  at  20  = 
28“,  38.5  and  43“  belong  to  both  50:3(^04)2  and  CePO^ . The  intensity  of  these  peaks 
decreases  with  increasing  calcination  temperature  and  peaks  disappear  with  an  1100  “C 
calcination  temperature.  A new  peak  appears  at  20  = 27“  with  a calcination  temperature 
of  800  “C  and  it  becomes  strongest  with  calcinations  temperature  1100  “C.  The  XRD 
pattern  of  the  1 100  “C  calcined  sample  successfully  matches  JCPDS  43-0643  (shown  in 
table  4-1)  except  a d spacing  of  1.9128  A.  However,  d spacing  of  1.9128  can  be  still  be 
indexed  as  (521),  as  shown  in  figure  4-1.  No  second  phase  is  detected.  The  difference 
between  the  experimentally  obtained  d spacings  and  those  in  JCPDS  43-0643  is  small. 
This  difference  may  be  caused  by  a different  synthesis  method. 

Figure  4-2  shows  the  calculated  XRD  profile  of  Ba^CeiPO^)^  using  the 
crystallographic  data,  lattice  parameter  and  CaRIne  Crystallography  V3.1  software.  The 
calculated  XRD  files  predict  all  of  the  d spacings  in  our  experimental  XRD  result.  It 
includes  reflection  (521)  at  20=  47.5“,  which  is  not  in  the  JCPDS  43-0643  but  appears  in 
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our  experimental  XRD  data.  It  suggests  single-phase  Ba^Ce(PO^)^  is  produced  during 
calcination  at  1100  °C. 

Infrared  spectra  of  precipitates  with  different  calcination  temperatures  are  shown  in 
Figure  4-3.  Table  4-2  gives  assignments  of  bands  with  different  calcination  temperatures. 
The  bands  at  1050  and  566  cm  ' of  Ba^Ce(PO^)^  belong  to  a stretching  mode  and 

bending  mode  of  POl~ , respectively.  The  band  at  2000  cm''  may  originate  from  the 
vibration  mode  of  bonds  M-O(P)  (M=  Ba,  Ce).  The  band  of  stretching  vibration  at  1100 
°C  is  broad,  which  can  be  narrowed  by  measuring  IR  at  lower  temperature.  It  may 
obscure  the  possible  splitting.  The  bending  vibration  modes  of  POl~  become  degenerated 
with  increasing  calcination  temperatures.  Compared  with  BCG5Pd  in  figure  4-4,  the 
stretching  and  bending  vibration  modes  of  synthesized  Ba-^Ce(PO^)^  move  to  slightly 
higher  wave  number.  However,  it  still  suggests  that  a similar  chemical  environment  of 
the  tertiary  phase  and  our  synthesized  Ba^Ce{PO^  )j . 

4.3.2  Particle  Size  Characterization 

Figure  4-5  (a)  is  a bright  field  image  of  Ba^Ce{PO^)^  particles,  which  revealed  a 
homogeneous  and  well-crystallized  Ba^Ce{PO^)^  with  a grain  size  around  several 
hundred  nanometers.  Figure  4-5  (b)  is  a selected  area  diffraction  pattern  of  Ba^Ce(PO^)^ 

with  incidence  [B9].  TEM  image  simulation  software  was  used  to  obtain  corresponding 
d-spacings,  which  are  listed  in  table  4-1.  They  match  those  obtained  by  XRD. 

4.3.3  Crystal  Structure  Characterization 

According  to  JCPDS  43-0643,  Ba^Ce(PO^)^  belongs  to  cubic  system  with  space 

group  1 43d  (220)  with  lattice  parameter  a=10.51 12  A.  According  to  Wyckoff  positions 
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of  group  220,  the  P atoms  are  fixed  at  12a  sites,  the  Ba  and  Ce  atoms  are  randomly 
distributed  on  the  16c  sites,  and  the  O atoms  are  at  48e  sites.  Without  neutron  diffraction 
experiments,  the  crystal  structures  of  Ba^Ce(PO^)^  can  only  be  approximately  illustrated 
based  on  the  available  information  of  the  Ba^La{PO^)^  [76],  which  has  a similar  crystal 
structure  with  Ba^Ce(PO^)^ . The  coordinates  of  atoms  in  Ba^Ce(PO^)^  are  based  on 
those  of  Ba^La{PO^);i . Software  “MS  modeling”  and  “CaRIne  Crystallography  V3.1”  are 
used  to  illustrate  the  crystal  structure  of  Ba^Ce(PO^)^ . The  results  are  shown  in  figure  4- 
6 and  4-7,  respectively. 

The  O atoms  have  two  sets  of  different  crystallographic  positions  in  the  eulytite- 
type  Ba^Ce(PO^)^ . However,  software  “MS  modeling”  could  not  plot  the  two  positions 
at  the  same  time.  Only  the  first  set  of  oxygen  sites  is  considered  when  using  the  ‘MS 
modeling”.  The  coordinates  of  P,  O and  metals  M (M=  Ba,  Ce)  in  Ba^Ce{PO^)^  are 
(0.375,  0,  0.25),  (0.466317,  0.361411,  0.707516),  and  (0.06025,  0.06025,  0.06025), 
respectively,  based  on  that  of  Ba^La(PO^)j  [76].  The  bond  lengths  and  angles  are 
measured  from  the  illustrated  structure.  The  measured  bond  lengths  of  M-0  and  P-0  in 
Ba^CeiPO^)^  are  2.818  and  1.427  A,  respectively.  The  angle  of  O-P-0  is  104.88“.  The 
[111]  directional  view  clearly  shows  the  metal  atoms  are  in  the  three-fold  axis. 

Figure  4-7  (a)  shows  the  crystal  structure  of  Ba^Ce(PO^)^  using  “CaRIne 

Crystallography  V3.1”  software.  Ba  atom  is  represented  by  blue  sphere.  Ce  atom  is 
represented  by  blue  sphere  with  red  in  the  center.  The  positions  of  Ba  and  Ce  are  random 
but  their  molar  ratio  of  3 is  fixed.  The  two  different  crystallographic  positions  for  the  O- 
atoms,  01  and  02,  are  represented  by  dark  and  shallow  green  spheres.  Figure  4-7  (b) 
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gives  [001]  directional  views  of  Ba^La{PO^)^ , clearly  showing  the  two  sets  of  oxygen 
sites  [76].  The  two  possible  positions  of  oxygen  atoms  cause  rotational  disorder.  The 
rotational  disorder  usually  changes  the  coordination  of  the  metal  atoms  from  6 to  9. 
Ba^Ce{PO^)-i  belongs  to  the  same  eulytite  structure  and  a rotational  disorder  of  the 
phosphate  group  is  also  expected  as  Ba^La{PO^)^ . Neutron  diffraction  data  are  needed  for 
further  discussion. 

Figure  4-8  (a)  and  (b)  gives  [010]  and  [1 1 1]  directional  views  of  Ba^Ce(PO^)^, 
respectively,  with  all  the  oxygen  atoms  in  the  same  crystallographic  position  02.  There  is 
a three-fold  channel,  which  may  be  good  for  the  proton  transfer. 

As  mentioned  in  the  beginning,  these  illustrated  crystal  structures  of  Ba^Ce(PO^)^ 
are  based  on  the  available  information  of  the  Ba^La{PO^)j . Therefore,  all  the  obtained 
information  is  just  approximation  for  the  real  structure.  Neutron  diffraction  experiments 
would  be  desirable  and  extremely  useful  in  a better  understanding  of  the  crystal  structure. 
4.3.4  Defect  Structure  Characterization 
4.3.4.1  Raman  characterization 

Raman  spectroscopy  was  used  to  study  the  proton  incorporation  and  oxygen 
vacancy  formation  in  Ba^Ce(PO^)^  under  dry  or  wet  atmospheres. 

Raman  Characterization  of  Dry  Atmosphere  Treated  Samples 

Figure  4-9  shows  the  Raman  spectra  of  Ba^Ce(PO^)^  after  separate  treatments 

with  dry  air,  nitrogen  and  hydrogen.  Raman  shifts  of  Ba^Ce{PO^)^  are  listed  in  table  4-3 
and  compared  with  Sr  doped  LaP04  after  dry  atmosphere  treatments  [42]. 
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PO4  ■ is  the  main  existing  form  in  Ba^Ce{PO^  under  dry  conditions.  It  has  four 

vibration  bands:  an  asymmetric  PO2  stretching  mode  at  1057  and  1020  cm'\  a symmetric 
stretching  mode  at  944  cm'\  the  O-P-0  p2-bending  mode  at  571  cm'^  and  the  O-P-0  E- 
bending  mode  at  463  and  439  cm'\  The  P04^'  Raman  shift  in  Ba^Ce(PO^)^  is  slightly 
different  from  that  in  Sr  doped  LaP04,  which  is  due  to  different  cations  in  the  lattice  [40], 
Raman  shifts  at  1057  and  1020  cm"'  in  the  dry  hydrogen  treated  Ba^Ce(PO^)^  can 

also  be  assigned  to  ^2^7  according  to  reference  [41],  Their  relative  intensities  are 
highest  in  dry  hydrogen  treated  sample  and  lowest  in  dry  air  treated  one,  which  is  shown 
in  figure  4-9  (b).  It  may  suggest  that  the  oxygen  vacancy  concentration  increases  with 
reduced  oxygen  partial  pressure  under  the  dry  condition. 

For  the  dry  hydrogen  treated  Ba^Ce{PO^  )j , proton  incorporation  is  expected. 

However,  Raman  did  not  detect  any  HPOl~  shift,  which  is  reported  to  be  strong  at  the 
2500-2150  and  1150-1000  cm'*  regions  [41]. 

Raman  Characterization  of  Wet  Atmosphere  Treated  Samples 

Raman  was  also  used  to  investigate  the  proton  incorporation  in  the  presence  of 
water  vapor.  Vibration  bands  in  the  regions  of  1240-1460  and  2330-2900cm'*  are 
reported  as  the  existence  of  POH  bonds  in  reference  [77].  Vibration  bands  at  2176-2787 
and  1 278-1 332cm‘*  are  reported  in  reference  [78]  as  the  stretching  and  bending  modes  of 
POH  in  acid  pyrophosphate  KMgQ^H 2P2O-1  H^O,  respectively. 

Figure  4-10  (a,  b)  is  Raman  spectra  of  Ba^CeiPOi^')^  after  separate  treatments  in 
wet  and  dry  nitrogen.  Table  4-4  lists  Raman  shifts  of  Ba^Ce^POi^^^,  compared  with  Sr 
doped  LaP04  treated  in  the  wet  atmosphere  [42].  Some  liquid  water  was  added  into 
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powder  Ba^Ce(PO^)^  to  show  the  stretching  mode  of  absorbed  water,  which  appeared  at 
3333  cm’’.  This  stretching  mode  of  water  molecule  did  not  appear  in  the  wet  nitrogen 
treated  Ba^Ce(PO^)^,  suggesting  that  no  water  molecules  are  inserted  in  the  lattice.  The 

bands  at  2881  and  2848  cm’’  are  assigned  to  the  stretching  modes  of  POH  in  the  HPO4 
group.  The  bands  at  1439,  1295  and  1131  cm’’  are  assigned  to  the  bending  modes  of 
POH.  They  suggest  that  protons  exist  in  the  form  of  {HPO^)po^  in  the  wet  nitrogen  treated 

Ba^Ce{PO^)^ . 

Raman  results  also  provide  some  information  about  oxygen  vacancy  concentration 
in  the  presence  of  water  vapor.  Figure  4-10  (c)  demonstrates  that  the  relative  intensities 
of  P20j~  vibration  modes  at  1057  and  1020  cm’’  decreased  in  wet  nitrogen.  This 
indicates  that  water  partial  pressure  suppresses  the  oxygen  vacancy  formation. 

Figure  4-1 1 is  the  Raman  spectra  of  Ba^CeiPO^)^  with  wet  hydrogen  treatment. 

The  broad  band  at  1350  cm’’  is  assigned  to  the  bending  mode  of  POH  bonds  according  to 
references  [78,  79].  Its  relative  intensity  decreases  after  drying  at  100  “C  for  1 h, 
suggesting  the  proton  concentration  decreases  with  drying  process.  It  is  very  difficult  to 
compare  the  relative  intensity  of  vibration  bands  at  1057  and  1020  cm’’  due  to  the 
overlap  effect  with  band  at  1350  cm’’.  Therefore,  it  is  unclear  about  the  water  vapor 
effect  on  the  oxygen  vacancy  concentration  in  hydrogen  treated  BajCeiPO^)^ . 

Raman  Characterization  Summary 

Raman  results  explain  the  proton  incorporation  in  the  presence  of  water  vapor. 
Without  free  hydrogen  gas,  in  the  presence  of  water  vapor,  bound  hydrogen  is  produced 
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and  oxidized  forming  protons.  The  formed  protons  attach  onto  the  O atoms  of  PO4  group 
and  are  incorporated  into  the  lattice. 

The  band  wave  numbers  of  incorporated  protons  under  wet  atmospheres  can  be 
summarized  as  follows:  2881,  2848,  1439,1295,1 135  cm"'  in  wet  nitrogen,  1350  cm’'  in 
wet  hydrogen.  The  wave  number  reflects  the  bond  strength.  Thus,  the  bond  strength  of 
POH  decreased  with  decreasing  oxygen  partial  pressure  under  wet  condition. 

Oxygen  vacancy  concentration  increases  with  reduced  oxygen  partial  pressure 
under  the  dry  condition.  Water  vapor  suppresses  the  oxygen  vacancy  formation  in 
nitrogen.  The  water  vapor  effect  on  the  oxygen  vacancy  concentration  in  hydrogen  is  not 
clearly  demonstrated  by  Raman  results. 

Unfortunately,  Raman  results  did  not  demonstrate  the  existence  of  protons  in  the 
dry  hydrogen  treated  Ba^Ce(PO^)^ . 

4.3.4.2  NMR  characterization 

According  to  reference  [42],  in  Sr  doped  LaP04,  chemical  shifts  at  5.4  and  1.2  ppm 
in  high-resolution  NMR  were  considered  due  to  adsorbed  water  molecules.  The 
chemical  shift  at  12.4  ppm  was  considered  due  to  dissolved  protons  and  it  had  a slower 
relaxation  time  than  that  at  5.4  and  1.2  ppm,  reflecting  that  the  dipole  interaction  between 
protons  was  much  weaker  than  that  between  protons  in  adsorbed  water  molecules. 

Figure  4-12  is  *H  NMR  results  of  Ba^Ce(PO^)^  after  different  atmosphere 
treatments.  Corresponding  chemical  shifts  of  'H  NMR  are  listed  in  table  4-5  and 
compared  with  Sr  doped  LaP04.  Dry  air  treated  Ba^Ce(PO^)^  has  a chemical  shift  at  5 

ppm,  which  is  assigned  to  the  attached  water  molecules  from  the  environment.  In 
addition  to  5 ppm,  the  chemical  shift  at  15  ppm  appears  in  the  dry  hydrogen  treated 
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Ba^Ce(PO^)^ . The  relative  intensity  of  the  chemical  shift  at  15  ppm  is  much  larger  than 

at  5 ppm.  This  clearly  demonstrates  proton  insertion  under  the  dry  hydrogen  condition 
despite  the  fact  that  it  is  not  observed  in  the  Raman  spectra.  The  inserted  protons  can  be 
expressed  in  the  form  of  (HPO^)pQ^  [42].  The  difference  between  the  chemical  shift  at  15 
ppm  in  Ba^Ce{PO^)^  and  12.4  ppm  in  Sr  doped  LaP04  is  caused  by  the  different 
chemical  environment. 

The  'H  NMR  chemical  shift  at  15  ppm  in  wet  nitrogen  treated  Ba^Ce(PO^)^, 

shown  in  figure  4-12,  suggests  again  the  proton  insertion  in  the  presence  of  water  vapor. 

Raman  spectra  did  not  indicate  whether  or  not  oxygen  vacancy  existed  in  the  lattice 
under  the  wet  hydrogen  atmosphere.  ^'P  NMR  can  illustrate  it  more  clearly.  The  chemical 
shift  of  orthophosphate  is  -^A,  pyrophosphate  is  -7.2  and  -8.4,  and  hydrogen  phosphate 
is  -30.6  ppm  according  to  reference  [42].  It  would  have  been  very  helpful  if  the  ^'P  NMR 
was  applied.  However,  it  was  not  applied  due  to  the  low  concentration  of  phosphorus  (10' 
/g)  in  the  Ba^Ce(PO^)^ . 

4.3.4.3  Defect  structure 

Based  on  the  local  structure  characterization  results,  a defect  structure  model  of 
Ba^Ce{PO^)^  was  proposed  in  low  ( pq^  < atm ),  low-intermediate 

(10’’^  < pq^  <\Q~'’ atm),  intermediate  (10“’  < pq^  <l0“’a?w)  and  high  >l0“‘a?w) 

oxygen  partial  pressure  regions.  Defect  concentration  dependencies  of  P02,  Ph2  and  Ph20 
were  derived  and  listed  in  table  4-6. 

Low  Oxygen  Partial  Pressure  Region 

Oxygen  vacancy  exists  in  the  phosphate  lattice  at  the  low  oxygen  partial  pressure 
region  and  plays  a similar  role  as  that  in  perovskites.  However,  it  exists  as  a different 
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form  in  phosphates.  Oxygen  vacancy  in  phosphates  can  be  considered  as  two  PO4  sites 
sharing  two  P and  seven  O atoms.  It  also  can  be  considered  as  two  PO4  sites  losing  one  O 
atom.  Oxygen  vacancy  is  expressed  as  {PiO-jyjpo^  in  this  work,  the  same  way  Amezawa 
did  in  his  paper  [42]. 

The  following  defect  equilibrium  always  exists  at  the  low  oxygen  partial  pressure 
region: 

->iP20y)2Po,  +2e'  . K1  (4-1) 

i 

^\-[iP2^7)2PO,]^^ PO2  ■ 

The  internal  electron-hole  pair  formation  follows  the  reaction: 

null  = e +h  . K2  (4-2) 

At  the  low  oxygen  partial  pressure  region,  hole  concentration  is  negligible  and  its 
contribution  to  electroneutrality  can  be  omitted. 

Four  cases  were  discussed  in  the  low  oxygen  partial  pressure  region.  Case  1 is  the 
low  water  partial  pressure  condition  without  free  hydrogen,  case  2 is  the  low  water  partial 
pressure  condition  with  free  hydrogen,  case  3 is  the  high  water  partial  pressure  condition 
without  free  hydrogen,  and  case  4 is  the  high  water  partial  pressure  condition  with  free 
hydrogen. 

Case  1:  low  water  partial  pressure  without  free  hydrosen. 

Under  low  water  partial  pressure  condition  without  free  hydrogen,  there  is  no 
proton  source  for  Ba^Ce(PO^)^ . Oxygen  vacancies  are  the  only  significant  positive 
defects  to  charge  balance  electrons.  The  electroneutrality  condition  is  simplified  as 


2[{P  20j)2POf]  — • 
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Therefore, 

1 

K,=  [nl2]n^p\. 

_i 

Electron  concentration  and  oxygen  vacancy  concentration  are  proportional  to  . 

There  is  no  water  partial  pressure  dependence  due  to  the  low  water  partial  pressure 
condition.  Without  a proton  source,  there  is  also  no  hydrogen  partial  pressure 
dependence. 

Case  2:  low  water  partial  pressure  with  free  hydrosen. 

Free  hydrogen  is  the  only  proton  source  in  case  2.  Free  hydrogen  can  be  oxidized 
forming  protons,  which  are  attached  onto  the  O atoms  in  the  PO4  tetrahedron  and 
incorporated  into  the  lattice.  The  corresponding  defect  equilibrium  is  written  as 

{PO.Ypo,  +1^2^  iHPO,)po,  + e K3  (4-3) 

K,=n[{HPO,)Po,]p2- 

A small  amount  of  water  vapor  possibly  exists  due  to  the  existence  of  hydrogen  and 
oxygen  formed  in  formula  (4-1): 

H2+^02^H^0  K4  (4-4) 

^4  = PH20PH2P0I  ■ 

Dry  hydrogen  belongs  to  low  water  partial  pressure  condition  with  free  hydrogen. 
For  dry  hydrogen  treated  Ba^CeiPO^)^,  a chemical  shift  at  15  ppm  in  'H  NMR 

spectroscopy  result  (shown  in  figure  4-12)  clearly  demonstrates  proton  incorporation. 
Figure  4-9  (b)  demonstrates  that  the  oxygen  vacancy  concentration  increases  with 
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reduced  oxygen  partial  pressure  under  dry  condition.  Therefore,  oxygen  vacancies  and 
incorporated  protons  are  the  positive  defects  in  the  lattice.  The  electroneutrality  condition 
is  written  as 

[(//P04)poJ  + 2[(P2C>7)/.oJ  = «- 

If  oxygen  vacancy  concentration  dominates, 

2[(/’2  07)roJ  = « 

2K,p~^^=2K,K,p„\oPh^  =«'• 

-i  i _i 

Therefore,  electron  and  oxygen  concentration  are  proportional  to  and 

They  are  proportional  to  if  water  partial  pressure  is  held  constant  in  experiments. 

If  proton  concentration  dominates, 

[(HPO,)Po^  = n 

_i  i _i  i 
^3^4^  Ph^oPq2  ~ ^^P  H2~  ^ 

-i  i i 

Therefore,  electron  and  proton  concentration  are  proportional  to  p^^  p^^^  and  . 

They  are  proportional  to  if  water  partial  pressure  is  held  constant  in  experiments. 

Case  3:  hieh  water  partial  pressure  without  free  hydroeen. 

Without  free  hydrogen,  water  vapor  is  the  only  proton  source.  Water  vapor  reacted 
with  oxygen  vacancies  in  perovskites  to  form  protons.  It  also  could  react  with 
pyrophosphates  to  form  hydrogen  phosphates  in  Ba^Ce{PO^)^ . Defect  equilibrium  can 


be  written  as 


74 


{P20^)jpo,+H^O{g)^2{HPO,)po^.  K5  (4-5) 

Ks=kIk;'k;' 

Water  vapor  also  could  first  produce  bound  hydrogen  by  reaction  (4-4).  Produced 
bound  hydrogen  is  then  oxidized,  forming  protons  and  attaching  to  the  O atoms  of  PO4 
group.  The  whole  process  includes  defect  equilibrium  (4-3)  and  (4-4).  Defect  equilibrium 
is  simplified  as 

Vpo,  +H^0-^  1{HPO,)po^  + 2e'  + 1 02 . K6  (4-6) 

k,=kIk;^  =k,k;^ 

The  electroneutrality  condition  is  written  as 

\.{HPO^^ ) pQ^  ] + 2[(/’2  0-j  )2POt  ] — ^ • 

Therefore, 

i -i 

Kyti~'pjj  +2Kpi~^  = «. 

_i  -i  ' i 

p ^ = n . 

4K^ 

It  is  difficult  to  predict  defect  concentration  dependency  of  oxygen  partial  pressure 
from  the  above  formula  without  the  information  about  equilibrium  constants  AT,  - . 

Case  3 is  simplified  into  two  conditions:  one  is  oxygen  vacancy  concentration 
dominance  condition  and  the  other  is  proton  concentration  dominance  condition. 

If  oxygen  vacancy  concentration  still  dominates  in  high  water  partial  pressure 
condition,  electroneutrality  condition  is  simplified  as 
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2[(/*2^7)2f>04  ] ~ ” • 


2K,p^l  =2K,K,p-„\oP„^  =n^ 


Proton  concentration  is  negligible  under  this  hypothesis.  The  electron  concentration 


and  oxygen  vacancy  concentration  are  proportional  to  p^^  and  p^^p^^  ■ 

If  high  water  partial  pressure  suppresses  the  oxygen  vacancy  formation  and  proton 
concentration  dominates,  the  electroneutrality  condition  is  written  as 

[{HPO,)po^]  = n. 

Therefore, 


The  protons  and  electron  concentration  are  proportional  to  pj"  pfj  ^ . 

Mostly  likely,  case  3 is  the  mixture  of  these  two  conditions. 

Case  4:  hish  water  partial  pressure  with  free  hvdrosen. 

Under  high  water  partial  pressure  condition  with  free  hydrogen,  both  water  vapor 
and  free  hydrogen  provide  proton  sources.  The  corresponding  defect  equilibriums  are 
written  as 


{PO.Ypo,  ^{HPO,)po,  +e 


K3 


(4-3) 


2(^04  Ypo,  + "2O  ^ 2{HP0,  )po,  + 2e  +-O2  . 


K6 


(4-6) 


The  electroneutrality  condition  can  be  written  as 


[{HPO,ypoY  + ^\iPiOi)2PoY  = ri  . 
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The  only  difference  between  case  3 and  case  4 is  that  free  hydrogen  is  involved  in 
the  proton  incorporation  in  case  4.  If  we  know  about  the  information  of  equilibrium 
constants  K^  - under  these  two  cases,  the  free  hydrogen’s  contribution  could  easily  be 
demonstrated.  Unfortunately,  no  such  fundamental  research  was  reported. 

Wet  hydrogen  belongs  to  case  4.  For  wet  hydrogen  treated  Ba^Ce(PO^)j,  Raman 

spectroscopy  result  in  figure  4-1 1 shows  a broad  band  at  1350  cm'’,  assigned  to  the 
bending  mode  of  POH  bond,  which  clearly  demonstrates  the  proton  incorporation.  Figure 
4-11  also  shows  that  POH  vibration  intensity  is  strong  in  wet  hydrogen  treated  sample 
while  it  is  negligible  in  dry  hydrogen  treated  one.  This  indicates  that  water  vapor  plays  a 
much  more  important  role  than  free  hydrogen  in  proton  incorporation.  Therefore,  without 
further  information  about  free  hydrogen  contribution,  all  the  derivations  in  case  3 can 
approximately  apply  in  case  4.  However,  the  proton  concentration  is  expected  to  be 
higher  than  in  case  3 due  to  the  free  hydrogen  contribution.  Therefore,  proton 
concentration  is  considered  to  be  dominant  in  case  4.  The  electroneutrality  condition  is 
written  as 

[(HPO,)po,]  = n. 

_i  i _i  i 

I 

The  electron  concentration  and  proton  concentration  are  proportional  to  . The 

_\_  \_ 

defect  concentration  dependence  is  proportional  to  p^^^ . 
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Low-Intermediate  Oxygen  Partial  Pressure  Region 

Oxygen  vacancy  also  exists  in  the  phosphate  lattice  at  the  low-intermediate  oxygen 
partial  pressure  region  and  plays  a similar  role  as  that  in  the  low  oxygen  partial  pressure 
region.  The  following  defect  equilibrium  also  exists: 

^^PO,Ypo,^{P20^)^po,+^e  +^^0^.  K1  (4-1) 

However,  oxygen  vacancy  concentration  is  smaller  at  this  low-intermediate  oxygen 
partial  pressure  region.  Hole  concentration  can  still  be  omitted.  The  total  defect 
concentration  is  lower  than  that  in  the  low  oxygen  partial  pressure  region. 

With  a small  trace  of  free  hydrogen,  the  oxygen  partial  pressure  will  decrease  to 
Pq^  < 10’'^  atm . Therefore,  case  2 (low  water  partial  pressure  condition  with  free 

hydrogen)  did  not  exist  in  the  low-intermediate  region.  Only  cases  1,  3 and  4 were 
discussed  in  this  region.  Case  1 is  the  low  water  partial  pressure  condition  without  free 
hydrogen,  case  3 is  the  high  water  partial  pressure  condition  without  free  hydrogen,  and 
case  4 is  the  high  water  partial  pressure  condition  with  free  hydrogen. 

Case  1:  low  water  partial  pressure  without  free  hydrosen. 

Under  low  water  partial  pressure  condition  without  free  hydrogen,  there  is  no 
proton  source  for  Ba^Ce(PO^)^.  Oxygen  vacancies  are  the  only  positive  defects  to 
charge  balance  electrons.  Therefore,  the  electroneutrality  condition  is  simplified  as 

2[(^2^7)2/>04  ] = ” • 

Therefore, 

I 

K,=  [n/2yp\. 
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The  electron  concentration  and  oxygen  vacancy  concentration  are  proportional  to 

pj . There  is  no  water  partial  pressure  dependence  due  to  low  water  partial  pressure 

condition.  Without  a proton  source,  there  is  also  no  hydrogen  partial  pressure 
dependence. 

Case  3:  hish  water  partial  pressure  without  free  hydrosen. 

Without  free  hydrogen,  water  vapor  is  the  only  proton  source.  Water  vapor  could 
react  with  oxygen  vacancies  to  form  protons.  It  also  could  first  produce  bound  hydrogen. 
Produced  bound  hydrogen  is  then  oxidized,  forming  protons  and  attaching  to  the  O atoms 
of  PO4  group.  The  corresponding  defect  equilibriums  can  be  written  as 

iPiOi)2Po,  +H^O{g)  ->  2{HPO,)po,  ■ K5  (4-5) 

^{PO.Ypo,  +H2O  ^ 2{HPO,)po,  +2e’  +^0^ . K6  (4-6) 

The  electroneutrality  condition  is  written  as 
{{HPO^) PQ^  ] + 21(^2 0-j )2POf  1 = ” • 

There  still  exist  two  possibilities:  one  is  proton  concentration  dominance  and  the 
other  is  oxygen  vacancy  concentration  dominance. 

If  oxygen  vacancy  concentration  dominates  in  high  water  partial  pressure 
condition,  then  electroneutrality  condition  is  simplified  as 
2[(P2  6^7)2^04 1 “ ” • 

2K,p2  =2K,K,pjj^oPH,  =«' 

Proton  concentration  is  negligible  under  this  hypothesis.  The  electron  concentration 

_i  _i  1 

and  oxygen  vacancy  concentration  are  proportional  to  p^^  and  PifpP^^  • 
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If  proton  concentration  dominates,  the  electroneutrality  condition  is  written  as 
[{HPO,)po,^  = n. 

_i  i -i  i 

_i  i 

The  proton  and  electron  concentration  are  proportional  to  . 

Without  free  hydrogen,  the  defect  concentration  dependence  of  hydrogen  partial 
pressure  was  not  further  derived  in  this  case. 

Case  4:  hish  water  partial  pressure  with  free  hydrosen. 

Under  high  water  partial  pressure  condition  with  free  hydrogen,  both  water  vapor 
and  free  hydrogen  provide  proton  sources.  Their  corresponding  defect  equilibriums  are 
written  as 

iP04)po,  {HPO,)po,  + e K3  (4-3) 

2{PO^Ypo^  +H^0^2{HPO^)po^  +2e  . K6  (4-6) 

The  electroneutrality  condition  can  be  written  as 

[{HPO,)po^]  + 2[{P^O,)^PO^]  = n . 

The  only  difference  between  case  3 and  case  4 is  that  free  hydrogen  is  involved  in 
the  proton  incorporation  in  case  4.  However,  the  free  hydrogen  concentration  is  not  high 
at  the  low-intermediate  oxygen  partial  pressure  region.  Its  contribution  can  be  omitted 
compared  with  water  vapor  contribution.  All  the  derivations  in  case  3 can  approximately 
apply  in  case  4. 

If  oxygen  vacancy  concentration  dominates,  the  electron  concentration  and  oxygen 

-1  i -i 

vacancy  concentration  are  proportional  to  and  pj/^Pff\o- 
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If  proton  concentration  dominates,  the  proton  and  electron  concentration  are 

i _i  1 

proportional  to  and  Po\Ph^o- 

Intermediate  Oxygen  Partial  Pressure  Region 

There  is  no  free  hydrogen  in  the  intermediate  oxygen  partial  pressure  region  and 
water  vapor  is  the  only  possible  proton  source.  Two  cases  were  discussed  in  this  region. 
One  is  the  low  water  partial  pressure  condition  and  the  other  is  high  water  partial  pressure 
condition. 

Casel:  low  water  partial  pressure 

Dry  nitrogen  gas  belongs  to  case  1 . Raman  spectroscopy  result  in  figure  4-9 
suggests  that  oxygen  vacancy  exists  in  the  dry  nitrogen  treated  Ba^Ce(PO^)^ . However, 

oxygen  vacancy  concentration  is  lower  in  dry  nitrogen  than  in  dry  hydrogen. 

The  corresponding  defect  equilibrium  in  the  dry  atmosphere  is  written  as 

2{PO,Ypo4  ^ {PiOYiPo,  +1^2  +2e'  K1  (4-1) 

The  electroneutrality  condition  is  written  as 

2[(P  2^l)lPO,^'^  P~^  ■ 

Despite  the  fact  that  the  electron  and  hole  concentrations  are  not  high  due  to  the 
phosphate  structure  and  the  intermediate  oxygen  partial  pressure;  however,  they  could  not 
be  omitted.  The  electron  concentration  restricts  the  oxygen  vacancy  and  hole 
concentration.  Therefore,  the  total  defect  concentration  in  this  region  is  expected  to  be 
low.  There  is  no  apparent  oxygen  partial  pressure  dependency  under  this  condition. 
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Case  2:  hish  water  partial  pressure. 

Wet  nitrogen  belongs  to  case  2.  Water  vapor  is  the  only  proton  source  in  wet 
nitrogen.  Raman  spectroscopy  in  figure  4-10  provides  some  information  about  the  proton 
incorporation  under  this  condition.  The  stretching  mode  of  absorbed  water  molecule  at 
3333cm''  did  not  appear  in  the  wet  nitrogen  treated  Ba^Ce(PO^)^,  suggesting  that  no 

water  molecules  were  inserted  in  the  lattice.  The  stretching  modes  at  2881,  2848cm''  and 
bending  modes  at  1439,  1295,  1 131cm  ' of  POH  in  figure  4-10  indicated  that  protons  are 
inserted  in  the  lattice.  This  indicates  that  water  vapor  produces  bound  hydrogen  and  the 
produced  bound  hydrogen  is  then  oxidized,  forming  protons.  The  protons  are  attached  to 
O atoms  of  PO4  group  and  incorporated  into  the  lattice. 

The  corresponding  defect  equilibrium  in  wet  atmosphere  is  written  as 

2{PO^)po^+H20^2iHPO^)po^+2e' +^0:^.  K6  (4-6) 

The  electroneutrality  condition  is  written  as 

[iHPO,)poJ  + 2[iP2 O2 )2Po,  ] + P = n. 

Based  on  defect  equilibrium  (4-6),  proton  concentration  and  electron  concentration 
increase  in  the  presence  of  water  vapor.  However,  electron  concentration  is  not  high  and 
almost  equals  hole  concentration  in  the  intermediate  oxygen  partial  pressure  region, 
which  restricts  the  proton  concentration  based  on  the  electroneutrality  condition. 
Therefore,  the  total  defect  concentration  will  not  change  much  in  the  presence  of  water 
vapor  in  the  intermediate  region.  Their  dependence  of  oxygen  partial  or  water  partial 
pressure  will  not  be  apparent  in  this  region. 
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High  Oxygen  Partial  Pressure  Region 

We  do  not  know  of  any  appropriate  negative  ionic  species.  Therefore,  electrons  are 
the  only  negative  defects  to  charge  balance  holes,  oxygen  vacancies  and  possible  protons. 
Oxygen  vacancy  concentration  could  be  omitted  in  this  region.  In  the  presence  of  water 
vapor,  equal  amounts  of  protons  and  electrons  were  incorporated  into  the  lattice  based  on 
reaction  (4-6).  However,  electron  concentration  is  usually  lower  than  hole  concentration 
under  high  oxygen  partial  pressure  region.  Extra  negative  defects  are  needed  to  charge 
balance  holes. 

The  following  defect  equilibriums  can  exist  in  the  high  oxygen  partial  pressure 
region: 

^PO.Ypo,  ^{PiOYiPO,  +2e'  . K1  (4-1) 

'^{PO^Ypo,  +f^20  ^ 2{HPO^)po^  +2e  . K6  (4-6) 

The  electroneutrality  condition  under  dry  atmosphere  is  written  as 

2[(P 2 ^7  ) 2PO,  ] + P ~ ” ■ 

The  electroneutrality  condition  under  wet  atmosphere  is  written  as 

[{HPO^) pQ^  ] + 2[(P2 0-j)2po^  ] + P = • 

However,  as  stated  above,  electrons  should  not  be  a primary  species  at  high  oxygen 
partial  pressure.  Therefore,  the  defect  structure  model  in  this  region  faces  a challenge  and 
more  studies  are  needed.  The  defect  concentration  dependencies  of  oxygen  partial 
pressure  are  not  further  derived  due  to  the  challenge  we  pointed  out. 
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4.3.5  Chemical  and  Thermal  Stability 

Figure  4-13  is  the  high  temperature  XRD  of  Ba^Ce(PO^)^  in  the  helium 

atmosphere.  The  peak  at  20=48”  came  from  the  result  of  the  silver  paste  used  as  an 
adhesive.  All  the  peak  positions  moved  to  lower  angles  with  increasing  temperature, 
suggesting  the  corresponding  lattice  expansion  during  heating.  No  second  phase  was 
detected,  suggesting  that  the  stability  of  Ba^Ce{PO^  )j  is  good  in  helium  up  to  600  ”C. 

The  XRD  in  Figure  4-14  shows  little  difference  between  Ba^Ce(PO^)^  samples 
after  separate  treatments  in  argon,  air,  and  hydrogen  atmosphere  at  600  ”C  for  10  h.  This 
can  be  regarded  as  an  indication  of  stability  of  Ba^Ce(PO^)^  in  air,  argon,  and  hydrogen. 

4.3.6  Thermal  Expansion 

The  thermal  expansion  behavior  of  Ba^Ce{PO^)j,  in  air  and  hydrogen  is  shown  in 
Figure  4-15.  The  thermal  expansion  coefficients  of  Ba^Ce(PO^)^  in  the  hydrogen  and  air 

atmospheres  were  18x10’^  /°C  and  15x10'^  /°C,  respectively. 

4.4  Conclusion 

Single-phase  Ba^CeiPO^)^  was  synthesized  successfully  by  a coprecipitation 
method.  Ba^Ce(PO^)^  is  stable  in  dry  air,  argon  and  hydrogen  atmospheres.  Its  thermal 
expansion  coefficients  in  dry  hydrogen  and  air  atmosphere  are  18x10'^  /°C  and  15x10'^ 
/°C,  respectively.  Raman  and  NMR  characterization  results  indicate  that  hydrogen 
phosphate  exists  in  the  Ba^Ce{PO^)^  under  wet  atmospheres  and  dry  hydrogen.  A defect 
structure  was  proposed  based  on  characterization  results.  The  relationships  between 
concentration  of  defects  and  oxygen  partial  pressure  or  water  partial  pressure  are 
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predicted.  The  conductivities  of  Ba^Ce{PO^)^  will  be  measured  in  the  next  chapter  to 
assure  the  proposed  defect  structure  model. 
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degBe(29) 


Figure  4-1.  XRD  patterns  of  calcined  precipitate  with  different  calcination  temperatures. 

(a)  XRD  pattern  between  20=10-80°,  (b)  magnified  XRD  pattern  between 
20=25-45°  and  (c)  XRD  peak  intensity  versus  calcination  temperature. 


intensity  (a.u.)  intensity  (a.u.) 
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JCPDS  43-0643 


Figure  4-2.  Compare  XRD  results  from  JCPDS,  experiment  and  calculation  by  CaRIne 
Crystallography  V3.1  software. 


Figure  4-3.  IR  of  precipitates  after  calcining  at  different  temperatures  for  4 h. 
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Figure  4-4.  IR  spectra  of  Ba^Ce(PO^)^  and  composite  BCGSPd. 


Figure  4-5.  Bright  field  image  and  selected  area  diffraction  pattern  of  Ba^Ce{PO^  )^ 
particle,  (a)  Bright  field  image  and  (b)  selected  area  diffraction  pattern. 
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Figure  4-6.  Crystal  structure  of  Ba^Ce(PO^)^  using  “MS  modeling”  software,  (a)  3-D 
view,  (b)  [001]  directional  view,  and  (c)  [111]  directional  view  of  the  crystal 
structure. 
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Figure  4-7.  Crystal  structure  of  Ba^Ce(PO^)^  using  “CaRIne  Crystallography  V3.1” 
software,  (a)  3-D  view  and  (b)  [001]  directional  view  of  Baj,La{PO^)^  (from 
ref  [76]). 
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Figure  4-8.  Directional  view  of  Ba^Ce(P0^  )j . (a)  [010]  and  (b)  [1 1 1], 
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Raman  shift  (cm'^) 


Figure  4-9.  Raman  spectra  of  Ba^Ce(PO^)^  after  dry  air,  N2  and  H2  treatment,  (a)  at  400- 
3600  cm"'  and  (b)  at  900-1100  cm"'  region. 


Raman  intensity  (a.u.) 
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Raman  shift  (cm  ') 


Figure  4-10.  Raman  spectra  of  Ba^Ce(PO^)^  after  wet  and  dry  nitrogen  treatment,  (a,  b) 

at  400-3600  cm''  and  (c)  at  900-1 100  cm'*  region.  The  band  intensities  in  (a) 
were  adjusted  to  show  the  band  position  in  (b). 


Figure  4-11.  Raman  spectra  of  Ba^Ce(PO^)^  after  wet  hydrogen  treatment. 
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Figure  4-12.  High-  resolution  NMR  spectra  of  Ba^Ce(PO^)^  after  different  gas 
treatment  at  550°C  for  24  h. 


Figure  4-13.  High-temperature  XRD  pattern  of  Ba^Ce(PO^)^ . 
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Figure  4-14.  XRD  of  Ba^Ce{PO^)^  after  different  atmosphere  treatments  at  600  °C  for 
10  h. 


Figure  4-15.  Thermal  expansion  behavior  of  Ba^Ce(PO^)^  in  hydrogen  and  air 
atmosphere. 
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Table  4-1 . Measured  spacing  d compared  with  those  obtained  from  TEM  result  and 
literature  value*. 


Coprecipitation  method 
a=l  0.4466  A 

Literature 
a=10.5112  A 

Coprecipitation  method 
a=l  0.4466  A 

Literature 
a=10.5112  A 

(hkl) 

dexperiment. 

dlEM. 

dliterature 

(hkl) 

dexperiment. 

dlEM. 

dliterature 

211 

4.2458 

4.2648 

4.2922 

631 

1.5449 

1.5403 

1.5496 

220 

3.6897 

3.6934 

3.7145 

444 

1.5128 

1.5078 

1.5174 

310 

3.2960 

3.3035 

3.3244 

721 

1.4275 

1.4216 

1.4305 

321 

2.7916 

2.7919 

2.8093 

642 

1.4017 

1.3960 

1.4046 

332 

2.2271 

2.2272 

2.2415 

730 

1.3768 

1.3717 

1.3804 

422 

2.1361 

2.1324 

2.1459 

651 

1.3322 

1.3267 

1.3349 

510 

2.0529 

2.0487 

2.0615 

800 

1.3120 

1.3058 

1.3138 

521 

1.9128 

1.9073 

653 

1.2542 

1.2486 

1.2562 

611 

1.7016 

1.6947 

1.7052 

660 

1.2379 

1.2311 

1.2388 

620 

1.6561 

1.6518 

1.6624 

830 

1.2209 

1.2144 

1.2219 

*In  coprecipitation  method,  coprecipitates  were  sintering  at  1100  °C  for  4 h; 

In  literature  JCPDS  43-0643,  mixture  were  sintered  at  1000  °C  for  one  day  and  1350  °C 
for  10  days 


Table  4-2.  Wavenumbers  (cm'')  of  the  absorption  bands  in  infrared  spectra  of 
Ba^Ce(PO^)^ 


ui(cm'') 

U2(cm'') 

U3(cm'') 

U4(cm'') 

Troom 

3465 

1635 

1093,995 

612,558 

T600°C 

3391 

1653 

1093,997 

610,558 

H 

00 

o 

o 

o 

O 

3391 

1097,  996 

618,558 

Tiioo‘>c 

1050 

566 

ui;  stretching  vibration  of  water 
U2:  bending  vibration  of  water 
U3:  asymmetric  stretching  vibration, 

U4:  asymmetric  bending  vibration,  S^^P02 


Table  4-3.  Raman  shift  of  Ba^CeiPO^)^  at  dry  atmosphere  compared  with  Sr  doped 
LaP04 


Functional 

Group 

Raman  shift  (cm'*) 
Sr-doped  LaP04  [42] 

Raman  shift  (cm'*) 
Ba^CeiPO^)^ 

POl 

1073,  1065, 1055  and 
1025 

asymmetric 

stretching 

1057,  1020 

991  and  967 

symmetric 

stretching 

944 

619,  589,571,559  and  537 

O-P-0  F2-bending 

571 

466,416  and  395 

O-P-0  E-bending 

463  and  439 

PiO^i 

1036,764,  728 

1057  and  1020  [41] 
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Table  4-4.  Raman  shift  of  Ba^Ce{PO^)^  at  wet  atmosphere  compared  with  Sr  doped 
LaP04 


Raman  (cm'‘) 

Raman  (cm'*)  Ba3Ce(P04)3 

Sr-doped  LaP04  [42] 

wet 

wet 

assignment 

N2 

H2 

POl~ 

1073 

asymmetric  stretching 

1065 

1055 

1059 

1025 

1020 

991 

symmetric  stretching 

967 

944 

944 

619 

asymmetric  bending 

589 

571 

571 

577 

559 

537 

466 

symmetric  bending 

416 

438 

438 

395 

HPOl~ 

stretching  of  water 

2881 

stretching  of  POH 

2848 

1439 

1350 

bending  of  POH 

1295 

1135 

1131 

POH  in  plane  bending 

672 

POH  out-of-plane  bending 

Table  4-5.  Chemical  shift  of  *H  NMR  of  Ba^Ce{PO^)^  after  different  atmosphere 


treatments  compared  with  Sr  doped  LaP04 


Sr-doped  LaP04  [42] 

Ba3Ce(P04)3 

dry  air 

dry  H2 

dry  N2 

wet  N2 

Water  molecules 
5.4  and  1.2  ppm 

5 ppm 

largest  intensity 

5 ppm 

5 ppm 

5 ppm 

Proton 
12.4  ppm 

15  ppm 

largest  intensity 

15  ppm 
(noise) 

15  ppm 
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Table  4-6.  Defect  concentration  dependence  of  , Po^  and  p„^o  ■ 


Po, 

Ph^o  5 Ph^ 

n 

[(^2^7)2^04  ] 

[(HPO,)po,] 

Ph^  > P02  ’ Ph^o 

Ph2  ’ P02  y PH2O 

Ph2  y P02  y PH2O 

Low 

Low  p^J^o 

1 

1 

no 

No  free  H2 

P02 

Pol 

Low  ph^o 
With  p„^ 

I 1 1 

Ph^  P H2O  ’ Pol 

1 1 1 
P^Puloy  Pol 

negligible 

1 1 1 
Ph2  ’ Pol  Ph^o 

negligible 

1 1 1 

Ph2^  P0IPH2O 

High  ph^o 
No  free  H2 

I 1 

Pc^^P^^O 

I 1 1 

Ph^  P h^o  ’ P^ 

negligible 

1 1 1 
Ph2  P H2O  y Pol 

1 1 

Pol  PH2O 
negligible 

High  pti^o 
With  p„^ 

1 I 1 

Ph2  ’ Pol  PhjO 

negligible 

1 1 1 
Ph2  ’ Pol  PH2O 

Low- 

intermediate 

Low/7^^0 
No  free  H2 

4 

4 

no 

High  ph^o 
No  free  H2 

1 1 

PoIPh^o 

1 I 1 

Ph^PhIo’’  Pol 

negligible 
1 1 1 

Ph2  P 4o  y P(4 

1 1 

P(4  P^o 

negligible 

High  Pfj^o 
With  p„^ 

1 I 1 

phi  ’ Pol  Ph^o 
1 1 1 

Ph^  P h\o  ’ Pol 

negligible 

I 1 1 

Ph2  P 4o  y Pol 

1 1 1 
Ph2  ’ Pol  PH2O 
negligible 

Intermediate 

Low  ph^o 

n,  p,  [{PiOiyiPoJ  are  low, 
no  apparent  dependence 

no 
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n,  P,  [iP202)2P0,] , [{HPO^^ypoJ  are  low,  no  apparent 

dependence 

High 

Electroneutra 

ity  condition  faces  a challenge,  no  further  discussion 

CHAPTER  5 

CONDUCTIVITY  MEASUREMENT  OF  Ba,Ce{PO,), 

5.1  Introduction 

The  defect  structure  of  Ba^CeiPO^)^  was  proposed  in  chapter  4 based  on 
characterization  results.  The  relationship  between  charge  carrier  concentration,  oxygen 
partial  pressure  and  water  partial  pressure  was  derived.  The  measured  conductivities  of 
Ba^CeiPO^)^  are  presented  in  this  chapter  to  test  the  proposed  defect  structure  model. 

5.2  Experiment 

Ba^Ce(PO^)^  powder  was  pressed  into  pellets,  cold  isostatic  pressed  at  200  MPa  and 
sintered  at  1100,  1200,  1300  ”C  for  8 h.  SEM  was  used  to  characterize  the  microstructure 
of  the  sintered  pellets.  Platinum  paint  was  applied  on  each  side  of  polished  Ba^CeiPO^)^ 
pellets  and  platinum  wire  was  fastened  to  the  electrode.  They  were  sintered  at  1050"C  for 
half  an  hour. 

Conductivities  were  measured  by  two-point  method  with  a Solartron  1260 
frequency  response  analyzer.  For  conductivity  measurements  in  wet  and  dry  oxygen,  air, 
nitrogen  and  4vol%  hydrogen  atmospheres  in  the  300-600  °C  region,  gas  flow  rates  were 
controlled  at  50  seem. 

For  conductivity  measurements  as  a function  of  oxygen  partial  pressure,  the  oxygen 
partial  pressure  range  10'^^~10''^  atm  was  achieved  by  mixing  argon  with  hydrogen.  The 
oxygen  partial  pressure  range  10“  ~10'  atm  was  achieved  by  purifmg  argon  gas  through 
titanium  charger  at  800°C  to  remove  oxygen.  The  oxygen  partial  pressure  range  10''‘’~10' 
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atm  was  obtained  by  bubbling  purified  argon  gas  through  water  at  room  temperature. 
The  oxygen  partial  pressure  range  10'^^~10'^  atm  was  also  obtained  by  bubbling  mixture 
of  hydrogen  and  argon  through  water  at  room  temperature.  The  oxygen  partial  pressure 
range  10"  ~10'  atm  was  achieved  by  using  an  oxygen  pump  (SGM5EL)  with  argon  gas. 
The  oxygen  partial  pressure  range  0.2 1~1  atm  was  obtained  by  mixing  argon  and  oxygen. 

An  oxygen  sensor  (SGM5EL)  with  operating  temperature  750  °C  monitored  the 
oxygen  partial  pressure.  The  actual  oxygen  partial  pressure  at  the  testing  temperature  was 
calculated  using  thermodynamic  equations. 

5.3  Results  and  Discussion 
5.3.1  Sintering  and  Density  Issue 

Figure  5-1  shows  the  relationship  between  the  sample’s  density  and  sintering 
temperature  with  a dwell  time  of  8 h.  Figure  5-2  shows  XRD  results  with  different 
sintering  temperatures  and  dwell  times.  Figure  5-3  shows  the  SEM  cross-section  view  of 
the  sintered  pellets  with  magnification  3500.  As  expected,  the  relative  density  increased 
with  sintering  temperature,  from  78%  at  1100  °C,  to  89%  at  1200  °C,  to  95%  at  1300  °C. 
However,  the  pellet  sintered  at  1300  °C,  changed  its  color  to  green  and  a glassy  phase 
was  observed.  The  XRD  result  (shown  in  figure  5-2)  shows  that  impurity  phase  Ba2P207 
was  produced  at  high  temperature.  For  pellets  sintered  at  1100  and  1200°C  for  8h,  there 
was  no  second  phase  detected.  This  is  similar  to  Sr  doped  LaP04.  Sr  doped  LaP04 
sintered  at  1350  ”C  has  a second  phase,  Sr2P207,  and  the  resulting  conductivity  of  the 
material  is  lower  than  that  sintered  at  1200  °C  [42].  Therefore,  we  chose  1200  °C  as  the 
sintering  temperature.  We  sintered  pellets  at  1200  °C  for  8 h,  16  h and  24  h,  respectively. 
The  relative  density  was  91%  with  a sintering  temperature  of  1200  °C  for  24  h.  However, 
XRD  in  figure  5-2  shows  the  existence  of  impurity  phase  Ba2P207,  which  would  lower 
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conductivity.  A dwell  time  of  8 h was  chosen  to  assure  no  second  phase  was  produced 
during  the  sintering.  The  porosity  decreased  extensively  with  increasing  sintering 
temperature  based  on  a SEM  cross-section  image.  There  were  no  interconnected  pores  in 
the  cross  section  with  a relative  density  89%. 

Samples  used  for  conductivity  measurements  were  batch  sintered  at  1200  °C  for  8 
h.  Figure  5-3  shows  the  cross-section  image  of  sintered  Ba^Ce(PO^)^  with  magnification 
3500.  Figure  5-4  shows  the  SEM  surface  image  of  a porous  Pt  electrode  on  a 
Ba^Ce(PO^)^  sample  with  magnification  1500. 

5.3.2  Conductivity 

5.3.2. 1 Conductivity  in  different  atmospheres 

Atmosphere  and  Water  Vapor  Effects  on  the  Impedance  Plots 

Figure  5-5  (a-d)  illustrates  the  atmosphere  and  water  vapor  effects  on  the 
impedance  plot  of  Ba^Ce{PO^)^  at  600°C.  The  impedance  discontinuity  at  frequency 
around  1x10^  FIz  is  an  artifact.  Figure  5-5  (e-f)  enlarges  the  impedance  plot  of 
Ba^Ce(PO^)^  in  wet  and  dry  4vol%  hydrogen  in  figure  5-5  (d)  to  show  the  electrolyte 

and  electrode  semicircle.  The  semicircle  with  high  frequency  belongs  to  the  electrolyte 
and  the  semicircle  with  low  frequency  belongs  to  the  electrode.  The  capacitances  of  the 
high  frequency  and  low  frequency  semicircles  are  on  the  order  of  pF  and  nF,  respectively. 

Under  pure  oxygen,  the  impedance  of  Ba^Ce(PO^)^  changed  little  by  the  presence 
of  water  vapor.  Under  pure  air  or  nitrogen,  the  impedance  decreased  slightly  in  the 
presence  of  water  vapor,  suggesting  the  existence  of  proton  conduction.  Under  4vol% 
hydrogen,  the  impedance  decreased  significantly  in  the  presence  of  water  vapor. 
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Impedances  of  Ba^Ce{PO^)^  also  decreased  significantly  under  4vol%  dry  hydrogen 

compared  with  those  under  dry  nitrogen  or  oxygen. 

The  proposed  defect  structures  in  chapter  4 are  used  to  explain  these  atmosphere 
and  water  vapor  effects  on  the  impedance  plot  of  Ba^Ce(PO^)^ . 

Oxvsen  Atmosphere  (Hish  Oxvsen  Partial  Pressure  Resion): 

The  defect  equilibria  under  oxygen  atmosphere  are 

2iPO,Vpo,  ^ iP20^)2Po,  +2e  +^02 . K1  (4-1) 

^iPO.Ypo,  +H20^2{HPO,)po,  +2e  +^0^ . K6  (4-6) 

The  electroneutrality  condition  can  be  written  as 
[(NPO^)  pQ^  ] + ^[(P2^7)2POt  ] + P = • 

Electrons  are  the  only  negative  defects  to  charge  balance  holes,  oxygen  vacancies 
and  possible  protons.  The  oxygen  vacancy  concentration  could  be  omitted  in  this  region. 
However,  electron  concentration  is  usually  not  higher  than  hole  concentration  under  high 
oxygen  partial  pressure  region.  Therefore,  the  proton  concentration  will  not  increase 
much  in  the  presence  of  water  vapor  based  on  the  electroneutrality  restriction  and  as 
indicated  by  the  impedance  data  of  Ba^CeiPO^)^  (fig-  5-5a). 

Nitrosen  Atmosphere  (Intermediate  Oxvsen  Partial  Pressure  Resion): 

Dry  nitrogen  belongs  to  the  intermediate  oxygen  partial  pressure  region.  The 
Raman  spectroscopy  result  in  figure  4-9  suggests  that  oxygen  vacancies  exist  in  the  dry 
nitrogen  treated  sample.  The  corresponding  defect  equilibrium  in  the  dry  atmosphere  can 
be  written  as 

2{PO^) pon  (^2*2? )2/>04  ■'■■^^2  +2e 


K1 


(4-1) 
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The  electroneutrality  condition  can  be  written  as 

Wet  nitrogen  also  belongs  to  the  intermediate  oxygen  partial  pressure  region. 
Raman  spectroscopy  result  in  figure  4-9  confirmed  the  proton  incorporation  by  the 
stretching  modes  at  2881,  2848cm’^  and  bending  modes  at  1439,  1295,  1 131cm''  of  POH 
bond.  The  corresponding  defect  equilibrium  in  wet  atmosphere  can  be  written  as 

2(PO,ypo,  +H^O-^2{HPO,)po,  +2e'  . K6  (4-6) 

The  electroneutrality  condition  can  be  written  as 
[{HPO^ ) pQ^  ] + 2[{P2  0-j  ) 2PO,  ] + /^  = ” • 

Protons  are  incorporated  into  the  lattice  and  make  a contribution  to  the  total 
conductivity  in  the  presence  of  water  vapor.  However,  the  electron  concentration  almost 
equals  the  hole  concentration  in  the  intermediate  oxygen  partial  pressure  region,  which 
restricts  the  proton  concentration  based  on  the  electroneutrality  condition.  Therefore,  the 
impedance  will  increase  but  will  not  increase  much  in  the  presence  of  water  vapor  in  the 
nitrogen  atmosphere,  which  is  demonstrated  by  figure  5-5  (c). 

Hvdrosen  Atmosphere  (Low  Oxvsen  Partial  Pressure  Resion): 

For  dry  hydrogen  treated  Ba^CeiPO^)^,  chemical  shift  at  15  ppm  in  the  'H  NMR 

spectroscopy  result  in  figure  4-12  clearly  demonstrates  proton  insertion  into  the  lattice. 
The  corresponding  defect  equilibrium  can  be  written  as 

(PO,  Ypo^  +^H2^  {HPO, )po^  + e K3  (4-3) 

For  wet  hydrogen  treated  Ba^CeiPOy)^ » the  broad  band  at  1350  cm''  in  Raman 
spectroscopy  result  in  figure  4-11  is  assigned  to  the  bending  mode  of  POH  bonds. 
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suggesting  the  proton  incorporation.  The  corresponding  defect  equilibrium  can  be  written 
as 

iPO.Vpo,  (HPO, )po,  + e'  K3  (4-3) 

2{PO^)po^+H20^2{HP04)po^+2e' +^02-  K6  (4-6) 

The  electroneutrality  condition  in  wet  or  dry  hydrogen  can  be  written  as 
[{HPO,)po^]  + 2[{P2 

^7)2^04]-'*  • 

With  the  presence  of  hydrogen  gas,  proton  and  electron  concentration  increase  due 
to  the  reaction  (4-3).  The  mobility  of  proton  or  electron  is  much  larger  than  that  of 
oxygen  vacancy.  Therefore,  the  impedance  of  Ba^CeiPO^)^  decrease  significantly  under 

dry  hydrogen  atmosphere  compared  with  dry  oxygen  or  dry  nitrogen. 

Free  hydrogen  is  the  only  proton  source  in  dry  hydrogen  while  both  water  vapor 
and  free  hydrogen  provide  protons  in  wet  hydrogen.  In  the  Raman  results  of  figure  4-11, 
the  POH  bending  band  intensity  is  very  strong  in  wet  hydrogen  treated  Ba^CeiPO^)^ 

while  it  is  almost  undetectable  in  dry  hydrogen  treated  sample.  The  intensity  difference 
indicated  that  water  vapor  plays  a much  more  important  role  and  proton  concentration  is 
much  higher  in  wet  hydrogen.  Therefore,  the  conductivity  contribution  from  protons  is 
expected  to  be  larger  in  wet  hydrogen  than  in  dry  hydrogen,  which  is  demonstrated  by 
figure  5-5  (d). 

Atmosphere  and  Water  Vapor  Effects  on  the  Arrhenius  Plots 

The  bulk  electrical  resistance  (R)  values  were  obtained  from  the  low  frequency 
intercept  of  the  first  semicircle  on  the  real  axis  of  the  impedance  plot.  The  conductivity 
( (T  ) can  be  expressed  as 
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L 

a = 

A*R 


L 

-nd^R 

4 


(5-1) 


where  A is  the  cross  sectional  area. 

The  electrical  transport  behavior  can  be  represented  by  the  equation; 

aT  = CTo  exp(-  -^)  (5-2) 

kT 

where  o-q  is  the  pre-exponential  term,  Ea  is  the  activation  energy,  k is  Boltzmann 
constant  and  T is  absolute  temperature.  The  plot  of  log(cr7’ ) versus  1/T  yields 
information  about  the  pre-exponential  term  and  the  activation  energy. 

Figure  5-6  shows  a Arrhenius  plot  of  log(o-r)  versus  1000/T  for  Ba^Ce(PO^)^  in 

dry  and  wet  oxygen,  air,  nitrogen,  and  4vol%  hydrogen  at  300-600  °C.  Table  5-1  lists  the 
electrical  conductivity  data  of  Ba^Ce(PO^  )^  in  different  atmospheres.  The  conductivity 
of  BajCe(PO^)^  at  600°C  was  8.93x10'^  S/cm  in  wet  air  and  5.8x10'^  S/cm  in  dry  air. 

This  conductivity  behavior  indicates  that  proton  conduction  exists.  The  activation  energy 
in  dry  and  wet  air  was  0.821  and  0.810  eV,  respectively,  which  is  similar  to  the  activation 
energy  reported  for  Ca  or  Sr  doped  LaP04  at  moderate  temperature  [27].  The 
conductivity  of  BajCe(PO^)^  at  600  °C  was  8.91x10'^  S/cm  in  wet  oxygen  and  8.78x10' 

^ S/cm  in  dry  oxygen.  The  total  conductivity  does  not  change  much  in  the  wet  oxygen 
atmosphere.  The  activation  energy  in  the  dry  and  wet  oxygen  was  0.781  and  0.776  eV, 
respectively.  The  conductivities  of  Ba^Ce(PO^)^  at  600  °C  in  the  dry  and  wet  nitrogen 

were  7.05x10'^  and  7.67x10'^  S/cm,  respectively.  The  activation  energies  were  0.80  and 
0.829  eV  in  the  dry  and  wet  nitrogen,  respectively.  The  conductivities  of  Ba^Ce(PO^)^ 
at  600  °C  in  dry  and  wet  4vol%  hydrogen  were  4.04x10'^  and  1.46x10''*  S/cm, 
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respectively.  The  activation  energies  were  0.932  and  0.815  eV  in  the  dry  and  wet 
hydrogen,  respectively.  This  is  the  first  time  that  proton  conductivities  were  studied  in  the 
low  oxygen  partial  pressure  region,  especially  in  the  hydrogen  atmosphere. 

Figure  5-7  compares  the  conductivity  of  Ba^Ce(PO^)^  with  other  phosphates  in 

the  wet  air  and  oxygen  atmospheres.  The  conductivity  of  Ba^Ce(PO^)^  was  8.91x10'^ 

S/cm  at  600  °C  in  wet  oxygen.  The  conductivities  of  un-doped  LaP04  and  YPO4  were 
8x10'^  S/cm  and  1.26x10'^  S/cm  at  600  °C  in  wet  oxygen,  respectively.  The  conductivity 
of  5 mol%  Sr  doped  LaP04  was  1.5x10’^  S/cm  at  600  °C  in  wet  oxygen.  Therefore,  the 
conductivity  of  Ba^Ce(PO^)^  is  higher  than  un-doped  LaP04  and  YPO4  but  lower  than  5 
mol%  Sr  doped  LaP04. 

For  Sr  doped  LaP04,  both  electrons  and  dopant  can  charge  balance  oxygen 
vacancies,  holes  and  possible  protons.  The  electroneutrality  condition  can  be  written  as 

[{HPO,)po,  ] + ] + P = [Sr^]  + n. 

For  un-doped  phosphates,  the  electrons  are  the  only  negative  defects  to  charge 
balance  oxygen  vacancies,  holes  and  possible  proton.  The  electroneutrality  condition  can 
be  written  as 

[(HPO^) PQ^  ]+  21(^26*7 )2POt  ] + p = n . 

The  proton  concentration  increases  in  doped  or  un-doped  phosphates  in  the 
presence  of  free  hydrogen  or  water  vapor;  however,  the  positive  defect  concentration 
needs  to  equal  the  sum  of  the  electron  concentration  and  the  dopant  concentration  in  the 
doped  phosphate  while  it  is  equal  to  the  electron  concentration  alone  in  the  un-doped 
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phosphate.  The  result  is  a total  conductivity,  which  is  much  higher  in  doped  phosphates 
than  un-doped  ones  due  to  dopant  effects. 

In  the  future,  Ba^Ce(PO^)^  will  be  doped  at  Ce  site  with  divalent  elements  and  its 
conductivity  behavior  will  be  studied.  Higher  conductivity  is  expected  in  the  doped 
Ba^Ce(PO^)^. 

S.3.2.2  Conductivity  as  a function  of  oxygen  partial  pressure 

The  defect  structure  model  in  chapter  4 derived  the  relationship  between  the  charge 

carrier  concentration  of  Ba^Ce(PO^)^  and  oxygen  partial  pressure.  To  test  the  proposed 
defect  structure,  the  conductivity  of  Ba^Ce(PO^)^  was  measured  as  a function  of  oxygen 
partial  pressure. 

It  should  be  pointed  out  that  the  conductivity  is  not  only  related  to  charge  carrier 
concentration  but  also  related  to  carrier  mobility;  thus,  the  conductivity  dependence  of 
oxygen  partial  pressure  can  only  approximately  reflect  the  charge  carrier  concentration 
dependence  of  oxygen  partial  pressure. 

Figure  5-8  shows  the  conductivity  of  Ba^Ce(PO^)^  as  a function  of  oxygen  partial 

pressure  at  600  °C  and  550  °C.  The  figure  clearly  demonstrates  the  total  conductivity 
dependence  in  the  whole  oxygen  partial  pressure  region.  Oxygen  partial  pressure  range  is 
divided  into  four  regions:  high  (l~0.21atm,  IV),  intermediate  (0.21-10'^ atm.  III),  low- 
intermediate  (10'^~10'^^atm,  II)  and  low  (10'’^~10'^^atm,  I).  The  conductivity 
dependences  of  oxygen  partial  pressure  in  these  regions  of  figure  5-8  are  shown  in  figure 
5-9  to  5-15.  They  were  compared  with  defect  concentration  dependences  of  oxygen 
partial  pressure,  which  was  derived  from  the  defect  structure  model.  The  results  are 


shown  in  table  5-2. 
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Low  Oxygen  Partial  Pressure  Region  (P02=  atm) 

Four  cases  were  discussed  in  the  proposed  defect  structure  model.  Case  1 is  the  low 
water  partial  pressure  condition  without  free  hydrogen,  case  2 is  the  low  water  partial 
pressure  condition  with  free  hydrogen,  case  3 is  the  high  water  partial  pressure  condition 
without  free  hydrogen,  and  case  4 is  the  high  water  partial  pressure  condition  with  free 
hydrogen. 

Case  3,  however,  could  not  be  studied  experimentally  in  this  work  due  to 
equipment  limitation.  The  proposed  defect  structure  model  is  only  tested  in  cases  1 , 2 and 
4. 

Case  1:  low  water  partial  pressure  condition  without  free  hydrosen 

When  an  oxygen  partial  pressure  is  obtained  by  purifying  dry  argon  gas,  it 
corresponds  to  the  low  water  vapor  partial  pressure  condition  without  the  existence  of 
free  hydrogen.  The  experimentally  obtained  oxygen  partial  pressure  range  is  10'^°~10''^ 
atm. 

According  to  the  proposed  defect  structure  model,  oxygen  vacancies  are  the  only 
significant  positive  defects  to  charge  balance  electrons.  The  oxygen  vacancy  formation 
can  be  expresses  as  the  following  defect  equilibrium: 

^{PO.Ypo,  ^iPiOYiPO,  +2^'  +\o^  ■ K1  (4-1) 

The  electroneutrality  condition  is  simplified  as 

21(^2^7)2^04  ] ~ • 

\_ 

K,=  {nl2]n^pl^ 
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Electron  concentration  and  oxygen  vacancy  concentration  are  derived  to  be 
1 

proportional  to  . The  experimentally  obtained  total  conductivity  is  proportional  to 
--  -i  i 

Pol  600°C,  at  550°C  and  p^^  at  500“C,  which  is  shown  in  figure  5-9.  Therefore, 

the  conductivity  dependency  of  oxygen  partial  pressure  supports  the  proposed  defect 
structure  under  low  water  partial  pressure  condition  without  free  hydrogen. 

Case  2:  low  water  partial  pressure  condition  with  free  hydrogen 

When  an  oxygen  partial  pressure  is  obtained  by  mixing  dry  hydrogen  and  dry 
argon,  it  corresponds  to  the  low  water  vapor  partial  pressure  condition  with  free 
hydrogen.  The  experimentally  obtained  oxygen  partial  pressure  range  is  10'^^~10’'^  atm. 
Defect  equilibrium  under  this  condition  can  be  written  as 

(PO4  {HPO,)po,  + e K3  (4-3) 

There  exists  small  amount  of  water  vapor  due  to  the  existence  of  hydrogen  and 
oxygen  formed  in  formula  (4-1): 

K4  (4-4) 

^4  = Ph^oPh-^Po]  • 

Oxygen  vacancies  and  inserted  protons  are  the  positive  species  in  the  lattice.  The 
electroneutrality  condition  is  written  as 

[{HPO,)po,]  + 2[{P^O,)po,]  = n . 

If  oxygen  vacancies  dominate,  then. 


2[(^2f^7  )/>04  ] ~ ^ 
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lK,p^l=n\ 

-i 

Electron  and  oxygen  concentration  are  derived  to  be  proportional  to  . 

If  proton  concentration  dominates,  then, 

[{HPO,)po,]  = n 

_i  1 -1 

wpIopo:  =«' 

I 

Electron  and  proton  concentration  are  derived  to  be  proportional  to  p^^  p^  ^ . Water 
partial  pressure  is  low  and  can  be  considered  as  constant  in  case  2;  therefore,  the  electron 

and  proton  concentration  are  most  likely  to  be  proportional  to  p^^ . 

As  we  mentioned,  case  2 was  obtained  experimentally  by  mixture  of  dry  argon  and 
dry  hydrogen.  Lower  oxygen  partial  pressure  is  obtained  by  increasing  hydrogen 
concentration  in  the  gas  mixture,  which  provides  more  proton  sources  at  the  same  time. 
Therefore,  more  protons  are  incorporated  into  the  lattice  and  proton  concentration 
gradually  becomes  dominant  when  oxygen  partial  pressure  becomes  lower  in  case  2. 
According  to  the  proposed  defect  structure,  electron  and  proton  concentration  are 

proportional  to  p^^  when  proton  concentration  dominates,  which  is  smaller  than  p^^ 
when  oxygen  vacancy  concentration  dominates. 

The  experimentally  obtained  conductivity  is  initially  proportional  to  p^^  and 
conductivity  dependence  decreases  when  oxygen  partial  pressure  becomes  lower,  which 
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is  shown  in  figure  5-10  and  may  indicate  the  transition  from  oxygen  vacancy  dominance 
to  proton  dominance. 

Therefore,  the  conductivity  behavior  of  Ba^Ce(PO^)^  is  considered  as  a support  to 

the  proposed  defect  structure  under  low  water  partial  pressure  condition  with  free 
hydrogen. 

Case  4:  hieh  water  partial  pressure  condition  with  free  hvdro2en 

When  oxygen  partial  pressure  range  10''^~10'^^  atm  is  obtained  by  flowing  a 
mixture  of  hydrogen  and  argon  through  water,  it  corresponds  to  the  high  water  partial 
pressure  condition  with  free  hydrogen. 

The  corresponding  defect  equilibriums  can  be  written  as 

(^^4  + y ^2  ^ +e  K3  (4-3) 

2{HPO, )po,  + 2e'  + 1 • K6  (4-6) 

The  electroneutrality  condition  can  be  written  as 
[{HPO^ ) pQ^  ] + 2[(P2  Oj  ) 2po^  ] = w • 

With  free  hydrogen  and  water  vapor  contribution,  proton  concentration  is  expected 
to  be  dominant.  The  electroneutrality  condition  is  simplified  as 

KHPO,)po^]  = n. 

_i  i _i  i 

_l_  I 

Proton  and  electron  concentration  are  derived  to  be  proportional  to  P(Y q.  They 
are  proportional  to  if  water  vapor  is  kept  constant.  The  obtained  conductivity 
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dependence  of  oxygen  partial  pressure  is  (shown  in  figure  5-11),  which  is  almost 

negligible  when  oxygen  partial  pressure  decreased  below  10'^^  atm. 

Oxygen  partial  pressure  change  is  very  small  in  the  low  oxygen  partial  pressure 

region  and  water  partial  pressure  may  play  a more  important  role.  The  water  partial 

1 

pressure  was  held  constant  in  the  experiment;  therefore,  ^ did  not  change  much.  It 

is  difficult  to  tell  whether  or  not  the  conductivity  behavior  in  case  4 supports  the  defect 
structure  model  in  the  corresponding  region. 

Low-Intermediate  Oxygen  Partial  Pressure  Region  (P02=  10'^~10'*^  atm) 

In  the  proposed  defect  structure  model,  cases  1,  3 and  4 were  discussed  in  the  low- 
intermediate  oxygen  partial  pressure  region.  Due  to  the  equipment  limitation,  only  cases 
3 and  4 were  tested  in  this  work.  Case  3 is  the  high  water  partial  pressure  condition 
without  free  hydrogen  and  case  4 is  the  high  partial  pressure  condition  with  free 
hydrogen. 

Case  3:  hish  water  partial  pressure  condition  without  free  hydrosen 

When  an  oxygen  partial  pressure  is  obtained  by  flowing  purified  argon  through 
water,  it  corresponds  to  the  high  water  vapor  partial  pressure  condition  without  free 
hydrogen.  The  experimentally  obtained  oxygen  partial  pressure  range  is  10‘''’~10‘"  atm. 
Defect  equilibrium  is  written  as 

^(PO^Vpo,  +H^O^  2{HPO^ )po^  +2e'+j02.  K6  (4-6) 

The  electroneutrality  condition  is  written  as 


[(^PO^)pQ^  ] + 2[(P2  D7  )2/>o^  ] — ” • 
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If  oxygen  vacancy  concentration  dominates,  electroneutrality  condition  is 
simplified  as 

2[(^2  ^7  )lPOi  ] ~ ^*  • 

Electron  concentration  and  oxygen  vacancy  concentration  are  derived  to  be 
proportional  to  . 

If  proton  concentration  dominates,  the  electroneutrality  condition  is  written  as 
_i  i 

WpkoP'i 

I 1 

Proton  and  electron  concentration  are  derived  to  be  proportional  to  Po"Ph,o- 

Only  three  conductivity  data  points  were  obtained  due  to  the  equipment  limitations. 
The  conductivity  dependency  of  oxygen  partial  pressure  was  obtained  by  extending  the 
data  trend  to  oxygen  partial  pressure  = 10’^  atm.  The  conductivity  data  trend  shows  a 

Pq'^  dependence  at  600“C  and  a p^^^  dependence  at  550  and  500“C,  which  is  shown  in 

figure  5-12.  We  could  not  tell  whether  or  not  the  conductivity  dependency  of  oxygen 
partial  pressure  supports  the  proposed  defect  structure  under  this  condition  due  to  the 
limited  data. 

Case  4:  hish  water  yartial  pressure  condition  with  free  hydrosen 

When  an  oxygen  partial  pressure  range  at  lO  '^-lO’’  atm  is  obtained  by  flowing  a 
mixture  of  hydrogen  and  argon  through  water,  it  corresponds  to  the  high  water  vapor 
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partial  pressure  condition  with  free  hydrogen.  The  corresponding  defect  equilibriums  are 
written  as 

iPO,  Vpo,  {HPO.Ypo,  + e'  K3  (4-3) 

2{PO^Ypo^+HjO->2{HPO,)po,+2e  +^0^.  K6  (4-6) 

The  electroneutrality  condition  can  be  written  as: 

[{HPO,)po^]  + 2[(P^O,)^PO^]  = n . 

A small  trace  of  free  hydrogen  in  the  gas  mixture  decreases  oxygen  partial  pressure 
a lot.  Therefore,  the  free  hydrogen  concentration  is  not  high  at  the  low-intermediate 
oxygen  partial  pressure  region  and  its  contribution  to  the  proton  incorporation  can  be 
omitted  compared  with  water  vapor  contribution.  All  the  derivations  in  case  3 also 
approximately  apply  in  case  4. 

If  oxygen  vacancy  concentration  dominates,  electron  concentration  and  oxygen 

vacancy  concentration  will  be  proportional  to  . If  proton  concentration  dominates,  the 

-i  i 

proton  and  electron  concentration  will  be  proportional  to  Ph  o- 

The  experimentally  obtained  conductivity  is  proportional  to  , which  is  shown  in 

figure  5-13  and  may  suggest  that  oxygen  vacancy  dominates  in  this  region.  This  also 
suggests  that  the  conductivity  dependence  of  oxygen  partial  pressure  in  case  4 supports 
the  proposed  defect  structure  model  in  the  low-intermediate  oxygen  partial  pressure 


region. 


115 


Intermediate  Oxygen  Partial  Pressure  Region  ('P02  = 0.21-10'^  atml 

Intermediate  oxygen  partial  pressure  range  (10'^~10'‘  atm)  was  obtained  by  using 
an  oxygen  pump  (SGM5EL)  with  dry  or  wet  argon  gas. 

The  defect  equilibriums  can  be  written  as 

2{PO,Ypo^  ^iPiOi)2Po,  +^02  +2e’  K1  (4-1) 

2(PO,  Ypo.+HzO^  2{HPO,  )po^  + 2e’  + i . K6  (4-6) 

The  corresponding  electroneutrality  condition  under  dry  atmosphere  is  written  as 

2[{P lOYipo,  ]+  P - n . 

The  corresponding  electroneutrality  condition  under  wet  atmosphere  is  written  as 

\{HPO^) PQ^  ] + 2[(P2 O7  )2POt  P — ^ • 

Despite  the  fact  that  the  electron,  hole  and  oxygen  vacancy  concentrations  usually 
are  not  high  due  to  the  phosphate  structure  and  the  intermediate  oxygen  partial  pressure; 
they  could  not  be  omitted.  Electron  concentration  almost  equals  hole  concentration, 
restricting  proton  concentration  based  on  the  electroneutrality  condition.  The  oxygen 
vacancy  concentration  is  also  low  in  this  region.  In  the  intermediate  oxygen  partial 
pressure  region,  the  total  defect  concentration  is  expected  to  be  low.  No  apparent  oxygen 
partial  pressure  or  water  partial  pressure  dependency  is  expected  under  this  condition. 

The  experimentally  obtained  conductivity  in  this  region  is  low  and  independent  of 
oxygen  partial  pressure,  which  is  shown  in  figure  5-14.  The  conductivity  dependence  is 
just  as  the  defect  structure  model  predicted.  The  conductivity  behavior  supports  the 
proposed  model  in  the  intermediate  oxygen  partial  pressure  region. 
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High  Oxygen  Partial  Pressure  Region  (P02  = 0.21~1  atm) 

High  oxygen  partial  pressure  range  (0.21~latm)  was  obtained  by  using  a mixture  of 
argon  and  oxygen  gas. 

The  following  defect  equilibriums  exist  in  the  high  oxygen  partial  pressure  region: 

^(PO,Ypo,-^^PiOn)ipo,+2e  K1  (4-1) 

2{PO,Ypo^+H^O--^2{HPO^)po^+2e' K6  (4-6) 

The  electroneutrality  condition  under  dry  atmosphere  can  be  written  as 

2[{P iO-j)2po^  ^ + P = n . 

The  electroneutrality  condition  under  wet  atmosphere  can  be  written  as 

\{HPO^) PQ^  ] + 2[(/’2  0-i)2po^  ] + P = ” ■ 

As  we  pointed  out  in  chapter  4,  the  above  electroneutrality  conditions  face  a 
challenge:  Electrons  are  the  only  negative  defects  to  charge  balance  holes,  oxygen 
vacancies  and  possible  protons;  however,  electron  concentration  is  usually  lower  than 
hole  concentration  in  the  high  oxygen  partial  pressure  region.  Extra  negative  defects  are 
needed  to  charge  balance  holes.  The  defect  concentration  dependencies  of  oxygen  partial 
pressure  were  not  derived. 

Figure  5-15  shows  the  conductivity  of  Ba^CeiPO^)^  as  a function  of  oxygen 

partial  pressure  in  the  oxygen  partial  pressure  region  0.21-1  atm.  The  conductivity 
increases  slightly  with  increasing  oxygen  partial  pressure,  suggesting  p-type  electronic 
conduction.  The  conductivity  of  Ba^CeiPO^)^  is  slightly  higher  in  a wet  atmosphere 
than  in  a dry  atmosphere,  indicating  the  existence  of  proton  conduction.  In  the  presence 
of  water  vapor,  equal  amounts  of  protons  and  electrons  were  incorporated  into  the  lattice. 
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The  electron  concentration  restricts  the  proton  concentration.  The  result  is  that  the  total 
concentration  increases  slightly. 

Based  on  the  conductivity  data  of  Ba^Ce(PO^)^,  it  is  still  unclear  which  extra 
negative  defect  is  charge  balancing  the  holes. 

S.3.2.3  Conductivity  as  a function  of  water  vapor  partial  pressure 

Conductivity  dependence  of  water  vapor  partial  pressure  was  not  studied  in  the 
low-intermediate  oxygen  partial  pressure  region  due  to  equipment  limitation.  It  was  also 
not  studied  in  the  high  oxygen  partial  pressure  region  due  to  the  challenge  we  mentioned 
before. 

Conductivity  dependence  of  water  vapor  partial  pressure  was  only  studied  in  the 
low  and  intermediate  oxygen  partial  pressure  regions.  The  results  are  listed  in  table  (5-3) 
and  compared  with  defect  concentration  dependence. 

Low  Oxygen  Partial  Pressure  Region 

Wet  40vol%  hydrogen  belongs  to  case  4 in  the  low  oxygen  partial  pressure  region. 
Proton  concentration  dominates  in  wet  40vol%  hydrogen.  Proton  concentration  and 

_i  1 

electron  concentration  are  derived  to  be  proportional  to  based  on  the  proposed 

defect  structure  model.  When  water  partial  pressure  is  changed  in  40vol%  hydrogen, 
oxygen  partial  pressure  changes  slightly.  If  oxygen  partial  pressure  change  is  omitted,  the 

proton  concentration  and  electron  concentration  is  proportional  to  pf,  o- 

Figure  5-16  shows  conductivity  of  Ba^Ce(PO^)^  as  a function  of  water  partial 
pressure  in  40vol%  hydrogen  in  the  temperature  range  of  500-600  °C.  The  conductivity 

I j_ 

of  BajCe(PO^)^  is  proportional  to  at  550°C  and  at  600°C.  The 
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experimentally  obtained  conductivity  dependence  of  water  partial  pressure  is  smaller  than 
the  defect  concentration  dependence  in  the  proposed  model. 

It  was  pointed  out  at  the  beginning  of  chapter  5.3. 2.2  that  the  conductivity  is  not 
only  related  to  charge  carrier  concentration  but  also  related  to  carrier  mobility.  Thus,  take 
into  account  of  mobility  effect,  this  conductivity  dependence  of  water  vapor  partial 
pressure  still  approximately  reflects  the  defect  concentration  dependence  of  water  vapor 
partial  pressure  in  the  low  oxygen  partial  pressure  region. 

Intermediate  Oxygen  Partial  Pressure  Region: 

A nitrogen  atmosphere  lies  within  the  intermediate  oxygen  partial  pressure  region. 
Since  electrons  are  present  in  minor  concentration,  the  proton  concentration  is  expected 
to  be  low  and  will  not  increase  much  in  the  presence  of  water  vapor  due  to  the 
electroneutrality  restriction.  The  water  vapor  effect  on  defect  concentration  will  not  be 
apparent. 

Figure  5-17  shows  the  conductivity  of  Ba^Ce{PO^)^  as  function  of  water  partial 
pressure  in  the  temperature  range  of  500-600  °C  for  a nitrogen  atmosphere.  The 
experimentally  obtained  conductivity  is  approximately  proportional  to  at  500  °C, 

which  is  negligible  compared  with  that  in  wet  hydrogen.  The  conductivity  behavior  of 
Ba^Ce(PO^)^  as  a function  of  water  partial  pressure  supports  the  proposed  defect 
structure  in  the  intermediate  oxygen  partial  pressure  region. 

S.3.2.4  Conductivity  as  a function  of  hydrogen  partial  pressure 

Conductivity  dependence  of  hydrogen  partial  pressure  was  studied  only  in  the  low 
oxygen  partial  pressure  regions.  The  result  is  listed  in  table  (5-4)  and  compared  with 
defect  concentration  dependency. 
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Figure  5-18  shows  the  conductivity  of  Ba^Ce{PO^)^  as  a function  of  hydrogen 
partial  pressure  at  550-600  °C  with  =0.03 12  atm.  Wet  hydrogen  belongs  to  case  4 in 
the  low  oxygen  partial  pressure  region.  Proton  concentration  dominates  in  wet  hydrogen. 

Proton  concentration  and  electron  concentration  are  derived  to  be  proportional  to  Ph,- 

The  experimentally  obtained  conductivity  is  proportional  to  at  600°C,  much  smaller 

than  the  derived  pi  from  the  model. 

"2 

One  possible  explanation  for  this  discrepancy  is  the  proton  saturation  in  the  lattice. 
No  more  protons  can  be  incorporated  into  the  lattice  after  the  solubility  limit  is  reached. 
The  other  possibility  is  that  electron-proton  associations  are  formed  at  the  high  proton 
concentration,  without  making  a contribution  to  conduction.  The  association  reaction  can 
be  written  as 

(//PO4  )po,  +e  ^ {{HPO,)po^  e)\ 

Therefore,  the  conductivity  dependence  of  hydrogen  partial  pressure  is  lower  than 
the  model  predicted  one. 

5.3,2,5  Conductivity  under  CO/CO2  mixture  gas 

During  conductivity  dependency  studies,  many  gas  mixtures  were  tried  to  obtain 
the  required  oxygen  partial  pressure.  We  find  that  the  same  oxygen  partial  pressure 
region  Ifr'^’-lO  '^atm  can  be  obtained  either  by  a CO/CO2  gas  mixture  or  by  a purified 
argon  gas.  Figure  5-19  demonstrates  the  conductivity  results  of  BajCe{PO^)^  in  the 

above  two  gases  in  the  same  oxygen  partial  pressure  region.  Conductivity  in  CO/CO2 
mixture  gas  is  lower  when  oxygen  partial  pressure  is  less  than  10''^atm.  It  suggests  that 
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either  CO/CO2  does  not  come  to  equilibrium  or  the  CO/CO2  gas  mixture  lowers  the 
conductivity  of  Ba^Ce{PO^)^  when  the  oxygen  partial  pressure  is  less  than  10'’^atm. 

5.4  Conclusion 

The  relative  density  of  Ba^Ce(PO^)^  increased  with  sintering  temperature,  from 
78%  at  1100  °C,  to  89%  at  1200  °C,  to  95%  at  1300  °C.  However,  impurity  phase 
Ba2P207  was  produced  with  a sintering  temperature  of  1300  °C  for  8 h or  with  a sintering 
temperature  of  1200  °C  for  24  h. 

The  conductivity  of  Ba^Ce{PO^)^  is  8.91x10'^  S/cm  at  600  °C  in  wet  oxygen, 
which  is  higher  than  undoped  LaP04  and  YPO4  but  lower  than  5 mol%  Sr  doped  LaP04. 
The  conductivities  of  Ba^Ce(PO^)^  at  600  °C  in  dry  and  wet  4vol%  hydrogen  are 

4.04x10'^  and  1.46x10'^  S/cm,  respectively.  This  is  the  first  time  that  proton 
conductivities  were  studied  in  the  low  oxygen  partial  pressure  region,  especially  in  a 
hydrogen  atmosphere.  However,  the  total  conductivity  of  Ba^Ce(PO^)^  is  not  as  good  as 
that  of  perovskite  in  the  measured  temperature  region,  which  is  usually  10'^- 10'^  S/cm  at 
600  °C. 

The  Ba^Ce(PO^)^  conductivity  dependence  of  oxygen  partial  pressure  and  water 
partial  pressure  data  supported  the  proposed  defect  structure. 

In  the  future,  Ba^Ce(PO^)^  will  be  doped  at  Ce  site  with  divalent  elements  and  its 
conductivity  behavior  will  be  studied.  Higher  conductivity  will  be  expected  in  the  doped 
Ba,Ce{PO,\. 
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sintering  tenperature  (°C) 


Figure  5-1.  Relative  density  with  different  sintering  temperatures  for  8h. 


Figure  5-2.  XRD  of  Ba^Ce{PO^)^  with  different  sintering  temperatures  and  dwell  times. 
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Figure  5-3.  SEM  images  of  fracture  cross  sections  of  Ba^Ce{PO^)^  pellets  with  a dwell 

time  of  8 h and  different  sintering  temperatures,  (a)  1 100  "C,  (b)  1200  “C  and 
(c)  1300  °C. 
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Figure  5-4.  Porous  Pt  electrode  on  the  surface  of  a Ba^Ce(PO^)^  pellet. 
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Figure  5-5.  Impedance  plot  of  Ba^Ce(PO^)^  in  different  atmospheres,  (a-d)  is  the 
impedance  plot  of  Ba^Ce{PO^)^  in  oxygen,  air,  nitrogen  and  4vol% 
hydrogen,  respectively,  (e-f)  is  the  enlarged  impedance  plot  of  (d). 
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Figure  5-6.  Arrhenius  Plot  of  Ba^Ce{PO^)^  under  different  atmospheres. 
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Figure  5-7.  Conductivity  comparison  of  Ba^Ce(PO^)^  with  reported  data  from  literature. 
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Figure  5-8.  Conductivity  of  Ba^Ce{PO^)^  as  a function  of  oxygen  partial  pressure  in  the 
whole  oxygen  partial  pressure  region:  (a)  at  600“C  and  (b)  at  550“C. 
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Figure  5-9.  Conductivity  of  Ba^Ce(PO^)^  as  a function  of  oxygen  partial  pressure  in  the 
case  1 of  low  oxygen  partial  pressure  region. 
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Figure  5-10.  Conductivity  of  Ba^Ce(PO^)^  as  a function  of  oxygen  partial  pressure  in 
the  case  2 of  low  oxygen  partial  pressure  region. 
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Figure  5-11.  Conductivity  of  Ba^Ce(PO^)^  as  a function  of  oxygen  partial  pressure  in 
the  case  4 of  low  oxygen  partial  pressure  region. 


Figure  5-12.  Conductivity  of  Ba^Ce(PO^)^  as  a function  of  oxygen  partial  pressure  in 
the  case  3 of  low-intermediate  oxygen  partial  pressure  region. 
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Figure  5-13.  Conductivity  of  Ba^Ce(PO^)^  as  a function  of  oxygen  partial  pressure  in 
the  case  4 of  low-intermediate  oxygen  partial  pressure  region. 
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Figure  5-14.  Conductivity  of  Ba^Ce(PO^)^  as  a function  of  oxygen  partial  pressure  in 
the  intermediate  oxygen  partial  pressure  region. 
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Figure  5-15.  Conductivity  of  Ba^Ce{PO^)^  as  a function  of  oxygen  partial  pressure  in 
the  high  oxygen  partial  pressure  region. 
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Figure  5-16.  Conductivity  of  Ba^Ce{PO^)^  as  a function  of  water  partial  pressure  in 
hydrogen. 
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Log(PH^O)  (atm) 

Figure  5-17.  Conductivity  of  Ba^Ce(PO^)^  as  a function  of  water  partial  pressure  in 
nitrogen. 
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Figure  5-18.  Conductivity  of  Ba^Ce(PO^)^  as  function  of  hydrogen  concentration  with 
/>/^^o=0.0312atm. 
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Figure  5-19.  Conductivity  of  Ba^Ce(PO^)^  under  the  CO/CO2  gas  mixture  compared 
with  that  under  purified  argon  gas. 
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Table  5-1.  Electrical  conductivity  data  of  Ba^Ce{PO^)^  under  different  atmospheres. 


oo 

Ea  (eV) 

<^600-0 

O550“C 

dry  air 

2.8747E2 

0.821 

5.8E-6 

3.24E-6 

wet  air 

3.9546E2 

0.810 

8.93E-6 

5.17E-6 

dry  nitrogen 

2.6224E2 

0.800 

7.05E-6 

3.99E-6 

wet  nitrogen 

4.1068E2 

0.829 

7.67E-6 

4.08E-6 

dry  oxygen 

2.6020E2 

0.781 

8.78E-6 

5.08E-6 

wet  oxygen 

2.4434E2 

0.776 

8.91E-6 

5.13E-6 

dry  4vol%  H2 

7.7983E3 

0.932 

4.04E-5 

1.83E-5 

wet  4vol%  H2 

7.191E3 

0.815 

1.46E-4 

8.69E-5 

Table  5-2.  Conductivity  dependence  of  oxygen  partial  pressure. 


Oxygen  partial  pressure  region 

Defect  concentration 
dependence 

Conductivity 

dependence 

Low  po^ 
region 

Low  p„^o 
No  free  H2 

4 

1 

p^l  at  550"C 

Low  ph^o 
With  free  H 2 

1 

Pq®  if  vacancy  dominates 
1 1 

Pol  Ph  0 proton  dominates 

1 

Initially  p^^  , then 
decreases  with 
decreasing  p^^ 

High  ph^o 
No  free  H2 

1 

if  vacancy  dominates 
1 1 

Po^Pho  if  proton  dominates 

Not  studied  due  to 
equipment  limitation 

High  pff^o 
With  free  //j 

1 I 

Pol  PH2O 

negligible 

Low- 

intermediate 
Po,  region 

Low  ph^o 
No  free  H2 

4 

Not  studied  due  to 
equipment  limitation 

High  ph^o 
No  free  H2 

1 

if  vacancy  dominates 

1 I 

Po^Pho  if  proton  dominates 

1 

p^^  at  500  and  550°C, 

however,  not  enough 
data 

High  ph^o 
With  free  H2 

1 

p^^^  if  vacancy  dominates 

1 I 

Pol  Ph-p  ii^ proton  dominates 

4 

Intermediate  region 

Independent 

Independent 

High  po^  region 

dependence  was  not  derived 
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Table  5-3.  Conductivity  dependency  of  water  partial  pressure. 


Po^  region 

Defect  concentration 
dependence 

Conductivity  dependence 

Low  po^ 
region 

Low  p„^o 
No  free  H2 

Not  derived  due  to  low  water  partial  pressure 

Low  p„^o 
With  free  H 2 

High  pff^o 
No  free  H2 

I 1 

P^Pko  if  oxygen 

vacancies  dominate; 
1 1 

Po]Ph,o  if  protons 
dominate 

Not  studied  due  to 
experimental  condition 

High  p„^o 
With  free  H 2 

1 1 

Pol  Phjo 

1 1 

(6OOOC),  pJl^  (550»C) 

Low- 

intermediate 

Low  ph^o 
No  free  H2 

Not  derived  due  to  low  water  partial  pressure 

Po^  legion 

High  pfj^o 
No  free  H2 

11  11 
PoIpJi^o  ~ Ph^Pn^o 

Not  studied  due  to 
experimental  condition 

High  ph^o 
With  free  H2 

11  11 
PqIPh^o  ~ Pn-^PHfi 

Intermediate  pg^  region 

Independent 

Independent 

High  po^  region 

Electroneutrality  condition  faces  a challenge, 
dependence  was  not  derived 

Table  5-4.  Conductivity  dependence  of  hydrogen  partial  pressure. 


Low  pq^  region 

Defect  concentration  dependence 

Conductivity 

dependence 

Low  p,^^Q , Without  free  H2 

Not  derived  due  to  no  free  hydrogen 

Low  pff^Q  , With  free  H 2 

1 1 

PhjPh^o  if  oxygen  vacancy 

Not  studied 

1 

dominates;  p^^  if  proton 
dominates 

High  pfj^o , Without  free  H2 

Not  derived  due  to  no  free  hydrogen 

High  pff^o , With  free  H2 

4 

1 

Ph, 

CHAPTER  6 

TRANSFERENCE  NUMBER  MEASUREMENT 

6.1  Introduction 

Chapter  5 presented  the  measured  conductivity  of  Ba^Ce(PO^)^ , which  includes 

the  contribution  from  protons,  electrons  and  oxygen  vacancies.  Information  about  proton 
transference  number  is  needed  to  assess  the  proton  conductivity  and  to  determine  if 
Ba^Ce{PO^  )a  is  a mixed  conductor  or  pure  proton  conductor.  If  it  is  a mixed  conductor, 

hydrogen  flux  will  be  measured  directly  in  single  phase  Ba^Ce{PO^  )3 . If  it  is  a pure 

proton  conductor,  it  will  be  mixed  with  Pd  or  Ce02  to  form  a composite  and  hydrogen 
flux  through  the  composite  will  be  measured. 

The  transference  number  of  Ba^Ce(PO^)^  in  hydrogen  atmosphere  was  studied  by 


employing  the  following  cells: 

II,  100%  H2,  ?t/Ba,CeiPO,),/?X,  4%  Hj,  I (cl) 

II,  50%  H2+balance  He,  Pt/ Ba^Ce{PO, /Pt,  4%  H2, 1 (c2) 

II,  25%  H2+balance  He,  Pt/ Ba,Ce{PO, ),  /Pt,  4%  H2, 1 (c3) 

II,  12.5%  H2+balance  He,  Pt/ Ba,Ce(PO,),  /Pt,  4%  H2, 1 (c4) 

where  I and  II  correspond  to  the  two  sides  of  the  cell. 


The  transference  number  of  Ba^Ce(PO^)^  in  wet  hydrogen  atmosphere  was 

studied  by  flowing  hydrogen  gases  on  both  sides  of  the  above  cells  bubbled  through 
water  at  room  temperature. 
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In  a hydrogen  partial  pressure  gradient,  the  theoretical  electromotive  force  can  be 
expressed  as 


If  Ba^Ce{PO^)^  is  a mixed  conductor,  the  experimental  value  of  the  electromotive 

force  will  be  smaller  than  the  theoretical  one. 

The  experimental  electromotive  force  in  a hydrogen  partial  pressure  gradient  can  be 
expressed  as 


The  experimental  electromotive  force  in  a hydrogen  and  water  partial  pressure 
gradient  can  be  expressed  as 


due  to  the  possibility  of  both  proton  and  oxygen  ion  conductivity. 

To  measure  the  transference  number  between  300  and  600  °C,  a gas  tight  seal  must 
be  formed  between  the  sample  and  the  supporting  alumina  tube.  Rubber  0-ring  seals  can 
only  be  used  below  400  °C  due  to  their  heat  resistance  limits.  Metal  ring  seals  are  used  at 
high  temperature  but  their  thermal  stress  during  heating  or  cooling  causes  sample 
cracking  [80].  Glass-based  seals  have  better  spreading  properties  and  thermal 
compatibility  with  the  ceramic  membranes  and  supporting  tube.  The  sealing  material 
must  meet  the  following  requirements  [81-82]:  (1)  it  should  have  good  wettability  with 
the  membrane  and  the  supporting  tube,  (2)  it  should  fit  and  stay  in  the  gap  between  the 
membrane  and  tube  at  the  testing  temperature  range,  (3)  it  should  have  no  reaction  with 


r 

^ theoretical 


(6-1) 


P — f P 

^ exp  erimental  * ion  ^ theoretical  * 


(6-2) 


(6-3) 
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the  membrane  and  tube,  (4)  the  thermal  expansion  coefficient  of  the  seal  should  be 
compatible  with  the  membrane.  Because  of  the  intermediate  temperature  range  and  large 
thermal  expansion  coefficient  (15xlO'^/°C)  of  Ba^Ce(PO^)^ , conventional  sealants  do 
not  work  well.  The  thermal  expansion  coefficient  of  SZP-1  glass  is  in  the  range  10- 
12xlO'^/°C  [83].  In  this  chapter,  it  was  used  to  seal  Ba^Ce(PO^)^  membrane  and  AI2O3 
tube. 

The  transference  numbers  obtained  from  the  above  four  cells  should  be  identical  at 
the  same  temperature  if  sealant  is  perfect.  If  sealant  is  not  good,  the  true  hydrogen 
concentration  on  the  outlet  of  gas  I (4vol%  H2)  will  be  detected  by  the  mass  spectrometer. 

6.2  Experiment 

The  experimental  setup  for  Ba^Ce{PO^)^  transference  number  measurement  is 

illustrated  in  Figure  6-1.  Platinum  wire  was  attached  to  the  sample  and  carefully 
connected  with  electrode  leads  inside  the  setup.  SZP-1  glass  powder  was  pasted  around 
the  sample  and  tube.  The  whole  setup  was  heated  to  450  °C  at  3 °C/min  and  held  for  half 
an  hour  to  seal,  followed  by  a slow  heating  to  the  testing  temperature.  Gas  flow  rate  was 
controlled  at  50  seem.  The  open  circuit  voltage  was  measured  with  a Keithley  2000 
multimeter.  The  true  hydrogen  concentration  on  the  outlet  of  gas  I was  monitored  by  the 
QIOOMS  Dycor  Quadlink  mass  spectrometer. 

6.3  Results  and  Discussion 

The  total  ionic  transference  number  of  SrCeo.gsYo.osOa-s  in  dry  hydrogen 
concentration  cells  was  considered  the  proton  transference  number  by  Guan  [13].  This 
assumes  no  oxygen  ionic  transport  and  is  the  approach  we  are  using. 
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In  the  proposed  defect  structure  model  of  Ba^Ce{PO^  )j , a dry  hydrogen 
atmosphere  is  considered  the  low  water  partial  pressure  condition  with  free  hydrogen. 

The  defect  concentration  is  derived  to  be  proportional  to  if  oxygen  vacancy 
dominates  and  it  is  derived  to  be  proportional  to  if  proton  concentration  dominates. 

The  conductivity  of  Ba^Ce(PO^)^  is  initially  proportional  to  and  the  conductivity 

dependence  decreases  at  lower  oxygen  partial  pressure,  demonstrating  that  the  proton 
concentration  becomes  dominant  when  oxygen  partial  pressure  decreases.  The  oxygen 
partial  pressure  in  the  above  dry  hydrogen  concentration  cells  is  around  10'^‘^~10'^^  atm; 
therefore,  proton  concentration  dominates  in  these  cells  consistent  with  the  assumption 

hon  ~^OHq  ■ 

In  the  proposed  defect  structure  model  of  Ba^Ce{PO^  )j , a wet  hydrogen 

atmosphere  is  considered  the  high  water  partial  pressure  condition  with  free  hydrogen. 
Both  free  hydrogen  and  water  vapor  make  contributions  to  proton  incorporation; 
therefore,  proton  concentration  dominates.  Conductivity  of  Ba^Ce(PO^)^  remains 

unchanged  when  oxygen  partial  pressure  is  below  10'^^  atm,  which  confirms  that  the 
proton  concentration  dominates  in  wet  hydrogen  atmosphere. 

When  proton  concentration  dominates,  the  ionic  transference  number  could  be 
considered  proton  transference  number  because  the  mobility  of  protons  is  much  higher 
than  that  of  oxygen  vacancies.  It  can  be  expressed  as 
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where  c is  the  concentration  of  protons,  c ■ is  the  concentration  of  electrons,  is 
the  mobility  of  protons  and  is  the  mobility  of  electrons.  The  ionic  transference 

number  is  not  only  related  to  the  concentration  of  charge  carriers  but  also  related  to  their 
mobilities. 

Figure  6-2  shows  the  open  circuit  voltage  of  Ba^Ce{PO^)^  as  a function  of  hydrogen 
partial  pressure  ratio  in  the  dry  and  wet  hydrogen  concentration  cells  at  600  and  550  “C. 
The  proton  transference  number  is  0.53  in  cell  4 under  dry  hydrogen  and  0.52  under  wet 
hydrogen  at  600  “C.  Figure  6-2,  however,  also  shows  that  the  proton  transference  number 
increased  with  decreasing  hydrogen  partial  pressure  ratio,  indicating  the  SZP-1  glass 
sealant  is  not  perfect.  In  order  to  eliminate  the  leakage  effect,  the  true  hydrogen 
concentration  on  the  outlet  of  gas  1 was  monitored  by  mass  spectrometer. 

Figure  6-3  is  the  calibration  curve  of  the  mass  spectrometer.  In  the  calibration, 
hydrogen  gases  with  different  concentrations  were  fed  into  the  mass  spectrometer  to 
achieve  the  corresponding  hydrogen  partial  pressures.  Hydrogen  partial  pressure  as  a 
function  of  concentration  was  plotted.  The  true  hydrogen  concentration  in  gas  I outlet 
during  the  experiment  can  be  obtained  from  its  partial  pressure  in  figure  6-3. 

The  true  hydrogen  concentrations  in  the  gas  I outlet  for  the  four  cells  are  listed  in 
table  6-1.  The  values  are  higher  than  the  applied  concentrations  and  increase  with 
increasing  hydrogen  partial  pressure  on  the  side  II,  indicating  either  hydrogen  permeation 
or  the  existence  of  a leakage.  Argon  gas  is  detected  on  side  II  by  the  mass  spectrometer 
when  it  is  fed  on  the  side  I,  indicating  the  existence  of  a leakage.  The  true  theoretical 
electromotive  forces  were  obtained  by  using  formula  (6-1)  and  the  true  hydrogen 
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concentration  ratio.  The  measured  voltages  in  the  hydrogen  cells  are  also  listed  in  table  6- 

1. 

Figure  6-4  shows  the  open  circuit  voltage  of  Ba^Ce(PO^)^  in  dry  hydrogen 

concentration  cells  at  600°C  after  concentration  correction.  With  mass  spectrometer 
correction,  the  proton  transference  number  showed  no  apparent  dependence  on  the 
hydrogen  partial  pressure  ratio.  For  cell  1,  the  proton  transference  number  is  0.77  in  dry 
hydrogen  at  600°C.  The  proton  transference  number  studies  in  wet  hydrogen 
concentration  cells  are  more  complex  due  to  the  hydrogen  and  water  partial  pressure 
effect.  No  further  calibration  or  analysis  was  done  in  wet  hydrogen  concentration  cells. 

Hydrogen  concentration  cell  1 was  also  used  to  investigate  the  relationship  between 
proton  transference  numbers  and  temperature.  The  true  hydrogen  concentration  on  the 
gas  I outlet  in  cell  1 remains  15.4%  at  500-600°C.  The  corresponding  theoretical  and 
experimental  electromotive  forces  are  listed  in  table  6-2.  Figure  6-5  shows  the  proton 
transference  number  of  Ba^Ce(PO^)^  in  the  dry  hydrogen  atmosphere  as  a function  of 

temperature.  The  proton  transference  number  remains  0.79  from  500  °C  to  550  °C. 

The  proton  transference  number  should  have  no  relationship  with  hydrogen  partial 
pressure  ratio.  Cells  1-4  correspond  to  the  same  oxygen  partial  pressure  region. 

Therefore,  proton  transference  number  is  the  same  in  these  cells.  No  further  analysis  of 
proton  transference  as  a function  of  temperature  was  done  on  the  cells  2-4. 

The  conclusion  drawn  from  this  data  is  that  Ba^Ce(PO^)^  is  a mixed  conductor 
with  dominant  proton  conduction.  The  proton  conductivity  is  approximately  3.23x10'^ 


S/cm  at  600  °C  in  dry  hydrogen  atmosphere. 
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6.4  Conclusion 

Ba^Ce{PO^)^  is  a mixed  conductor  with  dominant  proton  conduction.  The  proton 
transference  number  and  electronic  transference  number  are  around  0.8  and  0.2, 
respectively,  at  600  °C  in  the  dry  hydrogen  concentration  cells.  Therefore,  the  proton 
conductivity  of  Ba^Ce(PO^)^  is  around  3.23x10'^  S/cm  at  600  °C  in  dry  hydrogen 

atmosphere  and  electron  conductivity  is  8.1x10’^  S/cm. 
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Figure  6-1.  Schematic  illustration  of  the  experimental  setup  for  transference  number 
measurement. 


Figure  6-2.  Open  circuit  voltage  of  Ba^Ce(PO^)^  as  a function  of  Log(PH2"/PH2')  in  the 
wet  and  dry  concentration  cells  at  600  and  550°C. 
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H2  voI% 

Figure  6-3.  Hydrogen  concentration  calibration  of  the  mass  spectrometer. 


Log(PH2"/PH2‘) 


Figure  6-4.  Open  circuit  voltage  of  Ba^Ce(PO^)^  as  a function  of  Log(PH2"/PH2‘)  with  or 
without  MS  correction  in  dry  concentration  cells  at  600°C. 
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Figure  6-5.  Proton  transference  number  as  a function  of  temperature  in  the  dry  hydrogen 
atmosphere. 
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Table  6-1.  Corrected  hydrogen  concentration  in  cell  l-cell4. 


cell 

applied  hydrogen 
concentration  of  gas  I 

true  hydrogen 
concentration  of  gas  I 

600T 

dry  H2  cell 

Emeasurement 

(mV) 

true 

Etheoretical 

(mV) 

1 

4vol% 

15.4vol% 

-54.30 

-70.36 

2 

4vol% 

9.4vol% 

-48.62 

-62.86 

3 

4vol% 

5.6vol% 

-39.28 

-56.27 

4 

4vol% 

4.5vol% 

-26.67 

-38.43 

Table  6-2.  Proton  transference  number  of  Ba3Ce(P04)3  at  different  temperatures  in  cell  1. 


Emeasurement  (mV) 

Etheoretical  (mV) 

transference  number 

500“C 

-49.10 

-62.31 

0.79 

550”C 

-52.08 

-66.34 

0.79 

600”C 

-54.30 

-70.36 

0.77 

CHAPTER  7 

HYDROGEN  PERMEATION  MEASUREMENT 


7.1  Introduction 


Hydrogen  permeation  can  be  defined  as  a process  by  which  pure  hydrogen  is 
selectively  passed  through  a separation  membrane  from  a gas  mixture.  The  driving  force 
usually  is  a pressure  or  potential  gradient.  The  hydrogen  permeation  process  is  a kinetic 
process.  The  characteristic  length  L(.=D*  !k^,  where  ks  is  the  surface  exchange 
coefficient  and  D*  is  the  self  diffusion  coefficient,  determines  the  transition  from  bulk 
diffusion  limited  to  surface  exchange  rate  limited  transport  [84].  When  the  thickness  L is 
larger  than  , the  hydrogen  flux  is  governed  by  the  bulk  diffusion.  A 2 pm  film  still  was 

considered  as  bulk  diffusion  controlled  in  Yb  doped  strontium  cerate  [85].  The  thickness 
of  all  the  membranes  in  our  test  was  larger  than  1 mm,  so  they  should  he  bulk  diffusion 
controlled. 

The  ionic  transference  number  of  Ba^Ce(PO^)^  in  hydrogen  concentration  cells 
indicates  that  it  is  a mixed  conductor  with  dominant  proton  conduction.  Therefore,  its 
hydrogen  permeation  flux  could  be  measured  directly  on  single  phase  Ba^Ce(PO^)^ . 

Generally,  the  hydrogen  permeation  flux  J across  an  oxide  membrane  can  be 

calculated  using  the  Wagner  equation,  which  assumes  the  bulk  diffusion  to  be  the  rate- 
limiting  step  [35]. 


(1-17) 
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According  to  literature  [13],  the  hydrogen  permeation  rate  of  0.1 1 cm  thick 
SrCeo,95Yo,o503-5  membrane  is  only  0.006  cm^  (STP)  cm'^  min’’  at  600°C  while  its 
conductivity  is  around  10'  ~10'  S/cm  and  proton  transference  number  is  around  0.86  in  a 
dry  hydrogen  concentration  cell  (4%  H2/  0.488%  H2). 

The  total  conductivities  of  Ba^Ce(PO^)^  are  4.04x10'^  and  1.46x10''^  S/cm  at 

600°C  in  dry  and  wet  4vol%  hydrogen,  respectively.  The  proton  transference  number  was 
around  0.8  and  the  electronic  transference  number  was  around  0.2  at  600°C  in  dry 
hydrogen  concentration  cells.  Therefore,  a low  hydrogen  flux  of  Ba^Ce(PO^)^  is 
expected  due  to  its  low  total  conductivity.  Hydrogen  permeation  is  only  studied  in  dry 
hydrogen  atmosphere  in  this  work. 

As  mentioned  before,  two-phase  composites  can  provide  continuous  channels  for 
both  proton  and  electronic  conductivity  and  enhance  hydrogen  transport.  Hydrogen 
permeation  was  tried  on  the  composites  BajCe(PO^)^  -Ce02  and  Ba^Ce{PO^)^  -Pd,  in 

which  Ce02  or  Pd  was  used  to  increase  the  electronic  conductivity  in  the  composites. 

7.2  Experiment 

For  Ba^Ce{PO^)^  membrane  samples,  the  sintering  procedures  were  the  same  as 
that  for  conductivity  measurement.  For  composite  Ba^Ce(PO^)^  -Ce02  or 
Ba^Ce(PO^)^  -Pd,  powders  were  mixed  and  ball  milled  for  24  h and  then  sintered  at 

1450  °C  for  10  h.  Their  relative  densities  are  higher  than  90%. 

The  hydrogen  permeability  measurement  system  is  shown  in  Figure  7-1.  The 
sealant  is  glass  powder  V-1003 1 . The  feed  gas  is  a mixture  of  dry  helium  and  dry 
hydrogen  with  a flow  rate  of  50  seem.  The  sweep  gas  is  helium  and  flow  rate  is  20  seem. 
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A small  amount  of  argon  gas  was  flowed  on  the  feed  side  to  check  for  a gas  leakage  in 
the  system.  The  gas  partial  pressure  on  the  sweep  side  was  measured  by  a mass 
spectrometer  (QIOOMS  Dycor  Quadlink). 

7,3  Results  and  Discussion 

7.3.1  Hydrogen  Permeability  of  Ba^Ce(PO^)^ 

Figure  7-2  shows  the  partial  pressure  of  various  gases  at  the  sweep  side  while 
testing  Ba^Ce(PO^)^  at  600  °C.  In  the  first  30mins,  only  helium  and  a small  amount  of 

argon  were  fed.  Therefore,  the  partial  pressures  of  other  gases  such  as  hydrogen, 
nitrogen,  and  oxygen  could  be  considered  as  baseline.  The  sealing  of  the  system  appeared 
to  be  sufficient  because  the  argon  partial  pressure  was  lower  than  8x10'”  kPa  and  the 
partial  pressure  of  nitrogen  and  oxygen  was  lower  than  2.67x1 0'^*^  kPa.  After  collecting 
baseline  for  30  mins,  helium  was  shut  down  and  pure  hydrogen  was  fed.  The  hydrogen 
partial  pressure  at  the  sweep  side  did  not  increase.  This  indicates  that  Ba^Ce{PO^)^  does 
not  have  good  hydrogen  permeability  under  these  conditions. 

7.3.2  Hydrogen  Permeability  of  Composite  Ba^Ce{PO^)^  with  Ce02 

Ba^  Ce{PO^  )3  is  a mixed  conductor  with  proton  conduction  dominance.  However, 
based  on  results  above,  its  electronic  conductivity  and  proton  conductivity  do  not  appear 
high  enough  for  it  to  have  good  hydrogen  permeability.  As  mentioned  in  chapter  3,  CeOj 
and  Ba^Ce{PO^)y  were  observed  in  composite  BCGSPd.  That  BCGSPd  sample 
demonstrated  a higher  hydrogen  flux  than  solid  state  synthesized  composite  BCG35Pd 
when  the  operating  temperature  was  lower  than  550  °C.  There  is  a possibility  that  CeOj 


may  have  some  contribution  to  that  hydrogen  flux.  Therefore,  Ba^Ce(PO^)^  was  mixed 
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with  CeOj  to  check  whether  or  not  electron  hopping  between  in  Ba^Ce{PO^)^  and 
Ce"^"^  in  CeOj  facilitated  the  hydrogen  permeation. 

Figure  7-3  shows  the  XRD  result  of  composite  CQO2- Ba^Ce{PO^)^  before  and 
after  sintering.  Ce02  and  Ba^Ce(PO^)^  were  indexed.  No  tertiary  phase  was  detected  by 

XRD.  Figure  7-4  shows  the  cross  section  of  composite.  There  is  no  interconnected  pore 
structure.  Figure  7-5  shows  the  partial  pressure  of  gases  at  the  sweep  side  while  testing 
the  composite  CqOi-  Ba^Ce{PO^)^  at  600  °C.  Helium  and  a small  amount  of  argon  were 
fed  in  the  first  30mins.  No  leakage  was  detected.  Then  helium  was  shut  off  and  hydrogen 
was  feed.  The  hydrogen  partial  pressure  on  the  sweep  side  is  1.80x10“'°  kPa,  slightly 
larger  compared  with  baseline  1.52x10“'°  kPa  in  figure  7-2.  The  actual  hydrogen  partial 
pressure  caused  by  permeation  is  2.8x10“"  kPa.  The  hydrogen  permeability  of  composite 
is  almost  zero  when  considering  the  mass  spectrometer  resolution.  This  indicated  that  the 
electronic  conductivity  provided  by  the  hopping  between  Ce^“^  and  Ce'*'^  was  not  enough 
for  hydrogen  permeation. 

7.3.3  Hydrogen  Permeability  of  Composite  35voi%Pd-65voI%  Ba^Ce(PO^)^ 

Ba-^Ce(PO^)^  was  also  mixed  with  35vol%Pd  and  sintered  to  form  a composite. 
Figure  7-6  shows  the  XRD  results  of  composite  35vol%Pd-65vol%5a3Ce(/'04)3 
sintered  in  air  and  treated  in  argon  atmosphere.  No  tertiary  phase  was  detected  by  XRD. 
Figure  7-7  shows  the  surface  image  of  composite  Vd- Ba^Ce(PO^)^ . Pd  isolated  itself 

from  Ba^Ce(PO^)^  grain  due  to  the  difference  bond  type  of  metal  and  ceramic.  Figure  7- 


8 shows  the  EDS  results  on  a Pd  grain  and  a Ba^Ce(PO^)^  grain.  They  all  contained  Pd, 
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but  the  concentration  of  Pd  in  Ba^Ce(PO^)^  grain  is  low.  Figure  7-9  shows  there  are  no 

interconnected  pores  in  the  cross-section  of  composite. 

Figure  7-10  shows  the  partial  pressure  of  gases  on  the  sweep  side  while  testing 
composite  35vol%Pd-65vol%5ai3Ce(F’04)3 . The  argon  concentration  increased  very 
quickly  once  it  was  fed,  indicating  the  sealing  was  poor.  The  inherent  properties  of  the 
sealant  or  microcracks  in  the  sample  may  be  the  source  of  the  leak.  The  sample  was 
broken  after  the  furnace  was  cooled  down  to  room  temperature,  so  the  origin  of  the 
leakage  is  not  clear.  Composite  5vol%Pd-95vol%5a3Ce(P04)3  was  also  tried  and  had  a 
similar  leakage  problem.  Therefore,  while  these  composites  may  have  a high  hydrogen 
permeation,  inability  to  obtain  sealed  crack  free  sample  prevented  further  measurements. 

7.4  Conclusion 

Glass  powder  V- 10031  was  found  to  be  a good  sealant  for  Ba^Ce{PO^)^  or 
composite  CeOj- Ba^Ce{PO^)^ . However,  Ba^Ce(PO^)^  did  not  exhibit  any  significant 

hydrogen  permeability  with  a He  sweep  gas.  The  electronic  conductivity  provided  by  the 
electron  hopping  between  Ce^"^  and  Ce'*'^  was  not  enough  for  hydrogen  permeation  of 
composite  CqOi- Ba^Ce{PO^)^  under  these  same  conditions. 

However,  the  transference  number  results  in  chapter  6 indicate  that  Ba^Ce{PO^)^ 

has  an  electronic  transference  number  of  0.2  in  the  dry  hydrogen  concentration  cells  and 
therefore  should  exhibit  hydrogen  permeation.  One  major  difference  between  the 
permeation  and  transference  measurements  was  the  composition  of  the  sweep  gas.  When 
the  hydrogen  permeation  is  measured  with  a He  sweep  gas,  the  corresponding  oxygen 
partial  pressure  region  on  sweep  side  belongs  to  the  intermediate  region;  therefore,  the 
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electronic  conductivity  is  lower.  If  the  hydrogen  permeation  is  measured  in  the  low 
oxygen  partial  pressure  region  with  a hydrogen  sweep  gas,  the  electronic  conductivity 
may  be  sufficient  for  hydrogen  permeation.  Further  permeation  work  of  Ba^Ce(PO^)^ 
with  a hydrogen  sweep  gas  is  necessary. 
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Figure  7-1.  Hydrogen  permeability  measurement  system. 
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Figure  7-2.  Gas  partial  pressure  on  sweep  side  while  testing  Ba^Ce{PO^)^ . 


20(°) 


Figure  7-3.  XRD  patterns  of  composite  Ce02-Ba^Ce{PO^)j^  before  and  after  sintering. 
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Figure  7-4.  Cross-section  image  of  composite  CeOj-Ba^CeiPO^)^ . 


Figure  7-5.  Gas  partial  pressure  on  sweeping  side  while  testing  composite  CeOj  - 
Ba^Ce{PO^)^ . 


156 


Figure  7-6.  XRD  pattern  of  composite  35  vol%  Pd-65  vol%  Ba^Ce(PO^)^  after  separate 
treatment  in  air  and  argon. 


Figure  7-7.  Surface  image  of  composite  65  vo\%  Ba^Ce(PO^);^  -35  vol%  Pd. 
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Figure  7-8.  EDS  results  of  Ba^Ce{PO^)^  grain  and  Pd  grain  in  the  composite  65  vol% 
5a3Ce(POj3-35  vol%Pd. 


Figure  7-9.  Cross-section  image  of  composite  35  vol%  Pd-65  vol%  Ba^Ce{PO^)^ . 
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Figure  7-10.  Gas  partial  pressure  on  sweeping  side  with  composite  35  vol%  Pd-65  vol% 
Ba^Ce(PO^)^ . 


CHAPTER  8 
CONCLUSION 

In  previous  work  in  our  lab,  the  hydrogen  flux  of  composite  BCG5Pd  was  higher 
than  that  of  BCG35Pd  when  the  temperature  was  lower  than  550°C.  Additionally,  this 
hydrogen  flux  did  not  achieve  an  equilibrium  value  for  20  hours  at  650°C  until  feed  gas 
was  shut  off  [36],  The  second  phase,  Pd,  did  not  exhibit  such  a trend  under  isothermal 
conditions,  indicating  the  existence  of  a tertiary  phase.  Analysis  results  showed  that  the 
phosphorous  impurity  was  introduced.  In  electroless  deposition,  exposure  to  acidic  bath 
decomposed  part  of  BCGO  to  form  Ce02  and  the  absorbed  phosphorus  impurity  covered 
the  BCGO  particles.  They  reacted  at  higher  temperature  to  form  a BaO-rich  compound 
Ba^Ce(PO^)^  and  leave  a certain  Ce02  excess.  Ba^Ce(PO^)^  is  most  likely  the  desired 

tertiary  phase  in  the  composite  BCG5Pd. 

Single-phase  Ba^Ce{PO^)^  was  synthesized  successfully  by  a co-precipitation 
method  at  a lower  temperature  than  that  required  for  solid-state  reaction.  It  was  stable  in 
dry  air,  helium,  argon,  and  hydrogen.  The  conductivity  of  Ba^Ce(PO^)^  is  8.91x10'^ 

S/cm  at  600  °C  in  wet  oxygen,  which  is  higher  than  undoped  LaP04  and  YPO4  but  lower 
than  5 mol%  Sr  doped  LaP04.  The  conductivities  of  Ba^Ce(PO^)^  are  4.04x10'^  and 

1.46x10'^  S/cm  at  600  °C  in  dry  and  wet  4vol%  hydrogen,  respectively.  The  proton 
transference  number  was  around  0.8  and  the  electronic  transference  number  was  around 
0.2  at  600°C  in  hydrogen  concentration  cells.  Therefore,  Ba^Ce(PO^)^  was  a mixed 
conductor  with  dominant  proton  conduction.  However,  neither  its  proton  nor  electron 
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conductivity  was  high  enough  for  hydrogen  separation  with  a He  sweep  gas.  The 
electronic  conductivity  provided  by  the  electron  hopping  between  Ce^"^  and  Ce'*'^  was  also 
too  low  to  allow  hydrogen  separation  in  the  composite  CcOi- B a under  these 

conditions.  It  is  still  possible  to  obtain  a hydrogen  flux  with  hydrogen  on  both  sides. 

The  results  of  Raman  spectroscopy  and  magic-angle  spinning  nuclear  magnetic 
resonance  showed  the  local  structures  of  Ba^Ce(PO^)^  at  different  atmospheres. 

Protons.  Free  hydrogen  can  be  oxidized  forming  protons  and  incorporated  into  the  lattice. 
In  the  presence  of  water  vapor,  bound  hydrogen  is  produced  and  oxidized  forming 
protons.  The  formed  protons  attach  onto  the  O atoms  of  PO4  group  and  are  incorporated 
into  the  lattice. 

The  defect  structure  of  Ba^Ce(PO^)^  was  proposed  in  the  low,  low-intermediate, 
intermediate  and  high  oxygen  partial  pressure  regions.  The  conductivity  of 
Ba^CeiPO^  )3  as  a function  of  oxygen  partial  pressure  and  water  partial  pressure 
supported  the  proposed  defect  structure. 
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